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SECTION 1 


1.0 INTRODUCTION 

This document is the Users Manual for the Solid Rocket Booster Perfor- 
mance Evaluation Model (SRB-II). Descriptions of the model, the program 
options, the required program inputs, the program output format and the 
program error messages are presented herein. This Users Manual is con- 
tained in Volume II of a set of four volumes vvhich describe SRB-II. 

Volume I contains the engineering description; Volume III contains the 
sample case; and Volume IV contains the program listing. 

SRB-II is written in FORTRAN Language and is operational on both the 
IBM 370/155 and the MSFC UNIVAC 1108 computers. The model has been 
extensively overlayed in an effort to meet the MSFC UNIVAC 1108 System 
core storage restriction of 32,000 words at any given time. 

Both analytical and functional simulation techniques are employed in 
SRB-II to predict performance of the propulsion elements in' the Space 
Shuttle Solid Rocket Booster. The model also contains analytical methods 
for reconstructing performance of motors for both static test and flight 
using measured data as direct input, SRB-II calculates the oerformance 
of a wide variety of possible solid rocket booster configurations using 
input data derived from theoretical analysis, small scale tests, large 
scale tests, and flight motor configurations. Three analytical scaling 
techniques are included in the model for the purpose of predicting 
delivered characteristic velocity (c*) and vacuum specific impulse . 
(Isp) of the motor. The predicted Isp and c* are used in conjunction 
v;ith motor propellant data to predict booster propulsion performance. 
Motor performance parameters such as thrust, flowrate and operating 
pressures may be calculated with an internal ballistics technique 
or scaled from available test data. 

The SRB-II Source and Executable Tape has been placed on file in the 
MSFC Data Systems Laboratory Library as tape number 29312. This tape 
contains the operational version of SRB-II demonstrated on the MSFC 
UNIVAC 1108 system, and the symbolic, relocatable, and absolute modules 
of SRB-II. 


T-1 
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SECTION 2 


2.0 PROGRAM DESCRIPTION 

The Solid Rocket Booster Performance Evaluation Model is made up of 
analytical and functional techniques combined to simulate the pro- 
pulsion elements of a solid rocket booster. The model v/as developed 

by making maximum use of existing computer programs and simulation 
techniques. Three existing comouter programs incorporated into SRB-II 
are the Boeing Internal Ballistics Computer Program, the Lewis Thermo- 
chemical Computer Program, and the One-Dimensional Two-Phase Flow Loss 
Computer Program. 

SRB-II has user options to select simulations for flight and static 
test prediction or reconstruction of a wide variety of solid rocket 
motor configurations. The options selected are dependent on data 
available to the user; therefore SRB-II may be input with data derived 
from theoretical analysis, small scale test data, large scale test 
data, and flight motor configuration data. SRB-II input is described 
in Section 3.0. 

The user has several output options for the calculated oropulsion 
system performance parameters, Output includes major parameters (isp 
thrust, flowrate, mass accounting, and onerating pressures) as a ’ 
function of time, as well as single point scaled data. SRB-II output 
is described in Section 4.0. 

The remaining paragraphs of Section 2.0 include a description of 
each SRB-TI Module, a description of the groundrules used in naming 
program variables and constants, and a discussion of the overlay which 
was required to attemnt to meet the MSFC UMIVAC 1108 System storage 
restriction of 32,000 words at any given time. 

2.1 MODULE DESCRIPTION 

The following procedure was used to synthesize the SRB-II in modular 
form. 

(1 ) The best solid rocket motor simulation techniques available were 
selected utilizing existing simulations where possible. 

(2) Computer routines were written to describe functions of individual 
components not already programmed. 

(3) The simulations and routines were grouped into modules. 

(4) The modules v/ere linked together with an Executive Module such 
that the path through th- model can be selected or controlled by 
the user. 

The modularized SRB-II program flowchart is shown in Figure 2-1 . All 
the modules are linked together with an Executive Module which is 
used to control the overall progfam for each simulation case. 
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2.1 (Continued) 

The following paragraphs will discuss each module in detail. 

2.1.1 Executive Module (MAIN) 


access to the other modules such that 

Vi 

various modules is controlled in the Executive Module by control 

Indicators are variables whose values a^fSelermined 

For^Sxamnip’"^fh“^^’" through the Centralized Incut Module 

cTiTiiii L user could either utilize or bvnass the BATES Isn 
?JSiLto?' appropriate value for the bX^E^Isp Module 

2.1.2 Centralized Input Module (INPUT) 

Ihp Module processes input data to the SRB-II usino 

the following four methods; BLOCK DATA NAMELIST qpt •Fnv.m=+ 

lu InpTAodXsS^o M„:° ; T^^cssses 

d k transfers the card input dat^^f 

reauire^the directly from.disk by the routines which 

^uir no thP directly by the routine re- 

reouirprl^hfra^^^' access by the program routine is 

quired because of core storage limitations set by the MSEC UNTVAC 

sI?tionT^^'^^ ^ complete discussion of input data is given in 


2.1.3 BATES Isp Module (BATES) 

Itn and characteristic 

velocity (c ), using the Air Force devised BATES method and is the 
preferred one of three modules provided for these calculations The 
Isp Scaling and Contractor Data Modules are considered alternate 
methods and are discussed in separate paragrachs. 

The Air Force devised the BATES method of accurately measuring Iso 

? hlilli Pnecision .test motorind closely 

fo^Sll1stifLs^cihT"^*’*®^®'■?■‘^®^ ' is an acronym 

I IS tic Test ^valuation ^calinp, Extensive use and modi first inn 

Drpdi^tfnn^^^ twelve years have extended this system 'to allow precise 
prediction of propellant performance in the expected "end item" or 
final production rocket motor. McuLeu ena isem or 

Procedure employs the use of a careful! v designed 
motor. Propellant samples are prepared' bv the 
ocket motor contractor in separate cartridges which are subsequentlj 


2-3 
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2.1.3 (Continued) 

loaded into a test motor much the same as a cartridqe is loaded into 
a gun. The motor is placed in a calorimeter capable of precise heat 
loss measurement and fired. An extremelv accurate test stand is 
manned by Air Force personnel who actually conduct the test, Air 
Force personnel also accomplish the data reduction, thus insurinq the 
consistency of all test procedures. Over 400 firinqs have been made 
to accurately determine the specific impulse of a number of propellants 
of interest to the Air Force. 


The BATES philosophy involves correcting an emnirically derived Isp 
by means of analytical corrections. The corrections required are 
described in the following paragraphs , 

2. 1.3.1 BATES No Loss Delivered Specific Impulse 

BATES "ideal" or no-loss delivered Isp is obtained by adding three major 
Isp correction factors to the PATES delivered Isp. These are heat 
flow, two-phase flow, and divergence losses. 

Heat flow loss is measured directly usinn the test stand calorimeter. 
These heat flow loss effects, as measured in the BATES program, can 
result in a loss of up to 5-1/2 seconds of Isp. 

A typical space shuttle solid rocket booster propellant contains 
sixteen percent aluminum resulting in an exhaust flow containing 
approximately thirty-percent solid particulate aluminum oxide (AlpDg). 
piese solid particles lag the gaseous flov/ and cause a performance 
decrement known as tv<o-phase flov' loss. The primary parameters 
required to calculate this flow loss are size and flov/ rate of the 
particles. A correlation for particle size and throat diameter has 
been adapted from the literature and programmed. This correlation is 
shown in Figure 2-2 and has been incorporated into the BATES Module. 
Particle flow rate will he calculated usino output from the Lewis 
Thermochemical Program (Reference 3) v/hich is incorporated as a sub- 
routine in the RATES Isp Module. Once the particle size and flow 
rate are known, the IDN/^ZL nrooram developed by the Air Force 
(Reference 4) is used to calculate the Isp decrement due to two-nhase 
flow loss. This program has been incorporated into the BATES Isp 
Module as Subroutine IND0ZL. 

The final Isp correction to the BATES small motor data, the divergence 
loss, accounts for the fact that all the exhaust velocities are not 
parallel to the axis of the nozzle. For a simple conical nozzle this 
loss is given by: 

IsPdiv ” - T/2 (cos ct) j Isn 

where Iso - theoretical specific impulse 

and a = nozzle half angle. 
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2. 1.3.1 (Continued) 

This divergence loss calculation has been incoroorated in the BATES 
Module, 


The three losses, when combined with the measured BATES Isp, result 
in the BATES no-loss delivered specific impulse. The BATES module 
coinputes a correlation with the theoretical Isp (Lewis Prograin) for 
use as a back-up. This provision allows performance extrapolation of 
similar propellants for which BATES data is not available by providing 
an estimated no-loss Isp. 

2. 1.3. 2 "End Item" Delivered Specific Impulse and Characteristic 
Velocity 


The basic assumption of the BATES method is that the delivered no-loss 
Isp of the small test motor is equal to the delivered no-loss Isp 
of a large motor for a given propellant composition. Therefore, the 
"end item" delivered Isp is obtained by deducting the large motor isp 
losses from the BATES no-Toss delivered Isp and then making any^necessary 
expansion ratio corrections. The major large motor Tosses are tv/o-phase 
flovj Toss, heat loss, divergence Toss and submergence Toss, ihe 4 . 1 , 
two Dhase flow and heat loss are accounted for in the same manner as the ^ 
previously described BATES motor losses. Submergence loss is also a required 
input. When the large motor has a conical nozzle the divergence losses 
are calculated as described for the BATES motor. ^ If the larp motor has 
a* contoured nozzle the divergence loss can be calculated with either of 
the following methods. The apDroximate method uses the following equation. 


where ISP = theoretical specific impulse 

a = equivalent conical angle of the contoured nozzle 
and 9 = angle of nozzle contour at the nozzle exit plane measured 

from the nozzle centerline 


The method of characteristics method uses the following eouati on 
ISP 


[1 - 1 (cos»] [L] [ISP] 


div ’’ ‘-2 2 

where ISP = theoretical specific impusTe 

rt = equivalent conical angle of the contoured nozzle 

and L = ratio fo the theoretical specific impulse of _a . 

nozzle to the theoretical specific impulse of an equivalent 
conical nozzle derived v/ith the ATPAP computer program. 

"End Item" motor characteristic exhaust velocity (C'*) is caTcu^ 
hv correctinq the theoretical c* at the nozzle throat. The corpction 
hL the form^of the ratio of the End Item motor actual ISP at the nozzle 
throat to the "end item" motor theoretical Isp at the nozzle throat. 
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2.1.4 Isp Scaling Module (SISCAL) 

T^n InH Module calculates “end item" Isp by calculatinq theoretical 

=l?ernlto ta’Sh’e'lATEl‘';Bchn1que!'''' ^ 

The previously mentioned Lewis computer program is incorporated into thp 
Isp Scaling Module as the means of calculatinq theoretiSi ^Jp^ 

is‘"cJrrectedYor^d^^^^^^ velocity. The theoretical imLlse 

Mn-f-nS T for divergence losses as previously described in 2 1 3 2 
inefficiencies will be taken into account by using the NASA* 
data {from Reference 5) shown as Figure 2-3 . This fiaurrshnS fhP 1 

types of prope?^anr^anrindic^^ various size motors with different 

spulation in the Isp Scaling Module. The derivation of delivered 
specific impulse used in the Isp Scaling Module is summarizedbelow; 

IsPdei = (A) (n^) (n^) (n^) (IsRj.^^^) 

where a = nozzle divergence loss factor 

Hp ~ nozzle efficiency factor 

rij, = c* efficiency factor 

Og = nozzle submergence loss factor 

^^Pdei “ delivered specific impulse 
^^*^theo ~ 'theoretical specific impulse 

"End item" motor characteristics exhaust velocity fc*l is ralniiatoH hw 
correcting the_ "end item" motor nozzle throat th^reti cal c* witli^the^^ 
motor c* efficiency Pc and the nozzle submergencrStor ns- 

2.1.5 Contractor Data Isp Module (CDS!) 

The Contractor Data Isp Module is used to calculate Iso whenever 

1. The Lewis program is used to find the theoretical snecifir 
for two different conditions: 

(a) The contractor test conditions (Isp 1) 
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2.1.5 (Continued) 

2. The contractor test motor divergence loss (Dl) and the "end 
item" divergence loss (02) are calculated. 


3. The data from Figure 2-3 are used to find the Isp efficiency 
for the contractor test motor mass flow rate (ETAl) and "end 
item" mass flow rate (ETA2). When the "end item" nozzle is 
submerged, ETA2 is corrected for a nozzle submergence factor. 

4. The "end item" vacuum Isp is predicted as follows: 

\ 


'=Pvac ' ’=Ptest • 


Isp 2 Isp 2 - 02 ETA2 
Isp 1 • Isp 1 - 01 • ETAl 


The "end item" motor c* is calculated by correcting the contractor test 
motor theoretical c* by the ratio of the Isp efficiency (ETA2/ETAT). 


2.1,6 Internal Ballistics Module (IBM) 


The Internal Ballistic Module is a major portion of SRB-II. IBM is used 
to calculate significant SRB propulsion performance parameters. The 
basis for this module is the Boeing Internal Ballistics Program 
(Reference 6). 

Jhe Internal Ballistics Module will: 

(a) ■ Generate nominal flow rate, chamber pressure, and thrust 

histories as a function of time. 

(b) Evaluate propeTlant mass and volume. 

(c) Utilize actual grain dimensions. 

(d) Evaluate burning rate scale-up effects. 

(e) Evaluate the changes in nozzle throat area as a function 
of time. 

(f) Utilize nozzle efficiency factors as a function 
of time. 

(g) Generate moments of inertia and center of gravity histories of 
solid propellant grains. 

(h) Match burning rate with measured head-end pressure for buildup 
prediction and performance reconstruction. 


2-9 
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2.1.6 (Continued) ' '■ 

= l Module can be used in conjunction with the 

as rfunSinn sffects Of "inert" mass exoenditures 

with thp TnSL ■ I'^. '^^snersion Module can ‘be used in conjunction 

he Internal Ballisi,rc.s ilodule to qenerate performance dispersions 

selected performance variables. IBM 
IS used in conjunction with the Reconstruction Module to analyze 
performance data from flight and static test. ncuyze 

The module calculations are based on either of the motor confiaurations 
shown^ln Figure 3-13 This figure shows examnles of the aeometriLl 
models ' or a monolothic, tanered grain motor and a seomented o^'ain 
motor. The motor configuration is divided into three sections: the 

head-end section (fonvard dome), the cylindrical section, and the aft- 
nead section (aft-dome). The grain oeometrv is described with input 
reference planes within the cylindrical section. The cylindrical 
section IS further divided into a number of increments or mass addition 
regions by the location of increment dividing planes. These dividing 
planes are located at each reference plane, and at specified intervals 
(A Z) from each reference plane, until either a segment slot interface 
I'eference plane is oassed. During the computer solution 
of the internal ballistics, the port perimeter, nort cross sectional 
area, and moments of inertia are determined at each increment dividing 
linear interpolation between adjacent reference planes. Mass 
addition is assumed to occur as a step orocess between two successive 

increment dividing Planes. 

In general, the major inputs required by the Internal Ballistics Module 
are dimensions and geometry factors to describe the motor case , nozzle and 
propel lantgrain, and the propellant characteristics. Propellant charac- 
teristics include density, specific heat ratio, gas constant, molecular 
weight, characteristic exhaust velocity, and coefficients and exponents 
for the burn rate equation which is an expanded form of 

R|3 =• a(Pp". 


where 

a 

n 

and P 


propellant burn rate 
propellant burn rate coefficient 
propellant burn rate exponent 
motor chamber pressure 


MO 
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2.1.6 (Continued) 


The module is concerned with the burning characteristics and orogression 
of the propellant surface, and the dynamics of the aas generation and 
qas flow in the interior of solid propellant rocket motors. The solution 
of the internal ballistics requires the geometric analysis of the grain 
design to determine the burning surface area and Dropellant volume, 
and the physical analysis of the gas flow through the port cavity and 
nozzle to determine motor nerformance. 

The Internal Ballistics Module can be used to analyze a large variety 
of propellant grain configurations. Also, numerous simulation options 
are available. Figure 2-4 presents some of the more important caoabiliti 
of the Internal Ballistics Module. 

2.1.7 Thrust Scaling Module (FSCAL) 

The Thrust Scaling Module provides an alternate performance prediction 
technique which can be utilized when sufficient data are not available 
for running the more sophisticated Internal Ballistics Module, or 
when a quick, cost effective, relatively simole technique of determining 
motor performance is desired. The user of this technique must realize 
that it is only an approximate method of deriving the parameters and 
is only valid for small perturbations about the nominal. The' 

Thrust Scaling Module will scale a nominal performance versus time 
curve to compensate for variations in burn time and propellant mass 
for a particular solid propellant motor. 

2.1.8 "Inert" Hass Module (INERT) 

The "Inert" Hass Module will predict "inert" mass flow rate and "inert" 
mass thrust contribution. "Inert" mass accounting is required because 
the total mass expended by a solid rocket motor exceeds the propellant 
mass. The excess is considered expended "inert" mass and is made uo 
primarily of insulation and ablative material contained in the motor 
interior and nozzle throat. This expended mass has a slight influence 
on thrust and a significant influence on motor mass versus time. SRB-II 
contains an "Inert" Mass Module which scales the "inert" mass remaining 
onboard by proportioning "inert" mass to propellant mass remaining. 
"Inert" mass flowrate is calculated from the mass overboard schedule, 
and then "inert" thrust contribution is calculated from the "inert" 
mass flowrate and the input "inert" Isp. "Inert" effects are added 
to appropriate parameters from either IBM or FSCAL before output in 
the Output Module. 
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2.1.8 (Continued) 

Required inputs are SII, AMIRR(I), AMPRR(I), NP(IO), MIT0T, AND MPT0T 
(If NF=0). Data for these variables is provided through BL0CK DATA 
or NAMELIST INPUTl of the C0NTR0L DATA PACKAGE. For definitions, see 
Table 3-XIII. 


2.1.9 Flight and Static Test Reconstruction Module (REC0N) 

The Reconstruction Module will accurately analyze the data resulting 
from both static tests and flight tests of the* Soace Shuttle SRB 
motors. Reconstruction of both flight and static test data is 
based on the technique which forces the prediction routine to produce 
a match between calculated and measured values of certain critical 
parameters. 

2.1.9.1 static Test Reconstruction 

Reconstruction of the SRB motor static tests is accomplished through the 
use of a progressing time steo algorithm which requires the calculated 
burning area to be consistent with the measured chamber pressure 
using the Internal Ballistics Module to force the match at each time 
step. As a result of this matching process, a history of motor mass 
flow rate and mass overboard is produced. 

2-. 1.9. 2 Flight Reconstruction 

The flight reconstruction is accomplished through the use of a 
progressing time step algorithm v/hich requires the calculated burning 
area to be consistent with measured chamber pressure using the Internal 
Ballistics Module to force the match. As a result of this matchina 
process, a history of motor mass flowrate and mass overboard is produced. 

Reconstruction of flight parameters using acceleration as the basis has 
been established as the most accurate method for evaluating solid booster 
flight performance. For the Shuttle, acceleration matching is complicated 
by simultaneous burning of the two solid boost motors and the three liquid 
main engines. Reconstruction is further complicated due to the aerodynamic 
forces generated by the complex shuttle configuration. In order to 
establish the contribution of the solid booster motors to the instantaneous 
vehicle acceleration, and thus establish the SRB total thrust, it is 
necessary to apply Newton's second lav/ to the vehicle in some specified 
direction. The resulting equation is then solved for total solid motor 
thrust' and may be expressed as: 

'^SRH ^ ^^^VEHICLE ^MEASURED^ “ ’^VEHICLE 
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2.1 .9.2 (Continued) 

where = The sum of thrust produced fay solid motors 

in the specified direction 

^^VeHICLE = Instantaneous total mass of the vehicle 

^MEASURED = f1easured acceleration in the SDecified 

direction 

^VEHICLE “ Summation of all forces acting on the 

vehicle in the specified direction 
excluding 

Instantaneous total mass, MveHICLEj is calculated from an input vehicle 
mass history which v/ill be the summation of the masses of all vehicle 
components, excluding the solid rocket motors, summed with the results 
of a mass balance performed on the solid rocket motors. This SRB mass 
balance is accomplished by algebraically summing the initial solid 
motor masses with the intenral of the existing mass flow rates cal- 
culated in the Internal Ballistics Module. Parts of the FveHICLE 
term are required input items and include liquid engine thrust, 
aerodynamic forces, an estimated thrust for the motor not being recon- 
structed in the run and any other forces acting on the vehicle in the 
specified direction. ‘ 

2.-1, 10 Dispersion Module (DISP) 

The Dispersion Module perturbs independent casual variables such 
that the effects of these perturbations are seen in the critical output 
quantities such as thrust, Isp, and mass flov/rate of the solid rocket 
motors. Perturbations are made on single IBM input variables such as 
nozzle erosion rate, propellant density or characteristic velocity 
while all other variables are held at their nominal values. Thus, a 
complete set of dispersion predictions can be generated by utilizing SRB-IJ, 

2.1.11 Output Module (0UTPUT) 

The Output Module is caoable of producing data in the follov/ing three 
forms: tabulated computer printout, seven and nine track-performance 

data tapes (IINIVAC 1108 version generates seven track tapes v/hile 
IBM-370 version generates nine track tapes) and plot tapes, and 
punched data cards. The last two items are optional outputs and 
their generation are determined by setting input indicators. The 
following paragraphs describe the three output methods available in 
this Module. 

2.1.11.1 Tabulated Computer Printout 

The tabulated computer printout, controlled by the Output Module, will 
provide the user v/ith tabulated data of performance parameters. The 
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2,1.11.1 ^(Continued) 

format of the printout displays key performance parameters. Variables of 
special interest may be obtained at the user s discretion by setting the 
appropriate indicators in the input data packages. 

A two part tabulated printout will be made for all computer '"uns. The 
first part will consist of both the input data entered through the 
Centralized Input Module and the results of all initial calculapons. 

The second part of the printout will contain performance parametp and 

variables for the entire simulation. The printout time increment 
will vary for the ignition transient, web time, and thrust decay __ 
portions of the simulation and will be determined by an input NAMELIST 
indicator. Typical output for 5RB-II is given in Section 4.0. 

2.1.11.2 Performance Data Tapes 

The user has the option to set an indicator which will require the model 
to deneratf a data tape. If the indicator is set, the^program automatically 
records the desired data on magnetic tape at predetermined time increments 
for the ignition transient, web time and thrust decay portions of the 
simulation. 

2.1.11 .3 Plot Tapes 

■plot tanes created by the Output Module are an intermediate step in 
obtaiSrg hfrd copy "plots. 

into a plot tape using existing Benson and Lerner (B/L) or MSFC SC 40^u 
Plot Programs. The hard copy plots can then be plotted on the B/L 
or SC 4020 plotter, 

2.1.11.4 Punched Data Cards 

The third Output Module option available is punched data cards, . 

nunched data cards are output when an indicator is set in the input data. 
Contents and format of the punched card output variables are discussed 

in Section 4.0. 

2.1.12 Mozzle Submergence and Contour Effects Module (NSCE) 

The Nozzle submergence and Contour Effects Hodul^ Si istics^an^^^^ 
of nozzle contour and submergence on the intern^ l^allistics and overai^ 
performance of a solid rocket motor. Contour effects are simulated^ 

Ssing functional modeling techniques which account for changes^ 

formance due to nozzle expansion section shape. Effects of submerging 
thrnSSlf entrance section in the motor grain are modeled using both. 
Analytical and functional simu 1 ation_techniques_to^accou^ 
in motor internal flow field and performance. Results from 

and submergence simulations are transmitted to As 

Module where computations 'fo'" critical output parameter ^ ^ 

specific ^mpulse^ flowrate, mass accounting and center of gravi^^^ 
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2.2 VARIABLE AND CONSTANT NAMES 

vvt^rti aeve I oped prior to Tncorporatron in the model. 

2.2.1 Common Storage 

isVhp’"LcLn”T°'^ wTlT be CWMON/BLKnnn/, where nnn 

steS "‘‘f and crogmstSo fn 

ComSon rard.' variable names are out on the'original 

2.2.2 FORTRAN Statement Numbers 

Subroutines have sequential statement numbers in increments of loj i .e. 

10 IF (X.LT.A) R0 T0 20 
B = C 
G0 Tf5 30 
2Q B = D 
C = F 

30 E = B * C 

corrections are to be inserted, statement numbers "sol it fhp 


2.2.3 


FORTRAN Nomenclature 




SYMBOL 

A 

C 

D 


DEFINITION 

altitude, (the last 

letter of A indicates the predicted value) 

coefficients, constants, specific head, G*, Cn CR 

. , F’ , ■ 

diameter, density, delta, distance or length 
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2.2.3 

SYMBOL 

E 

F 

G 

H 

I 

J 


(Continued) 


DEFINITION 


efficiency, expansion ratio 
thrust, force, fuel 
gamma 
height 

fixed point indicator 
fixed point indicator 



K fixed point indicator 

L fixed point indicator 

M M V'/ill be made REAL,(i.e., floating ooint), mixture 

ratio (MR), mass, mass flowrate (MF), moment, mach 
number 


N 

0 

P 

Q 

R 

S 

T 

V 
W 
X 

Y 
Z 


fixed point indicator 
oxidizer 

pressure, percent, propellant 
heat transfer 

radial distance, burn rate (RB) 
specific impulse (SI), stage 
temperature, time, total impulse 
volume, velocity 
weight 

the last letter of X means in the X direction 

the last letter of Y means in the Y direction 

the last letter of Z means in the Z direction or the 

previous value of a parameter 
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2.2.4 FORTRAN Constants 

The constants used in SRB-II have been assigned names according to the 
followina list. 


AK(I) 

CK(I} 

DK{I) 

IK(I) 

MC(I} 

MK(r) 

flK(I) 

PK(I) 

TK(I) 


curve fit: constants 
universal constants 
distance or dimension constants 
fixed point indicators 

motor constants and tag values (floating point) 
mass constants (floatino point) 
odd constants 

propellant characteristic constants 
time constants 


2.2.5 Program Indicator Concepts 

Indicators which are used specifically for one module will begin with 
I followed by a three letter abbreviation of thfi module and a counter 
number, i.e. an indicator for the BATES Module miaht be 

IBATOl 


Indicators which are used generally or in the Executive Module will 
with N and have a descriptive name such as 


begin 


NISP 
NKEC0N 
NDISP , 

2.3 SRB-II OVERLAY STRUCTURE 

The development of an overlay structure for SRB-II was required 
in an attempt to meet the MSFC UNIVAC 1108 core storage restriction 
of 32,000 words c.t any given time. 

The SRB-II overlay for the IBM 370 and UNIVAC 1108 versions of the 
program are shown i n Figures 2,5 and 2.6, respectively . One section 
(left side of tree shown in the figures) contains the Isp scaling 
modules and the necessary subroutines for these modules. The other 
section (right side of tree shown in the figures) contains the Internal 
BaHistic Module and necessary subroutines for this module. The 
routines. shown in the section with MAIN (Executive Module) remain in 
core at all times during a computer run. 
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FIGURE 2-5 SRB-II OVERLAY, TREE FOR THE IBM 370 COMPUTER SYSTEM 
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FIGURE 2-6 SRB-II OVERLAY TREE FOR THE MSFCUMIVAC 1108 COMPUTER SYSTEM 
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SECTION 3 


3.0 PROGRAfl INPUT OATA 

The input data for the SRB Performance Evaluation Model (SRB-II) 
consist of variable values compiled as an integral part of the program 
(BL0CK DATA), data input from punched cards, and data read from magnetic 
tape. All card inputs to SRB-II are arranged in input data packages 
and processed through the Centralized Input Module. Magnetic tape 
inputs are accessed by the subroutine requiring the data. BL0GK DATA 
subroutines are discussed in paragraph 3. T.l . Magnetic tape inputs _ are 
discussed in paragraphs 3.1.4 and 3.3. The remainder of Section 3 is 
devoted to a description of the punched card input required by SRB-II 
for successful operation of the program. 


3.1 DATA INPUT METHODS 

SRB-II uses four methods for transmitting data to the program modules. 

The data input methods are: 

1) BL0CK DATA 
2} NAMELIST 

3) Set Format Cards 

4) Set Format Magnetic TaDe 

3.1.1 BL0CK DATA 

The BL0GK DATA method consists of storing input data directly in the 
performance model using BL0CK DATA subroutines. The data in these sub- 
routines are stored in the computer memory during program compilation 
and are available to the module requiring the data. The BL0CK DATA 
subroutines form an integral part of the SRB performance model and 
consist principally of tabulated data for independent and deoendent 
parameters. Nominal motor performance parameter values, booster con- 
stants, program indicators, dispersion limits, and empirical coefficients 
are contained in the BL0CK DATA subroutines. The BLiiCK DATA method is 
augmented by the NAMELIST method. 

3.1.2 NAMELIST ■ 

The NAMELIST method is the primary method of entering input data into the 
model using punched cards. This method, when used in conjunction with ' 
BL0CK data' subroutines, is in most cases superior to field; fonaatted 
input methods because: strict arrangement of data is not required; 

the number of data cards is significantly reduced since only parameters 
subject to changes in value need to be entered; and field format 
violations are eliminated since both the name and value of the parameter 
are entered. NAMELIST can be used to override the assigned value of 
any BL0CK DATA yari able, thus avoiding the. need.to recompile the computer 
program for each BL0CK DATA change. 


D256-1 0020-2 


3.1 .2 (ContiriLied) 

A NAMELIST variciblG may be in the form of a subscripted array or a 
si ngle val ued variable. The val ue of the NAMELIST variabl e may be 
integer, real, literal, complex, hexadecimal, or logical. Examples 
of some of these variables are given belov/: 

Single valued variable: GAMMA =1.160, 

Single subscripted variable: THRlJST(l) = 24n00nn.0, or 

ISP (1) = 276.16, 

Multiple valued array; THRtlST(l) = 2400000.0, 2410000.0, 

In this example THRUST(2) = 2410000.0 

Logical variable: STATIC = .TRUE. 

3.1.3 Set Format Cards 

As noted in paragraph 3.1.2, the NAMELIST method for data input is 
superior to set format input in most cases. One example of a situation 
in which this is not true is the handling and sorting of large quantities 
of,alpha-numeric data. This situation occurs in inputs for the LEWIS 
subprogram; inputs consist largely of alpha-numeric data which identify 
the propellant chemical composition, ' The LEWIS input data are the 
only SRB-IT card input requiring set formats. 

A discussion of the data required for the LEWIS subprogram and the set 
format for these data are presented in paragraph 3.2.2'. 

3.1.4 Set Format Magnetic Tape 

Magnetic tape inputs to SRB-II are restricted to reconstruction data 
and LEWIS thermodynamic data. A discussion of data required for re- 
construction and the set format for this tape is presented in paragraph 

3.3. 

Execution of the LEWIS subroutine to obtain theoretical rocket perfor- 
mance parameters requires that the LEWIS THERMO data tape be available 
to the program. This tape is included in the MSFC tape library as 
tape number 23810. The format and contents of this tape are given in 
Appendix A. ' ' 

3.2 PROGRAM INPUT DATA PACKAGES 

Card inputs for SRB-II are arranged in data packages. These data packages 
are processed and transferred to disk by the Centralized Input Module, 
Incorporation of this data package concept gives the user flexibility 
in assembling and arranging the input data required for execution of’ 
the program. 
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3.2 (Continued) 

Data for SRB-II are arranged in as many as four data packages per case. 
These packages are shov/n graphically in Fioure 3 -I for a case data 
deck. The Identification Control Cards for each data package must 
be arranged in the data package in precisely the order shown; however, 
the data packages may be input in random order for each case. Data 
packages not required for a case may be excluded in the data deck^ 
for that case. It is also permissable to include a data package in 
the case data deck which is not required for that case (dataware 
ignored by the program). /' "CASE END" card (see Table 3-1 ) is required 
as the last card for each SRB performance case data deck. 

3.2.1 Control Data Package 

The Control Data Package is the input route by which data for the 
majority of SRB-II modules and subprograms are transmitted to the program 
(Exceptions: data for the LEWIS Subroutine, the One-Dimensional Nozzle 

Subroutine, the Nozzle Submergence and Contour Effects Module, and the 
Internal Ballistics Module are input in separate data packages). Program 
indicators which control the path followed through the program and options 
exercised in the program are also input through the Control Data Package. 
(Exception: Internal Ballistics and Nozzle Submergence and Contour Effects 

indicators input in the IBM Data Package). The Control Data Package must 
be included in every case . 

Control Data Package contains only the NAMELIST INPUT!. - Mo other 
NAMELISTS or set format data are included irf the data package. Variables 
input through the Control Data Package are identified and defined in 
the Program Options Section (paragraph 3.4). 

3.2.2 Lewis Data Package 

The Lewis Data Package is required for use when the option to exercise 
the LEVfIS Subroutine is chosen (NLEWIS > 0). The LEWIS Subroutine is 
used to calculate the theoretical values of Isp, c* and nozzle exit 
pressure, as well as several other theoretical solid rocket motor 
performance parameters. LEWIS Subroutine is a modification of the CEC71 
computer program developed by NASA Lewis Research Center (documented in 
Reference 3)'. The CEC71 program was modified to fit the requirements of 
the SRB performance model by remoying all options not require.d to analyze 
rocket problems. The program v/as then reduced in an effort to meet the 
core size requirements of the MSFC UNIVAC 1108 operating rules. 

The Lewis Data Package consists of set format input data cards. This 
data package is the only card input in the SRB model which is in set 
format rather than NAMELIST format. Formatted input data are required 
to provide the user with an efficient manner of supolying the needed 
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FIGURE 3-1 INPUT DATA PACKAGE ARRANGEMENT 
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TABLE 3-1 IDENTIFICATION CONTROL CARDS FOR DATA PACKAGES 


REQUIRED WORDS 


COMMENTS 



C0NTROLfjDATA 

LEWIS0DATA 

END 

IBM0DATA 

IDN0ZL0DATA 

CASE0END 


IDENTIFIES CONTROL DATA PACKAGE 

IDENTIFIES LEWIS DATA PACKAGE 

SIGNIFIES THE END OF LEWIS DATA PACKAGE 

IDENTIFIES THE IBM DATA PACKAGE 

IDENTIFIES THE ONE-DIMENSIOMAL NOZZLE DATA PACKAGE 

SIGNIFIES END OF DATA FOR A CASE 


AIT Identification Control cards start in card column 1. 
After the Required Words are punched, the remainder of 
'the card is available to the user for data identification 
and comments. 

NOTE: b is a Blank Character 

0 is the alphabetic letter 0. 
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3.2.2 (Continued) 

parameters to the program. Data cards in the package are identified as 
follov/s: 

Package Identification Card 
REACTANTS Card and Species Cards 
INSERT Cards 
0MIT Cards 

END Card ■ 

A typical Levn's Data Package is shown in Figure 3-2. 

Package Identification Card - The package identification card is required 
to let the program know that Lev/is subroutine data are available and to 
flag out this data to the input processor. The package identification 
card is the first card in the package and consists of the words LEWISji 
DATA punched in card columns 1 through 10. Card columns 11 through 80 
are available to the user for free format identification of the package. 


REACTANTS Card and Species Cards - The first card in this set contains 
the word REACTANTS punched in card columns 1 to 9. The last card in the 
set is blank. Between the first and last cards may be any number of 
cards (up to a maximum of 15), one card for each reactant species being 
considered. The cards for each reactant species must give the chemical 
formula, the relative amount of the reactant, and the initial enthalpy 
value. The format and. contents of the cards are summarized in. Table 3. II , 
A list of some REACTANTS cards is given in Table 3. HI. A set of example 
REACTANTS and Species cards are shov/n in Figure 3.3. 

The relative amounts of reactants may be specified in several v/ays. They 
may be specified in terms of moles, mole fraction, or mole percent (by 
keypunching M in card column 53) or in terms of weight, weight fraction, 
or weight percent (blank in column 53). A letter F is required in card 
column 72 of all reactant species cards. This indicates to the program 
that this is a solid rocket problem. 


Enthalpies are taken from the reactant species cards unless zeros are 
punched in card columns 37 and 38. For each card with the "00" code, 
an enthalpy will be calculated for the species from tape input data for 
the temperature given in card columns 64 to 71. The automatic calculation 
option is available only if tv/o conditions are met: 

(T) The reactant must also be one of the species in the set of THFRMD 
data. For example, NH 3 (g) is in the set of THERMO data and the 
enthalpy for this species can be calculated, NH 3(1 ) does not 
appear in the list and enthalpy cannot be calculated. 
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table 3- II 

REACTANTS CARDS 



CONTENTS 

FORMAT 

CARD COLUMNS 

REACTANTS 

3A4 

1 to 9 

One card for each reactant species 
(maximum 15), Each card contains: 

(1) Atomic symbols and formula num- 
bers (maximum 5 setsja 

5(A2, F7.5) 

1 to 45* 

(2) Relative v/eightb or number of 
moles 

F7.5 

45 to 52 

(3) Blank if (2) is relative weight 
or M if (2) is number of moles 

AT 

53 

(4) Enthalpy or internal energy^, 
cal /mole 

F9.5 

54 to 52 

(5) State: S,L,or G for solid, 
liquid or gas, respectively 

A1 

63 

(6) Temperature associated with 
enthalpy in (4) °K 

F8.5 

64 to 71 

(7) F if fuel or ^ if oxidant 

A1 

72 

(8) Density in g/cm^ (optional) 

F8.5 

73 to 80 


Blank 


Program will calculate the enthalpy or internal energy f4l 
for species in the THERMO data at the temperature (6) if 
zeros are punched in card columns 37 and 38. 

ill total fuels or oxidant in total 
oxidants. All reactants must be given either all in relative 
weights or all in number of moles: m reiai:ive 


*Atomic Symbols start in columns 1, 10, 19, 28 37. 




TABLE 3-HI LIST OF REACTAf ARDS FOR SOME OXIDANTS & FUELS 


CHEMICAL 

CHEMICAL FORMULA (CARD COlUMNS 1 to 46) 

Per- 

cent 

(cc 

54- 

62) 

ASSIGNED 
ENTHALPY, 
cal/mole 
(cc 54- 
62) 

(a) 

TEMPER- 

ATURE 

K 

(cc 64- 
71) 

(b) 

DENSITY, 
G/cm^ 
(cc 73- 
80) 

Acetonitrile 

Acetylene 

Air 

Aluminum 

Ammonia(g) 

Ammonia(l ) 

Ammonium percfil orate 

Aniline 

Argon 

Benzene 

Ben/llium 

Butane 

1-butene 

Chlorine(g) 

Chlorine(l) 

Chi ori ne trifl uoride(q) 
Chlorine tri fluoride (i ) 
Cyanogen{g) 

Cyanogen(l) 

Diborane 

Ethane 

Ethyl alcohol 
Ethylene 
Ethylene oxide 
Ethylene pol.ymer 
Fluorine(g) 

Fluorine(l) 

Granhite 

Helium 

Heptane 

Hydrazine 

C 2. H 3. Nl. 

C 2. H 2. 

N 1. 561 76>i .41959 AR.009324G .000.300 

ALl . ‘ 

HI. H 3. 

il l . H 3. 

fl 1. H 4. CLl . fl a. 

C 6. H 7, M 1. 

ARl. 

C 6. |i 6. 

BEl . 

C 4. H 10. 

C 4. H 8. 

CL2. 

GL2. 

CLl.; F 3. 

CLl. F 3. 

C 2. N 2. 

C 2 , fl 2. 

B 2. H 6. 

C 2. H 6. 

G 2. 116. 0 1. 

C 2. 11 4 . 

C 2. H 4. 0 1. 

c 1 . : , H 2 . 

■ f: 2 / ■ 

E 2 . . 

C 1.: 

HEl . 

C 7. II 16. 

N 2 . 1 ! 4 . 

100. 

loo. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

inn. 

100 . 

100 . 

100. 

100. 

100. 

100. 

100. 

100. 

12800. 

49270. 

-28.2 

0.0 

-10970. 

-17090. 

-70690 

7100. 

0.0 

+11718. 

0.0 

-36080. 

-5800. 

0.0 

-5391 . 
-39000. 
-45680. 
+73840. 
67655. 
4970. 
-25008. 
-66370. 
8100. 
-18840. 
-6100. 
0.0 

-3098. 

0.0 

0.0 

-536.30. 

12100. 

L 

L 

G 

S 

G 

L 

S 

L 

G 

L 

S 

L 

L: 

G 

L 

G 

L 

G 

L 

L 

L 

L 

L 

L 

S 

G 

L 

S 

G 

L 

L 

298.15 

192.60 

298.15 

298.15 

298.15 
239.72 

2.98.15 

298.15 
298.15 
298.15 
298.15 
272.65 

266.9 
298.15 

239.09 

298.15 
284.55 

2.98.15 
252.01 
180,59 
184.52 

298.15 
169.44 
283.72 
298.15 
298.15 

85.02 

298.15 

298.15 

298.15 

298.15 

F 

F 

fl 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

0 

0 

0 

0 

F 

F 

F 

F 

F 

F 

F 

F 

fl 

0 

F 

F 

F 

F 

.7857 

.610 

2.702 

.676 

1-95 

1.02173 

.8737 

1.85 

.6012 

.6263 

1.56 

1.8517 

.9537 
, 4.371 
.5464 
.7893 
.5683 
.8824 
.935 

1.505 

2.25 

.67951 

1.0036 


^Phaser S, solici; L, liquid; G, gas. 
Fuelj F; oxidant, j/f. 


Eased on the following molar percents: Ng = 78 . 0881 , 5(2 = 20.9495, Ar = 0.9324, Gf^2 = 0.0300 
Estimate based on paraffin Hydrocarbon series. 


D256-1 0020-2 






PACKAGES 




D256-1 0020-2 


3.2,2 (Continued) 


(2) The temoerature T must be in the ranqe Tin,.,/!. 2 T < T^,• k X 1 ? 
w ere to is the temoerature ranqe of the THFRMO^riata! 


^s^do?^nlentld^•n®L^fS well as the temoerature ranqe for the data 
cumented in Reference 3 and is presented in Apnendix A. 

If’ISERT cards may follow the REACTANTS 

snecips iS JL iih is ODtional. They contain the names of Particular 
species in the libran/ of thermodynamic data for the snecific nurooses to 

INSERf nrc^rd rn? contains the v/ord ?if1IT (in card columns 1 to 4) or 
columns 1 to 6) and the names of from one to four soecips 
starting in columns Ifi, 31, 4fi, and 61. The names must he exactly the 
same as they annear in the THERMO data (see Anpendix A). 


^MIT cards - These cards list soecies to be omitted from the THERMO 

data. If fiMlT cards are not used, the program will consider 
as possible species all those species in the THERMO data 
which are consistent with the chemical system heino con- 
sidered. Occasionally it may be desired to specifically 
omit one or more species from cqnsiderations as possible 
species. This may be accomplished hv means of OMIT 
cards. ' ' ' 


jNSERT cards - These cards contain the names of condensed snecies only. 

They have been included as options for tv/o reasons. 


more important reason for includinq the 
piSERT card option is that, in rare instances, it is 
impossible to obtain convergence for rocket problems 
ivithout the use of an INSERT card. This occurs when, 
by considering gases only, the temperature becomes ex- 
tremely low (say several degrees Kelvin). In these rare 
cpes, the use of an INSERT card containing the name of 
the reouired condensed species will eliminate this kind of 
convergence difficulty. When this difficulty 
following message is printed bv the orogram: 

ERATURE IMPLIES CONDENSED SPECIES SHOULD HAVE 
ON AN INSERT CARD". • 


occurs, the 
"LOW TEMP- 
BEEN INCLUDED 


The second and less important reason is that if one 

that one or several narticular condensed species 
will be present among the final equilibrium compositions 
for the first assigned point, then a small amount of 
computer time can besaved by using an INSERT card. 

Those condensed snecies v'hose chemical formulas are in- 
cluded on an piSERT card i-/ill he considered by the program 
during the initial iterations for the first assigned point. 
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3.2.2 (nontinued) 

If the INSERT card \/ere not used, only gaseous species 
v/ould he considered during the initial iterations. Hov/ever, 
•after convergence, the program vvould automatically insert 
the appropriate condensed species and reconverge. For 
all other assigned points the inclusion of condensed 
species is handled automatically by the program. There- 
fore, it usually is immaterial whether or not INSERT cards 
are used for the purpose of savinn computer time. 

END card - The END card signifies the end of the Lewis Data Package 
to the program and consists of 1 card v/ith the word FND 
punched in card columns 1 through 3. The END card is the 
last card in the Lewis Data Package. 

NOTE: The LEWIS THERMO data tape {MSFC Tape i¥23810) must be available 
to the program if the LEWIS Subroutine option is exercised. The 
Internal I/O Unit Number is 4. 

3.2.3 IDMrtZL Data Package 

The IDN(?ZL Data Package is a required input item v/hen the option to 
exercise a BATES Isp prediction is chosen (MSI = 3). The data contained 
in the IDN0ZL package is used in the One-Dimensional Nozzle Subroutine 
to calculate losses caused by solid particles in the motor exhaust 
products. Data input through the IDN0ZL Data Package are contained in 
the NAMELIST GiNEDNZ. Mo other NAMELISTS or set format data are included 
in the data package. Variables input through the IDN(3ZL Data Package 
are identified and defined in the Program Options Section (paragraph 3.4). 

3.2.4 IBM Data Package 

The IBM Data Package is a required input item when the option to predict 
motor performance using the Internal Ballistics Module is chosen (NF=1). 
The IBM Data Package is also required when a reconstruction is performed.- 
The data contained in the IBM Data Package are used in the Internal 
Ballistics Module and the Nozzle Submergence and Contour Effects Module 
to predict or reconstruct motor performance. Data input through the 
IBM Data Package are contained in the NAMELIST IBDATA. No other NAME- 
LISTS or set format data are included in this data package. Variables 
input through the IBM Data Package are identified and defined in the 
Internal Ballistics Options Section (paragraph 3.5). 

Multiple. Case Restriction: For multiple case runs in which the 

Internal Ballistics (IBM J option is executed, the IBM Data Pack- _ 
age must be included in the first case data deck, even if the option 
is not chosen in that particular case. The IBM Data Package 
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3,2.4 (Continued) 

must also he Included in the data deck for each case in ivhich the 
TDM Ballistics option i? chosen. The program writes the entire 

IBM Data Package from the first case on disk and then reads it in 
each subsequent case for which the IBM option is exercised. Thus 
data in the subsequent cases need only to update the data of the ’ 
first case (if no update to the first case data is required, at least 
one variable must be included in the IBDATA NAMELIST in each case). 

3.3 RECOMSTRUCTinN DATA TAPES 

The magnetic tapes used for reconstruction of fliqht or static test 
consist of one file and are required to have 10 variables per record 
written in binary format. There is no limit on the number of records 
in the file. Variables are required in the following order: 


1. 

TIME 

Time from ignition signal (SECONDS) 

2. 

PHEAD 

Measured head-end pressure (PSIA) 

3. 

FLDCL 

Force indicated by load cells in static test (LBp) 

4. 

AVEHI 

Measured vehicle lonoitudinal acceleration in 
flight (g's) 

5. 

MM2EST 

Mass of the motor not being reconstructed 
during this run (LBn^]) 

6. 

M0RB 

Mass of the orbiter excludino the external tank 
propellant during fliqht (LBf^) 

7. 

MET 

Mass of liquid propellants in the external tank 
during flight (LB[,^) 

8. 

FSSME 

Thrust of space shuttle main engines durino flinht 
(Total in LBp) 

9. 

DRAG 

Total shuttle vehicle drag durino fliaht (LBp) 

10. 

FM2EST 

Estimated thrust of the other solid motor not beino 
reconstructed during this run (LRp) 


User Action Required : Static test requirements are set to zero for 

flight reconstruction and flight requirements are set to zero for static 
test reconstruction. 


3.4 PROGRAM OPTIONS 

The SRB-II Performance Evaluation Model is constructed such that 

the user has options to choose the method of calculation for several of 

the intermediate calculation and output variables. Choice of calculation 
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3.4 (Continued) 

metnod is largely dependent on availability of data. This section of 
the document defines th-. available options and ahows the arranoement 
of the data for successful simulations using SRB-II. 

The program options are presented graphically as an Options Tree in 
Figures 3-4 and 3-5. The Option Tree is used to define the available 
options and to identify the program input indicators v/hich are required 
to obtain the desired option. Once the program indicators are defined, 
all the required input data for that option can be defined. For 
example if a reconstruction is required, the program indicator NRFC{!)N 
v/ould have to be input as NREC^N =1 and data required for a reconstruction 
would have to be made available to the program. Mote that the program 
automatically assumes one option by assigning a value to the indicator 
for that option is no value is input; all data for the automatically 
set option are required input in this situation. A further example 
of use of the Options Tree for a full problem follows: 

PROBLEM : Suppose that the user requires a nominal prediction of 

motor performance using contractor sunplied data to 
scale Isp and c*. The user desires to use an internal 
ballistics simulation for perfnrmanco prediction. 

SOLUTION : ^Tom the Options Tree, it is required that the indicator 

NREC0N be set to 0 (automatically set); the indicator 
NSI=2 for contractor data scaling prediction; the 
LEWIS subroutine is desired for the theoretical rocket 
performance parameters and the indicator NLEWIS=1; 
the motor nozzle is not submerged, indicator SILSB-1.0 
(automatically set); the divergence loss is to be cal- 
culated and CLAME=0.0 (automatically set); the motor 
nozzle is conical, indicators ALPHA=0.0, QEX=0.0, 

XLAMBD=1.0 (automatically set); since no dispersion 
prediction is required, indicator NDISP=0 (automatically 
set); and since an internal ballistics prediction is 
desired, NF=1 (set automatically). 

The input variables required for each option are identified in following 
paragraphs. Variables are identified and defined by option for each 
available option. Thus for the second example, it would he necessary 
to assemble the data required for a Contractor Data Scaling Prediction, 
paragraph 3.4.2, (the Contractor Data Scaling Prediction paragraph 
defines data requirements for each of the suboptions, LEWIS, submergence^ 
loss, and divergence loss) and the data required for an Internal Ballistics 
Performance Prediction, paragraph 3.5. 

The program indicators NREC0N and NCASES are input to the program in 
the Control Data Package and are contained in NAMELIST INPUT!. The 
variable TITLE is also input through the Control Data Package and 
NAMELIST INPUTl. TITLE is the variable input used to identify the run. 
TITLE may contain up to 80 alpha-numeric characters enclosed in blocks 
of six in Hollerith fields (6H//////,) and separated by commas. 
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3.4 (Continued) 

EXAMPLE : 

TITLE = 6HSRB P,6HERF0RM,6HANCE f1,6H0DEL C,6HASE AL,6HPHA , 
6mo/15/,6H73 

3.4.1 Isp Scalino Prediction 

The option to oredict the end item delivered Isp and c* usino the Isn 
Scaling technique is exercised by setting the program indicator MSI = 1 . 
Data required for an Isp Scaling prediction are listed in Table 3-IV. 

Data which are listed under the indicator MSI = 1 are required for all 
Isp Scaling predictions. Data listed under NLEWIS •= 0 are optional 
and required only if NLEWIS = 0. Likewise data listed under NLEWIS = 1 
are optional and required only if the indicator NLEWIS - 1, The 
variable CLAME may be calculated from the nozzle geometry (CLAME=0.0) 
or input (CLAME>0.0). The data package and NAMELIST in which the 
variables are input to the program are identified in the last two 
columns of Table 3-IV. 

Table 3-V lists definitions of all input variables identified in 
Table 3-IV. Variables are listed in alphabetical order to facilitate 
use. 

3.4.2 Contractor Data Scaling Prediction 

The option to predict the end item delivered Isp and c* usinn contractor 
supplied data is exercised by setting the program indicator NSI = 2. 
Inputs which are required for scaling contractor data are listed in 
Table 3-VI. As shown in paragraph 3.4.1, data required for a particular 
option are listed under the indicator for that option. Definitions for 
all contractor data input variables are listed alphabetically in 
Table 3-VII. 

3.4.3 BATES Prediction 

The option to predict the end item delivered Isp and c* usina BATES 
program data is exercised by setting the program indicator NSI = 3. 

Input data required for the BATES prediction are listed in Table 3-VIII, 
As shown in paragraph 3.4.1, the data required for a particular option 
are listed under’ the indicator for that option. Definitions for all 
BATES input variables are listed alphabetically in Table 3-IX. 

3.4.4 Dispersion Prediction 

The option to disperse the nominal prediction is exercised by setting 
the program indicator NDISP > (1. Input data required for a dispersion 
prediction are listed in Table 3-X. As shown in Paragraph 3.4.1, the 
data required for a particular option are listed under the indicator 
for that option. Definitions for all dispersion prediction input variabl 
are listed alphabetically in Table 3-XI. 
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3.4.5 Thrust Scalinq Prediction 

The option to predict the end item performance parameters as a function 
of time usinn the Thrust Scaling technique is exercised by setting the 
program indicator NF = 0. Input data required for a Thrust Scalinq 
prediction are listed in Table 3-XII. As shown in paragraph 3.4.1, 
the data required for^a particular option are listed under the indicator 
for that option. Definitions for all Thrust Scalinq input variables 
are listed alphabetically in Table 3-XIII. 

3.4.6 Reconstruction 

The option to perform a reconstruction of motor performance using 
measured test data is exercised by settino the program indicator 
NREC0M = 1. Input data required for a reconstruction are listed in 
Table 3-XV. As shovm in paragraph 3.4.1 , the data required for a 
particular option are listed under the indicator for that option. 
Definitions for all reconstruction input variables given in Table 3-XV 
are listed alphabetically in Table 3-XVI. 

The BIAS DATA are used to manipulate the reconstruction data after they 
are read from tape. The BIAS subroutine must be recompiled to include 
any bias factors required for correction of transducer drift and zero 
shift in a tape parameter. 

CONTROL SYSTEMS DATA are used to correct the measured or calculated 
•thrust for effects of thrust vector control. Thrust is corrected 
to nozzle centerline value by subroutine C0NT. 

The Reconstruction Module uses the Internal Ballistics Module to perform 
the required match on measured head-end pressure. Thusj the IRfI Data 
Package must be included for a reconstruction run. 

3.4.7 Optional Output 

SRB-II has several output options v/hich are keyed by program indicators 
input through the data packages. A complete description of these 
output options is given in Section 4.0. 

Theoretical rocket performance from the LEWIS subroutine may be obtained 
by setting the indicator IRKT01=1 in the Control Data Package. The 
results of the One-Dimensional Nozzle subroutine may be obtained by 
setting the indicator IDNPRT=1 in the IDN0ZL Data Package. Extended 
print for the One-Dimensional Nozzle may be obtained bv setting the 
indicator IDBUG = 1 in the IDN0ZL Data Package. Four Internal Ballistics 
Module print options are available and are controlled by the indicator 
PRTFLG (see paragraph 3.5 for definition and values). PRTFLG is input 
through the IBM Data Package. 
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TABLE 3-IV INPUT DATA REOIJIRED FOP AN ISP SOALINP PREniCTIOH 


KEY: NREC0N = 0 


VARIABLE 

NAME 

TYPE 

VALUE IF 
NOT INPUT 

REQUIRED BY 
INDICATOR 

INPUT DATA 
PACKAGE 

NAMELIST 

NSI = 1 

INTEGER 

1 


G0NTR0L DATA 

INPUTl 

PCAVE 

REAL 

494.0 

NSI = 1 

II 

II 

DTE 

}| 

71.0 

II 

II 

II 

ETACS 

It 

0.994 

II 

11 

II 

AKOO) 

1) 

0.8543 

11 

It 

II 

AK(n) 

n 

0.0244 

II 

II 

11 

AK(12) 

11 

-0.00125 

II 

II 

II 

ARETE 

It 

6,0 

II 

II 

11 

0K(3) 

II 

1.0 

II 

II 

II 

SILSB 

» 

1.0 

II 

II 

** 

NLEWIS =0 

INTEGER 

0 

_ 

tl 

il 

AK(1) 

REAL 

246.62 

NLEWIS=0 

II 

It 

AK(2) 

II 

5.406 

II 

It 

II 

AK{3) 

It 

-0.144 

■1 

It 

It 

MC(7) 

II 

5144.0 

II 

11 

II 

CSCCfEFCl ) 

12 

4735.23 

ft 

II 

11 

CSCdEF(2) 

II 

0.09546955 

It 

II 

II 

CSC(7EF(3} 

II 

-0,359902X10-4 

II 

II 

II 

NLEWIS = 1 

INTEGER 

0 

___ 

II 

II 

LEWIS INPUT 
DATA* 

ALPHA- 

NUMERIC 

— 

NLEWIS=1 

LEWIS DATA 

NONE 

CLAME =0.0 

REAL 

0.0 


C0NTR0L DATA 

INPUT! 

AHALFE 

11 

.3054327 

CLAME=0.0 

11 

It 

ALPHA 

u 

0.0 

11 

It 

II 

QEX 

11 

0.0 

It 

11 

11 

XLAMBD 

It 

1.0 

II 

11 

II 

CLAME >0.0 

■ 

II 

0.0 

v., 

II 

II 


*SEE SECTION 3.2.2: "LEWIS DATA PACKAGE" 
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TABLE 3-V 

VARIABLE 

AHALFE 

AKO) . 

AK(2) 

AK(3) 

AK(IO) 

AK(n) 

AK(12) 

ALPHA 

ARETE 

CLAME 

CSC0EF(1) 


VARIABLE DEFIMITIONS FOR AN ISP SCALING PRFniCTIDfl 

PEFINITinN 

End item nozzle divergence half angle. 
(RADIANS) 

Curve fit coefficient for theoretical vacuum 
Isp, in seconds, versus expansion ratio at 
average nozzle throat total pressure. Data 
derived from LEWIS Output. Isp==AK(l) + AK(2)* 
ARETE + AK(3)*ARETE**2 

Curve fit coefficient for theoretical vacuum 
Isp, in seconds, versus expansion ratio at 
average nozzle throat total pressure. Data 
derived from LEWIS Output. Isp=AK(l) + AK(2)* 
ARETE + AK(3)*ARETE**2 

Curve fit coefficient for theoretical vacuum 
Isp, in seconds, versus expansion ratio at 
average nozzle throat total pressure. Data 
derived from LEWIS Output. Isp=AK(l) + AK(2)* 
ARETE + AK(3)*ARETE**2 

Semi -log to base e curve-fit coefficient for 
Isp efficioncy versus flowrate, nrsp = 

AKfTO) + AK(11 )*1 n(m) + AK(12)*1 n{m)**2 

Semi-log to base e curve-fit coefficient for 
Isp efficiency versus flowrate, nlqp = 

AK(IO) + AK(ll)*ln(m) + AK(12)*ln (mP*2 

Semi-log to base e curve-fit coefficient for 
Isp efficiency versus flowrate, msp = 

AK(IO) + AK(n)*ln(m) + AK(12)*1 n(m)**2 

Equivalent conical divergence angle of the 
contoured nozzle. (RADIANS) 

End item nozzle area ratio. (A^j^jy/A-pj^j^^^-p) 

End item divergence loss factor x, where 
X = (■ ^^^2^*^ ) and a is the nozzle divergence 
half angle in degrees. 

Curve fit coefficient for end item propellant 
theoretical c*, in ft/sec, versus nozzle 
stagnation pressu'-e in psia. c* = CSC!!tEF(l ) + 
CSC?)EF(2)*PN + CSCF»EF(3)*PN**2 
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TABLE 3-V 

VARIABLE 

CSC0EF(2) 


CSC0EF(3) 


DTE 

ETACS 

nc{7) 

^^K{3) 

PCAVE 

SILSB 

QEX 

XLAMBD 


VARIABLE DEFINITinNS FOR AN ISP SCALING PREDICTION 
{Continued) 

DEFINITION 

Curve fit coefficient for end item propellant 
theoretical c*, in ft/sec, versus nozzle 
stagnation pressure in psia. c* - CSCf»EF(l) + 
CSC?)EF(2)*PN f CSCJ!iEF(3)*PN**2 

Curve fit coefficient for end item nronellant 
theoretical c*, in ft/sec, versus nozzle 
staqnation pressure in psia. c* = CSC?»EF(1) + 
CSC0EF(2)*PN + CSC0EF(3)*PN**2 

End item nozzle throat diameter. (IN) 

End item theoretical c* efficiency. 

Theoretical c* of the end item propellant 
at end item average chamber pressure. (FT/SEC) 

Empirical correction for end item c* in SISCAL. 

End item average chamber pressure. (PSIA) 

End item motor delivered specific impulse 
correction factor for the effect of nozzle 
submergence. 

Angle of the nozzle contour of the nozzle 
exit plane referenced to the nozzle centerline. 
(RADIANS) 

Ratio of the theoretical specific impulse of 
a contoured nozzle to the theoretical specific 
impulse of an equivalent conical nozzle derived 
with the ATPAP program. 
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TABLE 3-VI INPUT DATA REQUIRED FOR A CONTRACTOR DATA SCALINR PREDICTinN 
KEY: NRECtJN = 0 


VARIABLE 


VALUE IF 

NAhE 

TYPE 

NOT INPUT 

NSI = 2 

INTEGER 

1 

PCAVE 

REAL 

494.0 

DTE 

II 

71.0 

NC(23} 

II 

4886.0 

AK(IO) 

II 

0.8543 

AK 11) 

II 

0.0244 

Al((12) 

II 

-0.00125 

mz) 

« 

1.0 

ARET 

II 

9.8 

ARETE 

II 

6.0 

AHALFl 

a 

.2617994 

PCT 

II 

494.0 

S I caw 

II 

241.0 

WD2 

IT 

10.0 

SILSB 

II 

1 .0 

NLEWIS = 0 

INTEGER 

0 

AK(1) 

REAL 

246.62 

AK(2) 

II 

5.^06 

AK(3) 

n 

-0.144 

AK(13) 

II 

246.62 

AK(14) 

II 

5.406 

•AK(15) 

II 

0.144 

PAMT 

a 

14.7 

nc(6) 

a 

5144.0 

CSC/)EF(1 ) 

II 

4735.23 

CSCQEF(2) 

11 

0.095469 

CSCjOEF(3) 

a 

-0.359902X10-4 

NLEWIS = 1 

INTEGER 

0 

iMC(4) 

REAL 

inoo.n 

LEWIS INPUT 
DATA* 

AL.-NUM. 


CLAflE = 0.0 

REAL 

0.0 

AHALFE 

II 

.3054327 

ALPHA 

II 

0.0 

QEX 

II 

0.0 

XLAMBD 

a 

1.0 

CLAME > 0.0 ‘ 

ii 

0.0 



CWTR0L DATA INPUTl 


NLEWIS 



LEWIS DATA NONE 
CQfITRDL DATA IMPUTl 


*SEE SECTION 3. : "LEWIS DATA PACKARE 
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table 3-VII 

variable 

AHALFE 

AHALFl 

AK(1) 


AK(2) 

AK(3) 

AK(io) 

AK(ll) 

AK(12) 

AK(13) 


AK(14) 


VARIABLE DEFINITIONS FOR A CONTRACTOR SCALINO PREOICTinN 

definition 

FrADhS)™^^’® divsrqence half anqla. 

a™i\™®*(RAD*flN\)°‘"'‘ half, 

PTelsS?e''“jatrrf®''?™f ‘hroat total 

output. Isp. 

thMretl^arf*’’'"'®?* rtf"' "»tOT 

. Ta*?s^ua 

AK(1Ha1(2)*^®r1tem"(314M 

Semi-log to base e curve fit coefficient 

AK(lo^4^;f^f;^^(Vf^=(f,5»^®fli*4Isp = 

AK(IO) + AK(ll)*ln{tn)+ AK(12)*1 n(in)**2^^^ 

Semi-log to base e curve fit copffiV-!Qn+ 
AK(7of+^M?n?2T'^/M flowrate, njsp = 

akho) + /\K(ll)*ln(m)+ AK{^2)>^n(m)**?. 

‘'°®ffloient for contractor test 
motor theoretical vacuum specific impulse 
in seconds, versus expansion ratio on ’ 
average nozzle throat total pressur^ Dstci 
denved from LEWIS outout 
AK(14)*ARET + AK(15)*SeT«2 

Curve fit coefficient for contractor test 
motor theoretical vacuum specific impulse 
An seconds, versus expansion ratio T 

Pressure? Data 

derived from LEWIS outout j. 

AK(14)*ARET + AK(15)“MEf;*2 
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TABLE 3-VII 

VARIABLEeDEFINITIONS FOR A CONTRACTOR SCALING PREDICTION 

VARIABLE 

DEFINITION 

AK(T5) 

Curve fit coefficient for contractor test 
motor theoretical vacuum specific impulse, in 
seconds, versus expansion ratio on 
average nozzle throat total pressure. Data 
derived from LEWIS output. Isp=AK(13) -f 
AK{14)*ARET + AK(15)*ARET**2. 

ALPHA 

Equivalent conical divergence angle of the 
contoured nozzle. (RADIANS) 

ARET 

Contractor test motor nozzle area ratio. 

^^exit^^throat^ 

ARETE 

End item nozzle area ratio. 

(^exit/^throat^ 

CLAME 

End item divergence loss factor x, v/here 
X - ( — — ) and a is the nozzle diveraence 


half angle in degrees. 

CSC0EF(1) 

Curve fit coefficient for end item nronellant 
theoretical c* in ft/sec, versus nozzle 
stagnation pressure in psia. c* = CSC0EF(1) + 
CSC0EF(2)*PN + CSC0EF(3)*PN**2 

CSC0EF(2) 

Curve fit coefficient for end item proDellant 
theoretical c*, in ft/sec, versus nozzle 
stagnation pressure in psia. c* = CSC0EF(1) 
CSC0EF{2)*PN + CSC0EF(3)*PN**2 

CSC0EF(3) 

Curve fit coefficient for end item orooellant 
theoretical c* in ft/sec, versus nozzle 
stagnation pressure in psia. c* = CSC0EF(1) + 
CSC0EF(2)*PN + CSC0EF(3)*PN**2 

CSTAR 

Predicted end item characteristic velocity 
(c*). (FT/SEC) 

DTE 

End item nozzle throat diameter. (IN) 

MC(4) 

Typical BATES test motor nozzle throat total 
pressure, set to 1000 psia. .(PSIA) 

MC(6) 

Theoretical c* of the end item propellant 
at 1000 psia chamber pressure. (FT/SEC) 
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TABLE 3-VII 

VARIABLE 

0K(2) 

PAMT 

PCAVE 

PCT 

QEX 

SIC0N 

SILSB 

WD2 

XLAMBD 


VARIABLE DEFINITIONS FOR A CONTRACTOR SCALING PREDICTION 
(CONTINUED) 

DEFINITION 

Empirical adjustment for end item motor 
efficiency used in CDS I. 

Contractor test ambient pressure. (PSIA) 

End item average nozzle throat total pressure. 
(PSIA) 

Contractor test motor nozzle throat total 
pressure. (PSIA) 

Angle of the nozzle contour at the nozzle 
exit plane measured from the nozzle center- 
line. (RADIANS) 

Contractor test motor soecific impulse. 
(SECONDS) 

End item motor delivered specific impulse 
correction factor for the effect of nozzle 
submergence. 

Contractor test motor flowrate. (LBS/SEC) 

Ratio of the theoretical specific impulse 
of a contoured nozzle to the theoretical 
specific impulse of an equivalent conical 
nozzle derived with the ATPAP program. 
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TABLE 3-VIII INPUT DATA REDUIRED FOR A RATES PREDICTinN 


KEY: NREC/)N = Q 


VARIABLE 

NAME 

TYPE 

VALUE IF 
NOT INPUT 

REQUIRED BY 
INDICATOR 

INPUT DATA 
PACKAGE 

NAflELIST 

NSI = 3 

INTEGER 

1 


C0NTR0L DATA 

INPUT! 

AHALFE 

REAL 

.3054327 

NSI = 3 

It 

11 

AL 

tt 

0.943 

II 

II 

a 

ARETE 

II 

6.0 

II 

n 

II 

DTE 

II 

71.0 

II 

II 

II 

MC(4) 

II 

1000.0 

II 

II 

tl 

MC(8) 

II 

-1.0 

II 

II 

II 

MC{9) 

II 

2.9583 

II 

II 

II 

MC(IO) 

11 

2.9583 

tc 

II 

II 

0K(1) 

U 

1.0 

It 

II 

II 

QAIE 

II 

1450.0 

II 

It 

It 

QBT 

II 

4200.0 

II 

It 

II 

WINE 

II 

15081.0 

1! 

II 

" 

WPE 

II 

1676366.0 

II 

II 

a 

PSTAG 

II 

1000.0 

II 

IDN0ZL DATA 

0NEDNZ 

TSTAGK 

II 

3321 .0 

tl 

U 

n 

TEXITK 

II 

1957.0 

11 

II 

a 

CHEN 

II 

-433.11 

II 

II 

tl 

EXEN 

II 

-1225.4 

tl 

II 

a 

GMUS 

11 

0.00005 

II 

II 

a 

PR0P 

tl 

1.0 

II 

II 

II 

ALPHA 

II 

0.6 

11 

II 

II 

CPL 

II 

8503.3 

II 

II 

II. 

•CPS 

11 

8610.7 

it 

11 

a 

HLM 

II 

4.7944X107 

II 

1! 

It 

HSM 

II 

3.5978X107 

II 

II 

a 

TM 

II 

4167.0 

tl 

II 

N 

PARWT 

It 

101.94 

II 

1! 

a 

GMSSP 

II 

249.7 

a 

II 

a 

ATABL(l) 

11 

0.001 

a 

II 

a 

EPS! 

It 

0.1 

II 

II 

a 

EPSM 

11 

0.008 

II 

II 

a 

DXO 

II 

0.001 

II 

II 

a 

RTABL(l) 

II 

0. 000001 

II 

II 

a 

WTM0LG 

II 

19.20641 

II 

It 

a 

SIBML =0.0 

II 

0.0 

— 


IHPUTl 

ARETB 

II 

8.82 

SIBNL 

II 

a 

DTB 

It 

1.926 

II 

H 

a 

PCBT 

11 

1000.0 

II 

II 

a 

PAMBT 

II 

14.696 

It 
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TABLE 3- VI II INPUT DATA REniflRED FOP A BATES PPEPICTinM (CONTini'Fn) 


KEY: NREC0N = 0 

NSI = 3 

VARIABLE 

NAME 

TYPE 

VALUE IF 
NOT INPUT 

REQUIRED BY 
INDICATOR 

INPUT DATA 
PACKAGE 

NAMELIST 

SIBNL = 0.0 CONT. 






MC(3) 

REAL 

8.82 

SIBNL=0.0 

CflNTRjOL DATA 

INPUTl 

MC(5) 

It 

0.017 

II 

11 

a 

WPBT 

M 

69.83 

II 

11 

it 

SIBT =0.0* 

n 

0.0 


11 

ti 

ETABT* 

J1 

0.0 

SIBT = 0.0 

11 

It 

SIBT > 0.0* 

• 11 

0.0 

— 

II 

a 

DPB =0.0* 

II 

0.0 


II 

11 

AK(4)* 

II 

1.4939 

DPB = 0.0 

II 

II 

AK(5)* 

n 

. 

0.62403 

II 

IS 

a 

AK(6)* 

II 

-0.0030535 

II 

11 

a 

AK(7)* 

II 

-4.278X10-4 

II 

II 

II 

AK(8)* 

II 

8.6194X10-6 

. tr 

II 

a 

AK(9)* 

II 

-4.634X10-« 

II 

II 

a 

DPB > 0.0* 

II 

0.0 

— 

II 

II 

SIBNL >0.0 

M 

0.0 

— 

11 

ti 

OPE = 0.0 

II 

0.0 

... . 

II 

II 

AK(4) 

II 

1.4939 

DPE = 0.0 

II 

II 

AK(5) 

11 

0.62403 

II 

II 

a 

AK(6) 

II 

-0.0030535 

II 

II 

11 

AK{7) 

II 

-4.278X10-^ 

II 

II 

ii 

AK(S) 

II 

8.61 94X10-6 

II 

M 

II 

AK(9) 

11 

-A.63AX10-S- 

u 

II 

II 

DPE > 0.0 

II 

0.0 

— 

II 

II 

MLEWIS =0.0 

IFJTERER 

0.0 


II 

It 

AKO) 

REAL 

246.62 

fILEWIS = 0 

II 

II 

AK(2) 

It 

5.406 

II 

11 

a 

AK(3) 

It 

-0.144 

It 

It 

ts 

CSTARl* 

II 

5144.0 

II 

II 

n 

PAMBT* 

II- 

14.696 

II 

It 

II 

MC(2)* 

II 

5184.0 

II 

II 

a 

MC(6) 

II 

5144.0 

II 

II 

II 

MC(24)* 

II 

5144.0 

"11 

II ■■ ■■■■■.' 

II 

*REOUIRED ONLY If 

SIBNL = 

0.0 
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TABLE 3-VIII INPUT DATA REOIIIPED FOR A PATES PREDICTION (CONTINUED) 
KEY: NRECf)N = 0 NSI = 3 


VARIABLE 


VALUE IF 

REQUIRED BY 

NAME 

TYPE 

NOT INPUT 

INDICATOR 

NLEWIS = 0.0 CONT 




CSC0EF(1) 

REAL 

4735.23 

NLEWIS = 0 

CSC0EF(2) 

(t 

0.09546955 

11 

CSC0EF(3) 

It 

-0.359902X10-4 

II 

NLEWIS =1 

INTERER 

0 

— 

MC(15) 

REAL 

68.2 

NLEWIS =1 

LEWIS INPUT 
DATA* 

AL.NUM. 


II 

CLAME = 0.0 

REAL 

0.0 

— 

ALPHA 

n 

0.0 

CLAME = 0 

QEX 

It 

0.0 

It 

XLAMBD 

It 

1.0 

II 

CLAME >0.0 

II 

0.0 





■ - 



LEWIS DATA NONE 
C;)NTRpL DATA INPUT! 


* SEE SFCTION 3.2.2: "LEWIS DATA PACKARE 










TABLE 3-IX 

variable 

AHALFE 

AK(1) 

AK(2) 

AK(3) 

Ai((4) 

AK(5) 

AK(6) 

AK{7) 


!)2uu-KW2n-2 


VARIABLE DEFINITIONS FDR A BATES PREDICTION 

DEFINITION 

End Item nozzle diveroence half anqlp 
(RADIANS) 

Curve fit coefficient for theoretical 
vacuum po, in seconds, versus exoansion 
ratio at 1000 psia chamber pressure. Data 
denyed from LEWIS Output. Isn = AKfl) + 
AK(2)*ARETE + AK(3)*ARETE**2 ^ ^ 

Curve fit coefficient for theoretical 
vacuum Iso, in seconds, versus exnansion 
ratio at 1000 psia chamber pressure Data 
denyed from LEWIS Output. Iso = AK(1) + 
AK(2)*ARETE + AK(3)*ARETE**2 

Curve fit coefficient for theoretical 
vacuum Isp, in seconds, versus exnansion 
ratio at 1000 psia chamber nressure. Data 
LEWIS Output. Isp = AK(1) + 
AK(2)*ARETE + AK(3)*ARETE**2 

Curve fit coefficient for solid particle meai 
mass_particle diameter in microns versus 
nople throat diameter in inches. MDP = 
AK(4) + AK(5)*D + AK(6)*D**2 + AK(7)*D**3 
AK(8)*D**4 + AK(9)*D**5 ^ ^ 

Curve fit coefficient for solid particle mear 
mass narticle diameter in microns versus 
nozzle throat diameter in inches, MDP = 
AK{4) + AK{5)*0 + AK(6)*n**2 + AK(7)*n**3 
AK(8)*D**4 + AK(9)*D**5 ^ 

Curve fit coefficient for solid particle mean 
n^fiss particlG rliameter In niicrons vgtsus 
nozzle throat diameter in inches. MDP = 

AK(4) + AK(5)*D + AK{6)*D**2 + AK(7)*n**3 
AK(8)*D**4 + AK(9)*D**5 

Curve fit coefficient for solid particle mean 

ni^iss partTcie dianster in iriicrons VGrsus 
nozzle throat diameter in inches. MOP = 

AK(4) + AK(5)*D + AK(6)*D'**2 -h AK(7)*D**3 
AK(S)*D**4 + AK(9)*D**5 ' 
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TABLE 3- IX 

VARIABLE 

AK(8) 

AK(9) 


AL 

ALPHA 

ALPHA 

ARETB 

ARETE 

ATHBL 


CHEN 

CLAME 


CPL 

CPS 


VARIABLE DEFINITIONS FOR A RATES PREDICTION (Continued) 

DEFINITION 


Curve fit coefficient for solid oarticle mean 
mass particle diameter in microns versus 
nozzle throat diameter in inches. MDP = 

AK(4) + AK(5)*D + AK(6)*D**2 + AK{7)*D**3 
AK(8)*D**A + AK(9)*D**5 

Curve fit coefficient for solid particle mean 
mass particle diameter in microns versus 
nozzle throat diameter in inches. MDP = 

AK(4) + AK(5)*D + AK(6)*D**2 + AK(7)*D**3 
AK(8)*D**4 + AK(9)*D**5 

Propellant aluminum loading normalized to 
16"/: _ 5^AL 

16.0 

Temperature exponent in viscosity versus 
temperature equation of IDNOZL. Used only in 
IDNOZL Namelist. Set to 0.6 

Equivalent conical divergence angle of the 
contoured nozzle. / 

BATES test motor nozzle expansion ratio. 

(Aexit/Athroat^ 

End item nozzle area ratio, 

A dimensionless absolute tolerance on the 
gas and particle velocities used ‘in IDNOZL 
Runge - Kutta integration. Suggested value: 
0 . 001 . 

BATES test motor chamber enthalpy. (KCAL/ 

100 GRAMS) 

End item divergence loss factor X*, v^here 

1 t pric 

X is ( ■ 2 ' ) ,and a is the nozzle divergence 
half-angle in degrees. 

Specific heat of the BATES test motor com- 
bustion product liquid particles. (FT^/SEC^ - 
DEG R) 

Specific heat of the BATES test motor combustion 
product solid particles. (FT2/SEC2 - DEG R) 
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TABLE 3- IX 

VARIABLE 

CSC0EF(1) 


CSC0EF(2) 


CSC0EF(3) 


CSTARl 

DPB 

DPE 

DTB 

DTE 

DXO 

ESPM 


ESPl 

ETABT 

EXEN 


VARIABLE DEFINITIONS FOR A BATES PREDICTION (Continued) 

DEFINITION 

Curve fit coefficient for end item propellant 
theoretical c* in ft/sec, versus nozzle 
stagnation pressure in psia. c* = CSC0EF(1) + 
CSG0EF(2)*PN + CSC0EF(3)*PN**2 

Curve fit coefficient for end item propellant 
theoretical c* in ft/sec, versus nozzle 
stagnation pressure in psia. c* = CSC0EF(1) + 
CSC0EF(2)*PN + CSC0EF(3)*PN**2 

Curve fit coefficient for end item propellant 
theoretical c* in ft/sec, versus nozzle 
stagnation pressure in psia. c* = CSC0EF(1) + 
CSC0EF{2)*PN + CSC0EF{3)*PN**2 

End Item motor theoretical characteristics 
exhaust velocity at the nozzle throat (c*). 
(FT/SEC) 

BATES test motor combustion product solid 
particle mean mass particle diameter. (MICR0N) 

End item combustion product solid particle 
mean mass particle diameter. (MICRflMS) 

BATES test motor throat diameter. (IN) 

End item nozzle throat diameter. (IN) 

A dimensionless integration interval in 
IDM0ZL suggested value: O.nni. 

A dimensionless tolerance limit used as the 
denominator by ID0ZNL in solving the total 
momentum equation integration interval 
Suggested value; n.ODR. 

A dimensionless tolerance on gas velocity 
at left end of bracketing solutions in 
IDN0ZL. Suggested value: 0.1. 

BATES test motor Isn efficiency (actual/ 
theoretical) for standard conditions (1000 
to 14.7 Dsi a expansion). 

BATES test motor nozzle exit enthalov. 

(KCAL/100 GRAMS) 
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TABLE 3- 
VARIABLE 
GMSSP 

GMUS 

HLM 

HSM 

MC(2) 

MC(3) 

MC(4) 

MC(5) 

MC(6) 

MC(8) 

MC{9) 

MC(IO) 


IX VARIABLE DEFINITIONS FOR A BATES PREDICTION (Continued) 

DEFINITION 


BATES test motor combustion product particle 
mass density (LB^^/FT^). 

BATES test motor combustion products viscosity 
(LBp/FT - SEC) 

BATES test motor combustion product particle 
enthalpy at the particle meltinq point with 
the particle in the liquid state. (FT2 /sEC2) 

BATES test motor combustion product particle 
enthalpy at the particle melting point with 
the particle in the solid state. (FT2/SEC2) 

BATES test motor c*. (FT/SEC) 

Area ratio for theoretical optimum expansion 
from moo psia to 14.7 psia. Available from 
optional LEWIS Output. 

BATES test motor standard chamber pressure. 
(1000 PSIA) 

(1 - x) where X is the divergence loss factor 

1 cnsct 

2 and a is the nozzle divergence half- 
angle. 

Theoretical c* of the end item propellant at 
moo psia chamber pressure. (FT/SEC) 

End item nozzle entrance slope in the con- 
verging section. (See conical nozzle sketch 
below.) 

End item nozzle throat radius. (See conical 
nozzle sketch below). (FT) 

End item nozzle throat exit radius (See 
conical nozzle sketch below). (FT) 
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TABLE 3-IX VARIABLE DEFINITIONS FOR A BATES PREDICTION (Continued) 
VARIABLE DEFINITION 


CONICAL NOZZLE SKETCH 


SLOPF 

ric(R) 


CnMVEPRINR S; TION I 

I 



THROAT 


MC(15) 

Standard pressure ratio (1000 to 14.7 psia). 
Set equal to 68.2. 

MC(24) 

End Item motor theoretical specific impulse 
at the nozzle throat. (SEC) 

OK(l) 

Empirical adjustment for calculated end item 
c* from BATES Module. 

PAMBT 

BATES test ambient pressure. (PSIA) 

PARWT 

BATES test motor combustion product 
molecular v-^eiqht. (GRAMS/MOLE) 

PCBT 

BATES test motor steady state chamber 
pressure. (PSIA) 

PR0P 

BATES test motor propellant fuel indicator; 
0 or 1 for aluminum; greater than 1 for 
beryllium. 

PSTAG 

BATES test motor steady state nozzle entrance 
total pressure (PSIA). 

QAIE 

End Item expended inerts heat of ablation. 
(BTU/LBM) 

QBT 

BATES test motor heat loss. (BTU) 
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TABLE 3-IX 

VARIABLE 

QEX 

RATBL 


SIBT 

SILSB 

TM 

TSTAGK 

WINE 

WPBT 

WPE 

WTMflLG 

XLAMBD 


VARIABLE DEFINITIONS FOR A BATES PREDICTION {Continued) 

DEFINITION 


Angle of the nozzle contour at the nozzle 
exit plane measured for the nozzle centerline. 
(RADIANS) 


A dimensionless relative tolerance on the qas 
and particle velocity and enthalpy used in the 
IDNOZL Runge - Kutta integration. Suggested 
value; 1 X lO-fi 


BATES test motor specific impulse at standard 
conditions. (1000 to 14.7 nsia expansion). 
(SEC) 


End Item motor delivered specific impulse 
correction factor for the effect of nozzle 
submergence. 

BATES test motor combustion product solid 
particle melting point (DER R). 

BATES test motor ‘steady state nozzle entrance 
total temperature. (DEG K) 

End item mass of expanded "inerts". (LBfl) 

BATES test motor propellant mass. (LB[v]) 

End item mass of initial propellant. (LBM) 

BATES test motor combustion product gas molecular 
weight. Available from optional LEWIS Output. 
(GRAMS /MOLE) 

Ratio of the theoretical specific impulse of 
a contoured nozzle to the theoretical specific 
impulse of an equivalent conical nozzle 
derived with the ATPAP program. 
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TABLE 3-X INPUT DATA REQUIRED FOR A DISPERSION PREDICTION 


KEY: NRECON = 0 


VARIABLE 

NAME 

TYPE 

VALUE IF 
NOT INPUT 

REQUIRED BY 
INDICATOR 

INPUT DATA 
PACKAGE 

NAMELIST 

NDISP = 0 

INTEGER 

0 

— 

CONTROL DATA 

INPUTl 

NDISP = 1 


0 

II 

II 


DISLIM(I) 

REAL 

0.0C0615 

NDISP = 1 

II 


NDISP = 2 

INTEGER 

0 

— 

II 

II 

DISLIM(l) 

REAL 

-0.000615 

NDISP = 2 

II 


NDISP =3 

INTEGER 

0 


II 

II 

DISLIM{2) 

REAL 

0.00612 

NDISP = 3 

II 


NDISP = 4 

INTEGER 

0 

— 

II 


DISLIM(2} 

REAL 

-0.00612 

NDISP = 4 


II 

NDISP = 5 

INTEGER 

0 




DISLIM(3) 

REAL 

0.00268 

NDISP = 5 

II 

II 

NDISP = 6 

INTEGER 

0 

— 

II 


DISLIM(3) 

REAL 

-0.00268 

NDISP = 6 


Ji 

NDISP = 7 

INTEGER 

0 

— 



DISLIM(4) 

REAL 

0.000065 

NDISP = 7 


ir 

NDISP = 8 

INTEGER 

0 

— 

II 

II 

DISLi;i(4) 

REAL 

-0.000065 

NDISP = 8 


II 

NDISP = 9* 

INTEGER 

0 

— 



DISLIM{5) 

-eal 

0,00084 

NDISP = 9 

n 

1' 

NDISP =10* ^ 

INTEGER 

0 



II 

DISLIM{5) 

REAL 

-0.00084 

NDISP = 10 

It 

1* 

NDISP =11* ^ 

INTEGER 

0 

■ — — ' 


1* 

DISLIM{6) 

REAL 

0.00393 

NDISP = 11 

II 

i; 

NDISP. = 12* ^ 

INTEGER 

0 

— 

II 

V 

D1SLIM(6) 

REAL 

-0.a0393 

NDISP = 12 


I-- 

NDISP =13 ■ 

INTEGER 

0 



II 

DISLri{7) 

REAL 

0.05001 

NDISP = 13 

II 



*Hon -Autoinated Disoersion-User Action Required; Dis''erse innut data external 
to program. 

-fDoes not affect simulation - Used for run identifiction only. ^ 
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TABLE 3-X INPUT DATA REQUIRED FOR A DISPERSION PREDICTION (CONTINUED) 


KEY: NRECON = 0 


VARIABLE 

NAME 

TYPE 

VALUE IF 
NOT INPUT 

REQUIRED BY 
INDICATOR 

INPUT DATA 
PACKAGE 

NAMELIST 

NDISP = 14 

INTEGER 

0 

■ ■ 

CONTROL DATA 

INPUT 1 

DISLIM(7) 

REAL 

0.05001 

NDISP =14 

II 

II 

NDISP = 15 

INTEGER 

0 



It 

11 

DISLIM(8) 

REAL 

0.3861 

NDISP = 15 

II 

II 

NDISP = 16 

INTEGER 

0 


11 

It 

DISLIM{8) 

REAL 

-0.3861 

NDISP = 16 

IJ 

11 

NDISP = 17 

INTEGER 

0 


It 

ft 

DISLIM(9) 

REAL 

0.0006 

NDISP = 17 

M 

>1 

NDISP = 18 

INTEGER 

0 


M 

tl 

DISLIMO) 

REAL 

-0.0006 

NDISP = 18 

11 

. ■ 

II 

NDISP =19 

INTEGER 

0 

— 

II 

11 

DISLIM(IO) 

REAL 

5.000 

NDISP = 19 

II 

II 

NDISP = 20* 

INTEGER 

■ 

0 

. 

II 

11 

DISLIM(IO) 

REAL 

-5.000 

NDISP = 20 

II 

II 

NDISP = 21* ^ 

INTEGER 

0 


11 

M 

DISLIM(Tl) 

REAL 

0.01033 

NDISP = 21 

' n . 

ir 

NDISP = 22 

INTEGER 

0 

V, ■ ■ 

ti 

II 

DISLIMO 1) 

REAL 

-0.01033 

NDISP = 22 

• 

II 

' 

II 


flOTE: DISPERSION PREDICTIOF^S ALSO RECUIRE NOf'.INAL DATA DECK FOR INTERNAL 

BALLISTICS PREDICTION AS: INPUT. - v - 
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table XI 
VARIABLE 
DISLIMO) 

DISLIM(2) 

DISLIMO) 

DISLIM(4) 

DISLIM(5) 

DISLIM(6) 

DISLIM(7) 

DISLIM(8) 

DISLIMO) 

DISLIM(IO) 

DISLlM(n) 


VARIABLE DEFINITIONS FOR A DISPERSION PREDICTION 

DEFINITION 

Propellant density dispersion limit formed 
as the dispersion value of density divided 
by the nominal value of density. 

Burning rate law pressure exponent dispersion 
limit formed by the dispersion value of pressure 
exponent divided by the nominal value of 
pressure exponent. 

Burning rate law pressure coefficient dispersion 
limit formed by the dispersion value of pressure 
coefficient divided by the nominal value of 
pressure coefficient. 

Characteristic velocity dispersion limit 
formed by the dispersion value of characteristic 
velocity divided bv the nominal value of 
pressure coefficient. 

Propellant grain length dispersion limit 
formed by the dispersion value of grain length^ ' 
divided by the nominal value of grain length. 

Propellant grain web thickness dispersion 
limit formed by the dispersion value of grain 
web thickness divided by the nominal value of 
web thickness. 

Initial nozzle throat diameter dispersion 
limit. (INCHES) 

Initial nozzle exit diameter dispersion limit. 
(INCHES) 

Nozzle throat erosion rate dispersion limit 
based on the throat radius. (INGHES/SECOND) 

Propellant grain temperature dispersion limit. 
(DEGREES F) 

Initial "inert" mass dispersion limit formed 
by the dispersion value of "inert" mass divided 
by the nominal value of propellant mass. 
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TABLE XI VARIABLE DEFINITIONS FOR A DISPERSION PREDICTION (Continued) 

VARIABLE DEFINITION 

NDISP Indicator for dispersion orediction: 

NDISP=0, Nominal Prediction 
NDISP=1 , +propel1ant density dispersion 
NDISP=2, -propellant density disoersion 
NDISP=3, +burning rate law pressure exponent 
dispersion 

NDISP=4, -burning rate law pressure exponent 
dispersion 

NDISP=5, -fburning rate law pressure coefficient 
dispersion 

NDISP=6, -burning rate law pressure coefficient 
dispersion 

NDISP=7, +characteristic velocity dispersion 
NDISP=8, -characteristic velocity dispersion 
NDISP=9* +propellant grain length dispersion 
NDISP=I0*, -propellant grain length dispersion 
NDISP=11*, +propellant grain web/thickness 
dispersion 

NDISP=12*,' -propellant grain web/thickness 
disperson 

NDISP=13, +initial throat diameter dispersion 
NDISP=14, -initial throat diameter dispersion 
NDISP=15, +initial exit diameter dispersion 
NDISP=16, -initial exit diameter dispersion 
NDISP=17, +throat erosion rare dispersion (radius) 
NDISP=18, -throat erosion rate dispersion (radius) 
NDISP=19*,+Dropellant grain temperature dispersion 
NDISP=20*,-propellant grain temperature dispersion 
NDISP=21 5 +initial "inert" mass consumable 
dispersion 

NDISP=22, -initial "inert" mass consumable 
dispersion 


*Non-automated dispersion - user must perturb the input data in order to 
generate dispersion prediction data for this case. 
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TABLE 3-XII INPUT DATA REnUIRED FOR A THRUST SCALINP PPEPICTIOM 


KEY : NREC0N 0 


VARIABLE 

NAHE 

TYPE 

VALUE IF 
NOT INPUT 

REQUIRED BY 
INDICATOR 

INPUT DATA 
PACKAGE 

NAflELIST 

NF = 0 

INTEGER 

1 


cQntri!)l data 

INPUTl 

ATB(I)* 

REAL 

TABLE 3. XIV 

NF=0 

II 

II 

AFSRH(I)* 

II 

11 

II 

II 

II 

mpn?5m 

II 

1676366,0 

11 

II 

11 

NP(1)** 

INTEGER 

40 

II 

rr 

ir 

SII 

REAL 

190.0 

II 

II 

II 

MC(1) 

II 

0.1718X10-3 

II 

n 

II 

AMPRR{I)+ 

II 

TABLE 3. XIV 

II 

II 

II 

AMIRR(D + 

II 

II 

tr 

It 

II 

NP(10)++ 

INTEGER 

10 

II 

II 

II 

MIT0T 

REAL 

15081.0 

II 

II 

II 

NTB = 0 

INTEGER 

n 

— 

II 

II 

NTB > 0 

II 

0 

— 

(I 

II 

TB 

REAL 

150.0 

riTR>o 

11 

II 

NflP = 0 

INTEGER 

n 

— 

11 

II 

f!?'P > 0 

II 

0 

— 

11 

II 

HPT0T 

■ 

REAL 

1 676366. f’ 

NMP>n 

II 

II 


*liax1pium of 50 onints, +Hax1inum of 10 points. 

**fiP(l) must be 50. ++MP(10) must he ^ IP. 
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TABLE 3-XIII 

VARIABLE 

AFSRM(I) 

AMIRRd) 

AriPRR(I) 

ATB(I) 

MC(1) 

MIT0T 

MPT{9T 

NF 

NHP 

NP(1) 


VARIABLE DEFINITIONS FOR THRUST SCALING PREDICTIONS 

DEFINITION 


Array of nominal vacuum thrust used to generate 
a thrust seal 1 no oerformance prediction. Points 
in ATB{I) must correspond to points in AFSPM(I), 
i.e., ATB(3) is the time at which nominal thrust 
is AFSRM{3). (LBF) 

Array of non-dimensional inert mass remaininci 
for an inert mass calculation. Points in 
AMIRR(I) must correspond to points in AHPRR(I), 
i.e., AMPRR(5) must correspond to AMIRR(5). 

Array of non-dimensional ized propellant remaining 
for an inert mass calculation. Points in /\HPPR(I) 
must corresnond to points in /'MIRR(I), i.e., 

AriPRR(5) is the portion o-P propellant mass 
remaining which corresponds to AMIRR(5). 

AMPRR(I) is non-dimensional ized by dividing 
the mass remaining schedule by MPT?)T. 

Array of nominal times used to generate a 
thrust scaling prediction. Points in ATB(I) 
must correspond to points in AFSRM(I), i.e., 

ATB(3) is the time at which nominal thrust 
is AFSRM(3). (SECONDS) 

Constant for converting SRM thrust to SRM 
nozzle stagnation pressure for a thrust scaling 
prediction. MC(1) = 1.0/{CpAt) (IN-2) 

Total "inert" mass on board at time eoual zero. (LBM) 

Nominal propellant mass for a thrust scaling 
performance prediction. This mass must be 
compatible with ATB and AFSRM. dBH) 

Total initial propellant mass for a thrust 
scaling performance prediction. (LRM) 

Indicator used to establish v/hich module is 
to be used to generate propulsion performance 
predictions. NF = 0: Thrust Scaling Module. 

NF = 1 ; Internal Ballistics Module. 

Indicator used for thrust scaling prediction. 

NMP = 0: MPTF^T = MPNJ^'N, NMP > D: MPTJ^T is 

an input variable. 

Number of points in the ATR and AFSRM arrays. 

Maximum of 50 points. 


3-41 




TABLE 3-XIII 

VARIABLE 

NP(IO) 

NTB 

SII 

TB 


VARIABLE DpiNITIONS FOR THRUST SEALINR PRRniCTinNS 
(Continued) 


DEFINITIPN 

Number of points in the AMPRP and AMIPR 
arrays, flaximuro of 10 noints. 

Indicator used for thrust scalinq performance 
prediction. NTB = 0: TB = ATB (NP(1) - 1). 
NTB >0: TB was input by NAMELIST, 

Specific impulse of the inert material 
exhausted from the motor, (SECONDS) 

Motor burn time for thrust scalinq performance 
prediction. (SECONDS) 
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TABLE 3-XIV NOMINAL VALUES FOR USE IN THRUST SCALING PREDICTION 


POINT NO. 

ATB(I) 

AFSRM(I)* 

AMPRR(I) 

AMIRR(I) 

1 

0.0 

0.0 

0.0 

0.0 

2 

1.0 

2,969,184 

0.0 

0.0 

3 

4.0 

3,059,184 

0.0 

■ 0.0 

4 

8.0 

3,199,184 

0.090 

0.150 

5 

12.0 

3,299,184 

0.160 

0.235 

6 

16.0 

3,389,184 

0.255 

0.-325 

7 

20.0 

3,459,184 

0.740 

0.670 

8 

24.0 

3,549,184 

0.850 

0.770 

9 

28.0 

3,649,184 

0.928 

0.870 

10 

32.0 

3,839,184 

1.000 

1.000 

n 

36.0 

3,849,184 



12 

40.0 

3,929,184 



13 

44.0 

3,859,184 



14 

48.0 

3,799,184 



15 

52.0 

3,749,184 



16 

56.0 

3,699,184 



17 

60.0 

3,649,184 



18 

64.0 

3,559,184 



19 

68.0 

3,489,184 



20 

72.0 

3,439^,184 



21 

76.0 

3,439,184 



22 

80.0 

3,439,184 



23 

84.0 

3,439,184 



24 

88.0 

3,439,184 



25 

92.0 

3,439,184 



26 

96.0 

3,419,184 



27 

100.0 

3,409,184 



28 

104.0 

3,379,184 



29 

106.0 

3,359,184 



30 

108.0 

2,409,184 



31 

112.0 

3,199,184 



32 

116.0 

2,569,184 



33 

120.0 

1 ,149,184 



34 

124.0 

679,184 



35 

128.0 

559,184 



36 

132.0 

529,184 



37 

136.0 

479 .184 


• 

38 

140.0 

459,184 



39 

150.0 

0.0 



40 

200.0 

0.0 




*DATA DERIVED FROM 260-SL-l TEST 
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TABLE 3-XV iflPUT PATA RROUIREP FOR RFCnflSTPaJCTIPM 


KEY: NRECCJN = 1 


VARIABLE 

NAME 

TYPE 

VALUE IF 
NOT INPUT 

REQUIRED BY 
INDICATOR 

INPUT DATA 
PACKAGE 

NAMELIST 

NRECOli = 1 

INTEGER 

0 

— 

C0NTR0L DATA 

INPIITI 

BIAS DATA* 

REAL 

0 

NREC0N=1 

II 

It 

CONTROL SYS. 

II 

0 

II 

« 

II 

DATA** 





STATIC 

LOGICAL 

.TRUE. 

II 

II 

II 

NGfZPCfS 

INTEGER 

0 

It 

II 

II 

^K{4) 

REAL 

100.0 

II 

II 

II 

0K(5) 

II 

6.0 

It 

It 

II 

LISTl = .TRUE. 

LOGICAL 

.FALSE. 

— 

« 

II 

TSREC 

REAL 

0.0 

LISTl =. TRUE. 

n 

II 

TEREC 

II 

lono.o 

II 

II 

11 

CSBAR 

II 

4800.0 

II 

II 

It 

PBAR 

II 

1000.0 

It 

II 

II 

TWEB 

II 

1000. 0 

It 

II 

11 

PITH 

II 

looonoo.o 

II 

II 

II 

LISTl = .FALSE. 


.FALSE. 


II 

II 

PFLAG 

REAL 

0.0 

.IST1=. FALSE 

n 

a 

PADJC(l) 

II 

0.0 

It 

a 

II 

PADJC(2) 

II 

o.n 

II 

II 

a 

PHEPI 


100.0 

SI 

II 

11 

ATBAR 

II 

3959.0 ■ 

II 

II 

II 

MPTOT 

■ 

II 

1676366.0 

a 

II 

II 

NLEWIS = n 

INTEGER 

■ 0 


II 

II 

CSCflEFd) 

REAL 

4735.23 

NLEWIS=0 

II 

II 

CSCfiEF{2) 

It 

0.09546955 . 

IS 

n 

II 

CSCgEF(3) 

II 

•0.359902X10"^ 

II 

II 

a 

NLEWIS = 1 

INTEGER 

0 

— 

•II 

11 

LEWIS DATA 

AL.-NUM. 


NLEWIS=1 

LEWIS DATA 

NOME 

PACKAGEt 





NPH = 0 

INTEGER 

0 

— 

IBM DATA 

IBDATA 

Mfigri 

REAL 

2000000.0 

NPH=0 

C0NTP.0L DATA 

INPUTl 

MVEHI 

U 

5000000.0 

II 

II 

II 

MIfiST 

It 

5000Q.O 

II 

11 

|l 


*P.equirBS ComDilation of Subroutine BIAS and Possible Adriitinn of v'ariables to ' 
Common and Namelist. 


**Requires Comnilation of Suhrniitine CPFIT and Possible Addition of Variables to 
Common and Namelist. 

+See Section 3.2.2: "Lewis Data Packaae". 
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TABLE 3-XV INPUT DATA REnUIRED FOR RECOMSTRUCTIflN (CONTINUEP) 


KEY: NRECON = 1 


VARIABLE 

NAME 


TYPE 


VALUE IF 
NOT INPUT 


REQUIRED BY 
INDICATOR 


INPUT DATA 
PACKAGE 


NAMELIST 


NPH > Q 


INTEGER 


IBM DATA 


IBDATA 


PHST(I) 

TIMEPH(I) 


REAL 

U 


0 

0 


MPH > 0 

II 


II 

rr 


IN ADDITION 
A NOMINAL PN: 
NSCE MODULE 


[TO DATA It 
EDICTIOM 
ARE REQUIRED 


STED ABOVE, ALL 
JSING THE INTERNS 
FOR A RECONST; 


INPUT DATA 
L BALLISTIC 
■RUCTION. 


REQUIRED FOR 
5 MDDULF AND 
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TABLE 3-XVI VARIABLE DEFINITIONS FOR A RECONSTRUCTION 


VARIABLE 

ATBAR 

CSBAR 

CSC^EF(l) 

CSC?)EF(2) 

CSCJ^EF(3) 

LIST! 

MINST 

f1M?)Tr 

MPT0T 

MVEHI 

NLEWrs 

N0ZP0S 


DEFINITION 

Average throat area for reconstruction. 

Normany the arithmatic averaoe of initial 
and final throat areas. (IN2) 

Characteristic velocity used with PBAR to 
shift the theoretical c* curve for use in a 
reconstruction. (FT/SEC) 

Curve fit coefficient for end item oropellant 
theoretical c* (FT/SEC) versus nozzle stagna- 
tion pressures (PSIA). CSCPEF(I) is intercept 
at pressure = 0.0. (c* = + c^P^jS + c^P^2) 

Curve fit coefficient for end item propellant 
theoretical c* (FT/SEC) versus nozzle stagnation 
pressure. 

Curve fit coefficient for end item orooellant 
theoretical c* (FT/SEC) versus nozzle stagnation 
pressure. 

Logical indicator which, when true, shows the 
following variables were innut: TSRFC, TEREC, 
CSBAR, PBAR, TWEB and PITW. 

Mass of SRB inert structure, not including 
motors, (LBM) (FLIGHT ONLY) 

Initial total mass of the motor being recon- 
structed including propellant, case, insulation, 
nozzle, etc. (LBM) 

Total propellant mass loaded into the motor 
being reconstructed (LBM) 

Mass of Shuttle vehicle at launch. (LBM) 

(FLIGHT ONLY) 

Indicator for a Lewis Run. NLEWIS = 0: No 
LEWIS Subroutine execution, MLEWIS = 1: LEWIS 
Subroutine execution. 

Indicator for nozzle exhaust direction for a 
static test reconstruction. N(?iZP0S = -1 
Horizontal; Nf)ZP(as = n, lip; N07PJ?'S = 1 
Down, ’ 
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TABLE 3-XVI VARIABLE DEFINITIONS FOR A RECONSTRUCTION (Continued) 


VARIABLE 

DEFINITION 

NPH 

Number of points in the arrays of head-end 
pressure, PHST(I), versus time, TIMPH(I). 
Maximum of 70 points. 

MREC0N 

Indicator used to establish a prediction or 
reconstruction. NREC0N = 0: Prediction, 

NREC0M = 1: Reconstruction. 

0K{4) 

Characteristic velocity correction increment 
for Reconstruction. (FT/SEC) 

!3K(5) 

Iteration limit for a reconstruction. 

PADJC(l) 

Slope of pressure adjustment term in recon- 
struction. Used to adjust measured head- 
end pressure to nozzle stagnation pressure. 
Pns = Ph " (PADJC(1)*TIME + PADJC(2)) 

PADJC(2) 

Intercept of pressure adjustment term in 
reconstruction. Used to adjust measured head- 
end pressure to nozzle stagnation pressure. (PSIA) 
Pns = Ph “ (PADJC(l)niME + PADJC(2)) 

PBAR 

Time average nozzle stagnation pressure used 
with CSBAR to shift the theoretical c* curve 
for a reconstruction.' (PSIA) 

PFLAG 

Pressure level at which initial reconstruction 
calculations are started and ended. (PSIA) 

PHEPI 

Percentaae of head-end pressure integral which 
defines tWEB. (%) 

PHST(I) 

Array of head-end total pressure during a 
reconstruction. This is the dependent array; 
See TIMEPti(I) as independent array. (PSIA) 

PITW 

Head-end pressure integral used to calculate 
TWEB. (PSIA-SEC) 

STATIC 

Logical indicator used in a reconstruction to 
establish a static test or flight. STATIC * 
.TRUE.: Static Test. STATIC = .FALSE.: Flight. 

TEREC 

Time to end integration of pressures in a 
reconstruction. Set equal to time when the 
head-end pressure tail off reaches PFLAG. (SEC) 
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TABLE 3-XVI 

VARIABLE 

TIMEPH(I) 

TSREC 

TWEB 


VARIABLE DEFINITIONS FOR A RECONSTRUCTION (Continued) 

DEFINITION 


Array of time during a reconstruction. In- 
dependent array, PHST(I) is the deoendent array. 
(SEC) 

Time to start integration of pressures in a 
reconstruction. Set equal to time when the 
head-end pressure buildup reaches a value 
equal to PFLAG. (SEC) 

Web time of a grain for a reconstruction. 

Define as the time of occurance of a percentage 
of the total head-end pressure integral. (SEC) 
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3.4.7 (Continued) 

option^'s^exerdsed performance prediction 

indicator NCARD=1 thrcSrtrol olta1a?L"S!*"“ 

S?edio‘«on^by S0ttiJrth™iSdi-cSsVAPr'f?/r^’^ performance 

Unit 12) and NPL0T=1 fnlnt t^no NTAPE=1 (data tape, Internal I/O 

NPW a^e inpuMJrL^K'?h‘e*cS„=tr^rD«a^'p^(Lge!‘ 

3-5 INTERNAL BALLISTICS OPTIONS 

Xh??ai‘iyi„“ ??gi?f3:r *ui„j"„‘rs’ 

>E.!v':rs^i;r;.rars'C«“! 

3,5.1 Grain Geometry 

heL“?I^f’e thP ™otor is divided into three sections aft 

Ss 

arm?nv L n? ^ ‘ ^ cylindrical section may be divided into 

as many as one hundred mass addition regions. -v.uku moo 

sections, various options are available Tvnical 
motor configurations that can be analyzed are shown in Figure sT 
0 or designs with both internal and external taper, monolithic and 
segmented propellant sections can be analyzed. The cSoab it L 
imitations,^and inputs for the aft-head, fore-heSraSdcW^ 
section are discussed in the following paragrapSr ^^'indncal 
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FIGURE 3-6 INTERNAL BALLISTICS MODULE OPTIONS TREE 
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FIGURE 3-6 INTERNAL BALLISTICS MODULE OPTIONS TREE (CONTINUED) 
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INDICATORS ASSUHSO TO BE ZERO UNLESS SCT 
QTHERtllSEi ALSO. THE ORDER IN 
OPTIONS ARE PRESENTED DOES NOT FOLLOW A 
CHRONOLOGICAL SEQUENCE. 


FIGURE 3-6 INTERNAL BALLISTICS MODULE OPTIONS TREE (CONTINUED) 
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3, 5. 1,1 Aft-Head Section 

Two grain geometry options are available: a straight through grain, 

where the port cross-section is the same as that at the aft tangent 
plane, or an input table of burning area versus distance burned, 
ihe latter is generally used when a submerged nozzle or complex grain 
desips are considered. A typical straight through grain configuration 
IS shown in Figure 3-10. The required grain geometry” input at the aft 

discussed below under the heading "Cylindrical Section", 
section 3. 5. 1.3. Inputs required for simulating a submerged nozzle are 
discussed in Section 3.5.5. 

UD by a combination of elliptical 
and conical configurations. Figure 3-11 shows necessary inputs to 
describe the aft dome shape. ‘ ^ 

and 4) DLRF. The term 

moto^Jiq ellipse semi-axis perpendicular to the 

motor axis to the semi -axis coincident with the motor axis The case 

aMe^anale^bPtwLn%h'’®r*^^'^^^ maximum allow- 

LnteHile ellipse section and the motor 

centerline at the plane of the case opening is APNM To dpfinp thp 

end section, the program first computes the value of angle ANR 
This angle is defined by the tangent to the ellipse and the moW axis 

- value of ANR is less tLn o? 

teminatP^S? fh elliptic contour is assumed to 

flWD 3c ^ I the case opening diameter. If, however, the value of 

if temina?pfaJ value of Ap, then the elliptic contour 

thp pn3pf^^ between the tangent to 

the ellipse and the motor axis is equal to the value of APM. From 

S ^ frustum of a cone tangent to the ellipse is extended 

to the case opening diameter. 


Desired accuracy of the burning surface area and initial volume are 
gopned by the input DLRF. As this value approaches zero, thi 
calculations approach an exact solution at any value of thickness 
burned. A typical value is 0.01. 

Certain limiting requirements must be satisfied for the end dimensions 
as TO I lows: 


1. The value of DNI/2 must be less than or equal to the value of RF 
in the reference plane adjacent to the end section. 

2. The value of A0NM must be greater than zero. 

3. The value of BN must be either 

a. Equal to 1.00 (hemispherical end), 

b. Equal to or greater than 1.01 (oblate spheroid end) or 

c. Equal to or less than 0.99 (prolate spheroid end). 
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3. 5. 1.1 (Continued) 

configuration shown in’Figure 3-l" i1?ust?aSL and ^irSSn^• 

is feL'thai Sm 1' J-® ’■^=“'*’■"9 lenflth of the aft dome 

IS less than the maximum v/eb thickness of the aft tanoent olanp 

added Sf TL*'-^ difference between ?hesS fSs S]?l be 

thickneL Thi^irLS" it eoual to the maximum web 

sectiSn! Prevents the aft face from burning into the cylindrical 

JL^SenMiroH"" “Ptioo is used, BN, DNl , A?NM, and DLRF are 

d‘ J^P^r^nt S^SV^?1hral?fh^If=-^iSra"^^ ?h1 

thif5aL™alL"‘L™er^eJuired JnSS^ determined.- 

3, 5. 1,2 Fore-Head Section 

Three grain geometry options are available, Thev are strainht fhvnmmh 
grain burning area table, and a head-end with i%b! Thf? rsJ two 
are the same as discussed above in the aft-hLd section! 

Ilia straight through grain option are: 1) BH, 2) [IHl 

Definitions of these inputs and the correspondino 
limitations are identical to the equivalent aft-end section inputs 
and limitations. The fore-head increment sizes are set by the 
thp^fnloh ratio DLRF, -fche same as the aft-end increment sizes. In 

SLISI however, RF is the outside propellant radius of the forward 

iSh^/ calculation of the fore-head section 

JlSSeVlin geleJrJ ‘=™ent 

The third option is a head-end with web as shown in Figure 3-12 Re- 

«• S) DRW. and 6) AKK. 

T°£ d£fi£ition £fRS“£fr£f*°‘"' dP" PPly de an elliptical configuration. 
Thf + DH IS the same as that for a straight through grain 

The term BTA^E represents the ratio of the head-end with web ellipse ‘ 
semi -axis normal to the motor axis to the ellipse semi-axis which^lies 
on the motor axis. In addition to defining the propellant internal 
coutour, BTA/JE is used to transfer the program to the head-end "with 
web so ution.^ If the value BTAOE is non-zero, the head-end wit^ 
web solution s used. The term RIG is the radius of the ignUer 
opening, which is assumed to be inhibited on the internal surface 

to\e''cn?n!THLJ’'^ ^re always assumed 

be coincident. The reference fore-head length, HHR, is always measured 
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3, 5. 1.2 (Continued) 

from the center of the fore-head ellipse closure. A restriction 
to the head-end with web solution requires that the initial Dronellani- 
surface reach the igniter opening. If this restrictiSn is exceeded 
by using too small an input value for BTAf^E, the program will stop 
and print a statement that the head-end with web fraction is negative, 

initial volume calculations for the head- 
end vnth v/eb are based on approximate integrations of surface and 

Initial values of the required increment sizes, 

ORV and DL, are computed in the same manner as DL for a straiaht 

Subsequent calculations, the increment sizes are 
modified by the input constants AK and AKK, which adjust the increment 
sizes to account for changes in the configuration during burning. The 

Irrnl ^he primary factors in dkermining 

flc execution time for the head-end with web calculations. 

approach zero, the calculations approach exact solutions. 
Input value AK controls the accuracy of the burn surface area calcu- 

s the accuracy of the initial volume 
calculations. A recommended value for these terms is 0.50. 

lllTl P'f’onellant volume in both end sections at any thickness other 
tnan zero are computed by means of approximate integration of burning 

thickness burned. These integrated volumes’ 
are subtracted from the previous volume to give instantaneous volumes 

r.^1 i''T\ of this process, the final propellant weight 

at complete burnout is likely to be some value other than zero 
As the time increment input approaches zero, the calculation for in- 
stantaneous propellant weight aoproaches an exact solution. 

3. 5. 1.3 Cylindrical Section 


The cylindrical section is that section lying between the forward and 
• ^ planes. The cylindrical section length is designated as 

input HC^. Vanous grain geometry options exist as shown in Figure 3-13 
The grain may be either monolithic or segmented. 


The cylindrical section is described for input by dividing it into a 
number (maximum of 18) of geometrical reference planes, Figure 3-14. 
These reference planes are used in describing the grain geometry at 
a particular location. They are designated by AINCIN(1 ) . , .AINCIN(18) 
and are placed at desired distances from the forward tangent plane. 


■y 

r 


The value of AINCIN(l) must always be zero, thus placing the first plane 
at the forward tangent plane. Likewise, the AINCIN value of the last 
plane used roust be equal to the value of HC0, thus placing the last 
plane at the aft tangent plane. 
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3. 5. 1.3 (Continued) 


The general forked 

wagon wheel 

configuration for 

an input reference 

plane is shown in 

Figure 3-15. 

The following are 

required inputs: 

AN0(I) 

R4(I) 

ALSl(I) 

A01(I) 

RF(I) 

R5{I) 

ALS2(I) 

A02(I) 

TAUW(I) 

R6(I) 

ALA(I) 

A03(I) 

R2{I) 

R7(I) 

ALB(I) 

A04{I) 

R3(I) 

R8(I) 

ALE(I) 

A05(r) 


AN0 is the number of symmetrical sectors. The segment lengths ALC and 
ALD are not input. 

Figures 3-16 and 3-17 show the inputs and calculated plane constants 
for the general forked wagon wheel configuration. The calculated nlane 
constants allov/ the user a neans of checking if the orooer grain design 
is being analyzed. Although the program is designed to solve the 
general configuration shown in Figure 3-15, there are certain variations 
of this configuration that exceed the mathematical limits of the 
analysis. To obtain a program solution, all of the followina conditions 
must exist for each reference plane used: 

must be greater than or equal to zero 
must be less than or equal to T2M 
must be less than or equal to T12!'i 
must be greater than or equal to zero 

must be greater than or equal to zero 

must be greater than or equal to zero 

must be greater than or equal to zero 

must be greater than or equal to zero 

must be less than 90° and greater than or equal to A02 
must be greater than or equal to zero 
must be less than 90° and greater than or equal to A02 
must be greater than or equal to zero 
must be less than 9i;’ and greater than or equal to A04 
must be greater than or equal to zero 

If any of these restrictions are exceeded for any reference plane, the 
program will automatically stop and will print the reference plane 
dimensions along with a statement of the exceeded restriction. Although 
the program is self-checking for these mathematically invalid con- 
figurations, there are some physically invalid configurations " 
which there is no such check. For example, such a situation w.l i occur 
in a case where a value of an AN0 input is not a whole number, and the 
program will yield valid mathematical solutions that arc physically 
impossible. To avoid such problems, it is suggested that the user 
visually review a sketch of each reference plane input configuration 
before attempting to obtain a program solution. 

Special attention should be given to the manner in which input lengths 
ALA, ALB, and ALE are defined. The length ALA is measured along a common 


R1 

R3 

R7 

R9 

ALC 

ALD 

B71M 

B91M 

A01 

A0Z 

A03 

A04 

A05 

TH4 
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3, 5. 1.3 (Gontinued) 

tangent to the arcs generated by radii R2 and R3 (Figure 3-15). One of 
the points that defines this length is the point of tanqency vnth the 
arc generated by the radius R2. However, the other point that defines 
length ALA is not the point of tangency with the arc generated by 
radius R3. This end limit of length ALA is determined by bisecting 
the angle subtended at the center of the arc generated by radius R3 
and extending this bisector until it intersects the axis of symmetry of. 
the grain sector. From this intersection point, a perpendicular to 
line ALA is erected. The intersection of this perpendicular and the 
line ALA is the point which defines the end of length ALA. An expanded 
sketch showing the definition length of ALA is shown in Figure 3-18. 

Line lengths ALB and ALE are defined in a similar manner. 

Figure 3-19 shows the inputs required for the v/agon wheel grain 
configuration. 

Figure 3-20 shows the inputs required for the standard star grain 
configuration. 

Figure 3-21 shows the inputs required for the slotted-cone grain 
configuration. Burning area tables are always required for the fore 
and aft-head region with this grain conf^uration. ' In addition to the 
standard star input of Figure 3-20 an additional parameter TH?1 is required. 

Figure 3-22 shov;s the inputs required for a circular port grain con- 
figuration. 

The inert sliver option can be used for the wagon wheel, star, and 
slotted-cone grain configurations. If the inert sliver option is used, 
the distance burned on the inert sliver TSLVR(I) must be input. Inert 
slivers are considered only in t'le cylindrical section. TSLVR(I) is 
also required to be equal to or greater than the web thickness, TAUW(I). 

An input is required for the incremental distance burned before and 
after web time. These inputs are DTAU(I) and DTAUW(I) and are used in 
generating tables of geometry data versus distance burned for the 
reference planes and end sections. 

If there is a special grain configuration to be analyzed that the 
program options do not include, this can be analyzed by inputting a 
table of port perimeter versus distance burned in place of the grain 
geometry inputs. The initial port cross-sectional area is also required. 
When this option is used, the burning area versus distance burned is 
required for both the fore and aft-head sections. The required inputs 
are identical to the inputs that were described in the fore and aft- 
head sections. When this option is used, moments of inertia and centers 
of gravity for the fore-head, aft-head, and cylindrical sections will 
be bypassed unless they too are input for each reference plane and end 
section. 

The cylindrical section must also be described for input by dividing 
it by a number of increment dividing planes (a maximum of 101) where 
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3, 5. 1.3 (Continued) 

the region between adjoining increment dividing planes is called a mass 
addition region, Figure 3-14. Spacing of the increment dividing planes 
IS controlled by an input parameter DELZ. The increment dividing 
planes are located at each input reference plane and spaced DELZ from 
each reference plane until either a segment or the next reference 
plane is reached. This procedure is repeated and terminates at the 
aft tangent plane of the motor case. 

The location of the initial forward and aft segment faces should ‘be 
input if applicable. These inputs are SA(I) and SB(I), respectively. 

3.5.2 Propellant Pronerties 


4rDELF^ai?d''5rrsllf !"PUt 1) Tf), 2) GAW, 3} P„ 

S5?r£r."tF*^ 

Both T§ anH r™ nozzle flow rate and thrust coefficient, respectively 

propellant gas prooerties, arravs of 
temperature, y molecular wt., and c* versus’ procure can be innut 

Jnpri'^r determined by the reconstruction module or one of the 

specific impulse modules, in which case it need not be input 

3.5.3 Burning Rate Inputs 


Propellant burning rate, RB, at any increment dividing d 1 
computed in terms of the following values at the adiacent 
dividing plane in the upstream direction: 


ane, may be 
increment 


1 . 

Static pressure 

P 

2. 

Gas velocity 

U 

3. 

Mass velocity per unit area 

G 

4. 

Distance from stagnation point 

MRB 

5. 

Burning rate 

RBHI 

6 . 

Solid propellant density 

DELF 


Witir the exception of DELF, these val 
within the program. 


ues are not input but 


are calculated 
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3,5.3 (Continued) 

Fifty-one constants are available to define the burninp rate equation. 
Only the constants that are required for the particular burning rate 
equation to be used are input; the remaining constants are not input. 
These constants are as follows: 

AKG(l) through AKG{5) 

AKU(l) through AKU( 5) 

AKR(l) through AKRC39} 

AKSL0T(1) and AKSL0T(2) 

Prior to calculating burning rate, critical values of velocity, tJCR, 
and mass velocity per unit area, OCR, are obtained as follows: 

UCR = AKU(l) + AKU(2).p^^^^^) + AKII(4) 

GCR = AKG(l) + AKG(2)-P^*^*^"^^^ + AKG(4) 

Propellant burning rate is then calculated in one of tv/o ways. 

If G is greater than or equal to GCR or U is greater than or equal 
‘to UCR, the following relationship is used: . 

RB = AKR(l) + AKR(2).P^'^'^^^^ + AKR(4) + AKR(6) 

+ AKR(8).U^’^^^®^ + AKR{10)-G^‘^’^^^^^ + AKR(12).G^^^^^^^ 

+ AKR(20)*P'*^’^*’'^^^ ^ + AKR(23)-P^’^^^^^^.G^*^^^^^^ 

+ AKR(26)/[AKR(27)-P^’^^^^^'^ + AKR{29) 

The fore-head reference length (HHR) used in caTcul ati ng HRB must be 
input if fore-head burn data arrays are input. If the value of G is 
less than GCR and the value of U is less than UCR, then erosive burning 
does not exist and the following relationship is used, 

RB = AKR(35) + AKR(36) + AKR(38) 
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3.5.:i (Continued) 

simultaneously 

JreS tn^ 7 pJn ^"rP!'’ 3 .rate model one, of these values must always 
constants^m^ setting the values of ‘ 

AKUfir^AKuf? constants 

raiit AKU{4) equal to zero. The program will automati- 

cally stop If any of the terms AKG(l), AKG(2) or AKG(4) are not eaual 

equal to” 5 a?o°'''' 7 n'' ,*'’h ''.I!*'*’ ai-s n3t ’ 

that th» rrJ .’.So""'''' Program will print a statement 

that the GCR or UCR coefficients are invalid. If AKR(26) is not 

si™ taSenn:?’ "5' ''“'27) nor AKR(29) caS be nigaMvfor 
orSoiam w?if restriction is excLded, the 

program will stop and print a statement of invalid AKR(27) or AKR(29). 

The burning rate at the segment slot face is calculated by: 

RBSL0T=AKSL0T{1 ) 

where P is the static pressure in the port at the segment interface. 

For anisotropic burning, arrays of burn rate coefficients, AKRTAU(I) 
versus burn distance, TAUAKR(I) are required. Then the burn rate 
equation becomes 


RB=AKRTAU(I).P^^'^^^^ 




:) 


TS Simulated only during the buildup and tailoff 
intervals, and the same arrays may be used for both intervals. One 
option exists here for burning in the fore-head and aft-head sections 

TAliA^R^n AKRTAU(I) versus 

fore-head and aft-head sections. 

Ts desired that anisotropic burning not be considered in the 
end sections, AKRH and AKRN are set equal to AKR(2) and AKR(36). 

frequently a function of the distance burned 
radially for a given pressure. This may be due to a significant 
temperature gradient through the grain or' variations in the composition 
of the propellant constituents in the radial direction. An example of the 
latter occurs when the predominant alignment of oxidizer particles 
varies between the bulk of the propellant and the core interface 
or the case wall. 

Variation are NAKR, 

TBLAKR(I). A burn rate adjustment factor 
(AKRADJ) is determined from the TBLAKR(I) array as a function of the 
time array TIliAKR(I). The time points in the TBLAKR(I) array correspond 
0 radical burn distances for which burn rate variation data points 
are given. NAKEND determines whether or not AKRADJ will be applied in 
the end sections. 
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3.5.4 Nozzle Configuration 

The nozzle corifiguration , Figure 3’"23, is i‘np!^t fpr fv/o purposes, One, 
to calculate throat area which in turn is uspd t^f determine nozzle flow 
rate. The other is to provide the parai^ieters repfiirefi to calculate 
the nozzle thrust coefficient. ' 

Required inputs are DT, DE, and ANN (if Qther^t*^^'^ Also, one 
option in each of the following ^etp must be input; 

Ambient Pressure ; Input as a constant (PA) Qn as a variable with 
respect to time (NPA, T1!MEPA(I), and TB4PA(|)), 

I 

Nozzle Efficiency ; Input as a constant (CM) or ag a variable vn'th 
respect to time (NCH, XIt^^CM{I). , and TBLCM{I)), 

Nozzle Erosion; Input as a variable with respect to timg (NOT, TDELDT(I), 
and TJMEDTil)) or as a variable wUh respgct to np^zle ^tagnat)on 
pressure (ERBAR, P0NBAR,’ EREXP), If PP no;?z]q' enpslpt) i? desired, then 
no inputs are required. i 


Nozzle Contour ; Three possibilities exjst for describing the nozzle 
contour, The nozzle contpur must b® described ?p that a diyergejice 
flow loss factor can be determined for use in the thrust ealeulatipn. 

(i) For a conical nozzle, the nozzle cone halfrAnglP, AN2, is required. 
Also, ihe contoured nozzle indicator, JiC0r^T, must be input as zero. 


(2) For e contoured or bell-shaped noi^zle, thg nozzle half-angle can 
be approximated as the average of the initial and final nozzle 
angles. These inputs are AN2 and fHETEX, respectively, NC0NT 
must be input' as 1 . 


(3) If data is available for theoretical specific impulse fpr a 

and contoured nozzle, a more precise determination for the effective 
contoured nozzle halfr angle iS possible, This type of data is 
available from the Axi symmetric TworPhase Perfppt Gas Performance 
Program, ATPAP (Reference ^). This program can be used to calculate 
the'ratip of the theoretical specific impulse for the actual nozzle 
contour to the theoretical specific impulse for an equivalent conical 
nozzle, This ratio is called ^IRATf) and ^s inpuj along with the 
equivel^nt conical po^xle halfrangle, AN? ^nd f|C0NT?lf 

The form pf the delivered ppeffiplenii, CF, js; 


CFC 


2(GAMA) 

GAHA-f 


Z 1 


GAfiA+1 


■] 


/ GAW+T \ 
\ GAMAtI j 


1 
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3.5.4 (Continued) 


CF 


1 



/•GAMA-1 \ 

2 

. CFG 

r PE A GAMA J 

|_ PON^ 

1 

J [^.5+.5 C0S(AN2)J 




CM 


''PUn^ 


EPl 


+ 


where PE is the exit plane static pressure, PA is the ambient pressure, 
and EPl is the exit cone expansion ratio. 


3.5.5 Nozzle Submergence 

Solid propellant rocket motors are commonly restricted in total length 
for numerous reasons. Inserting or submerging the nozzle some distance 
into the motor chamber may be desirable for use on motors with a length 
restriction. This nozzle submergence permits a longer motor chamber 
with more propellant. 

Nozzle submergence introduces performance losses. These losses are 
primarily due to entrance effects, two phase flow effects, and flow 
field effects. Nozzle submergence also introduces additional con- 
siderations for the internal ballistics simulation. As shov^n in 
Figure 3-24, a portion of the propellant grain is located aft of the 
nozzle entrance plane and requires a reversed flow region. which allows 
the mass generated in the submerged area to exit the motor. In order 
to correctly simulate this phenomena, a different convergence scheme 
for the submerged nozzle case is required in the Internal Ballistics 
Module. This scheme requires the simulation of the stagnation regions 
and the flow field for the reversed flow. 

Required inputs for an internal ballistics simulation with a submerged 
nozzle are listed below. These input variables are defined in detail 
in Table 3-XVIII and several of these variables are shown in Figure 3-24 

(1) Nozzle submergence indicator, NSUBMG. Always input as one for a 
submerged nozzle simulation. 

(2) Submergence distance from aft tangent plane, HSUBMG. 

(3) Length of the cylindrical* portion of the aft dome, OELH- t^ot 
required for a nozzle that is only slightly submerged into the 
aft dome . 

(4) Maximum burn distances in Regions A and B, TAUMMA and TAUMNB, 
respectively. 

(5) Nozzle entrance coordinate data; TRSUB(I), TXSUB(I), and NPSUB. 

(6) Region B port area data; TA1B{I), TA2B(I), and TAUNB(I). 
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3,5.5 (Continued) 

be required These dahe aro hnv^n ^ odietry, additienal input data may 

?omp1efe lisVj?'’th’""* “lumr?o?1eg[uL""A""d ^"^fhe 

uZlTnslTn 5lta?r?n'?rble'U?r^= 

TAUN(I)?°'''"® functions of the independent array 

0) Burn area for each region, ABNA(I) and ABNB(I) 

( 2 ) 


(3) 

(4) 


S1Sb[;j: Tq^;f,°5 pM^iNAd) and 

gjenfK 


The following vari^bl^s are not fMnj^don? qf ptber y§riab)es' 

'\iSS“ =:"""?.S."S ,?:,“!;■ ■" 

3.5.6 Internal Gas Dynamics 

t 

Ilthe?‘st8aL’nr **'5 P^P^llant grain can be based on 

f£e|pnL"srt^sfiS,£ifii"^ 

The compulatla";i:i®]„“L?;al ij discu^sS' baip|' 

i“«r:ys:'s,s:",;‘.ss :: sy;.srs» 
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3.5.6 (Continued) 

terns to the momentum and continuity equations enables the program to 
be used for calculating ignition transients. 


Steady Flow 


In order to simulate steady flow internal gas dynamics, STDYST is 
set to 1.0 and STFLAG is set to 0.0. An initial guess for the fore- 

is required. The necessary time control inputs 
are DELTSS, DELTT0, and ANITW. These are discussed in Section 3.5.7. 


Non-Steady Flow 

In order to simulate non-steady flow, STFLAG is set to 1.0, and STDYST 
is set to 0.0. The initial guess for the fore-head pressure, PHI, 

IS input, however, unlike the case of steady flow, PHI must be input 
as twice the value desired. For example, if a start transient is 

with an ambient value of 14.7 for the fore-head pressure, then 
PHI=29.4. 


When exercising the option to input a buildup pressure trace and thus 
compute a" (in r=apb), the trace is input as the pressure array, 

PHST(I), versus the time array, TIMEPH(I). NPH is also input to specify 
the number of points in these arrays and PHI should be input as twice 
the value of PHST(l). AINCPL is the desired increment plane location 
at which the constant "a" is to be evaluated. 

The necessary time control inputs are DELTST, DELTSS, DELTT0, ANITW, 
and one of the start transient termination point inputs, TST, PST, or 
AITST. All of these time control inputs are discussed below. 

At the end of the start transient, the value of the computed constant 
a , or AKR(2) and AKR(36) of the input burn rate equation, may be used 
for the remaining calculations. If RBFLAG=0.0, AKR(2) and AKR(36) will be 
used during the web- time and tail -off calculations. If RBFLAG=1.0, then 
the last value of "a" will be used in place of AKR(2) and AKR(36) during 
the web- time and tail off calculations. 

3.5.7 Time Control Inputs 

Several options are available to define the calculation time interval. 
Output from the internal ballistics module is provided after each 
calculation cycle, therefore, time interval input d^ta simultaneously 
control module calculation and printout intervals. 

The time increments used during the start transient, steady-state or web- 
time, and tailoff intervals are DELTST, DELTSS, and DELTT0, respectively. 
If a start transient is simulated, DELTST must be non-zero. The start 
transient interval can be terminated by either inputs of time, TST, fore- 
head pressure, PST, or total impulse, AITST. The steady-state time 
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increment DELTSS is used during the web-time interval and is used until 
the tail-off interval begins. Tail -off starts when v/eb burnout occurs 
at a number of increment dividing planes which equals or exceeds ANITW. 

If ANITW is not input, a value of 1.0 is automatically assumed. The 
tail -off increment DELTT0 is used until run termination occurs. 

A simulation is normally terminated when the pressure ratio across 
the nozzle is insufficient to maintain sonic flow at the throat. 

The program can also be terminated earlier by the TIMAX flag, i.e., - 
v/hen time reaches or exceeds TIMAX. 

An option exists by which a special time increment may be used. This 
increment is DELTSP and will override any other program time increment 
betv/een the input time valuer pf TIMPJl and TIMpT2, 

3.5.8 Propellant Mass Correlation 

The Internal .Ballistics Module determines initial propellant mass from 
the inputs for propellant density and the propellant grain dimensions. 
Tolerances characteristic in measuring propellant dimensions and determining 
propellant density may be such that the mass calculated with these inputs 
may vary islightly from the mass determined by actually weighing a pro- 
pellant grain. These inaccuracies can be compensated for with a program 
option called propellant mass correlation. Figure 3-7 presents the 
■ logic flowchart for this correlation which- is performed in subroutine 
SETPH. 

Four possibilities are available for the program user. 1) Decide to 
make no correlation; 2) adjust the propellant mass by biasing the average 
propellant density, 3) adjust the propellant mass by biasing the grain 
dimensions, or 4) adjust the propellant mass by biasing the propellant 
density and the grain dimensions. Density adjustments are made in- 
ternally but grain dimension adjustments mus-^ be made through the 
input data. 

Related inputs are 1 ) WFACT, 2) RH0T0L, and 3) WFT0L, WFACT is the actual 
fuel weight which would be knov/n from grain v/eighing measurements. If 
WFACT is not input, no propellant mass correlation is performed. RH0T0L 
is the maximum amount allowed for a density adjustment. WFTOL is the 
minimum difference between WFACT and the calculated fuel v/eight which 
will be considered to be significant enough to perform a correlation. 




TALLE 3-XVII VARIADl.E DEFIMITIONS FOR flOTOR CASE AND ’’ROPELLAflT 
GRAIN GEOMETRY FACTORS 


VARIABLE 

AINCIN(I) 

AK 

AKK 

ALA(I) 

ALB(I) 

ALE(I) 

ALSl(I) 

ALS2{I) 

AN0(I) 

A0HM 

A0NM 

A01(I) 

A02(I) 


DEFINITION 


Reference plane location measured from forward 
tangent plane. See Figure 3-14. (IN) 

Adjusting factor w.hich determines the distance 
betvieen planes for the head-end web analysis and the 
end section analysis. A value of 0.5 is recommended. 
Required only when simulating a head-end with v;eb. 

Adjusting factor used in the head-end web block 2 
analysis to determine distance between planes. A 
value of 0.5 is recommended. Required onlv v/hen 
simulating a head-end with v/eb. 

Perimeter length of sector ALA for a reference plane. 
See Figures 3-15, 3-18 and ..-19. (IN) 

Perimeter lengt:h of sector ALB for a reference plane. 
See Figure 3-15. (IN) 


Perimeter length of sector ALE for a reference plane. 
See Figure 3-15. (IN) 

Length of web to inner grain point of a reference 
plane. See Figures 3-15, 3-19 and 3-20. (IN) 

Length of web to outer grain point of a reference 
plane. See Figure 3-15. (IN) 

Number of grain cross-sectional symmetrical parts 
for a reference plane. This works out to be half 
the total number of "pie slices" as shov^n in 
Figures 3-15 and 3-19 through 3-22. 


Maximum allowable angle between tangent to fore- 
head ellipse and motor centerline. See Figure 
3-11. (DEGREES) 

Maximum allowable angle betv/een tangent to fore- 
head ellipse and motor centerline. See Figure 
3-11. (DEGREES) 

Reference plane angle of side ALA. See Figures 
3-15 and 3-19 through 3-21. (DEGREES) 

Reference plane angle of side ALB. See Figures 
3-15 and 3-19 through 3-21. (DEGREES) 
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TABLE 3-XVII 

VARIABLE DEFINITIONS FOR MOTOR CASE AND PROPELLANT 
GRAIN GEOMETRY FACTORS (Continued) 

VARIABLE 

DEFINITION 



A03(I) 

Reference plane angle of side 
3-15 and 3-19 through 3-21. 

ALC. See Figure 
(DEGREES) 

A04(I) 

Reference plane angle of side 
3-15. (DEGREES) ‘ 

ALD. 

See Figure 

A05(I) 

Reference plane angle of side 
3-15. (DEGREES) 

ALE. 

See Figure 

BH 

Fore-head case ellipse ratio, 
be either 

The value of BH must 


a. Equal to 1.00 (hemispherical end) 

b. Equal to or greater than 1.01 (oblate spheroid 
end) 

c. Equal to or less than 0.99 (prolate spheroid 

end) < 

See Figure 3-12. 


BN 

BTA0E 


Aft-head case ellipse ratio. See BH definition 
above and Figure 3-11. 

Ratio of (head-end web ellipse axis normal to motor 
axis) to (ellipse axis parallel with motor axis). 
See Figure 3-12. Required only v'hen simulatino a 
head-end with web. 


DELZ 

DHl 

DLRF 


DNl 

DRVRF 


Maximum length of mass addition region. See Figure 


Fore-head case opening diameter. See Figure 3-11. (IN) 

aL/RF, where aL = increment size used for the burning 
surface area calculations, measured along the internal 
perimeter in the adjacent reference plane. A value 
of 0.01 is recommended. 

Aft-head case oppniog diameter. See Figure 3-11. (IN) 

aRv/RF, where aRv = thickness of each head-end with 
v>/eb volume increment, and RF = outside propellant 
radius of adjacent or first reference plane. A value 
of 0.01 is recommended. Required only when simulating 
a head-epd v/ith v/e)3. 
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TABLE 3-XVII 

VARIABLE DEFINITIONS FOR MOTOR CASE AND PROPELLANT 
GRAIN GEOMETRY FACTORS (Continued) 

VARIABLE 

DEFINITION 

DTAU(I) 

Increment used in generating geometry tables as a 
function of distance burned. This increment is 
used until the web has reached the wall, then 
DTAUW(I) is used. See definition of DTAUW{I) 
below. (IN) 

DTAUW(I) 

Increment used beyond web thickness (TAUW) in 
generating geometry tables as a function of 
distance burned. DTAU(I) and DTAUW(I) must be 
chosen such that the total number of points in 
the geometry tables does not exceed 50. (IN) 

HC0 

Length of motor's cylindrical section. See 
Figure 3-8. (IN) 

RF(I) 

Outer radius of propellant for a reference plane. 
See Figures 3-15 and 3-19 through 3-22. (IN) 

RIG 

Igniter opening radius. See Figure 3-12. (IN) 

R2(I) 

Radius of inner grain point for a reference olane. 
See Figures 3-15 and 3-19 through 3-21. (IN) 

R3(-I) 

Fillet radius between ALA and ALB for a reference 
plane. See Figures 3-15 and 3-19. (IN) 

R4(I) 

Fillet radius betv/een ALB and ALC for a reference 
plane. See Figure 3-15. (IN) 

R5(l) 

Fillet radius betv^een ALC and the web for a reference 
Diane. See Figures 3-15 and 3-19 through 3-21. (If!) 

R6(I) 

Fillet radius between the web and ALD for a reference 
plane. See Figure 3-15. (IN) 

R7(I) 

Fillet radius betvieen ALD and ALE for a reference 
plane. See Figure 3-15. (IN) 

R8{I) 

Fillet radius between forks of forked wagon wheel 
for a reference plane. See Figure 3-15. (IN) 

SA(I) 

Location of slot forv'/ard interface measured from 
forward tangent plane. See Figure 3-14. (IN) 
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TABLE 3-X.VII 


VARIABLE DEFINITIONS FOR MOTOR CASE AND PROPELLANT 
GRAIN GEOMETRY FACTORS (Continued) 


VARIABLE 

SB(I) 

TAUW(I) 

TH0(I) 

TSLVR(I) 


DEFINITION 

Location of slot aft interface measured from forward 
tangent plane. See Figure 3-14. (IN) 

Reference plane web thickness. See Fioures 
3-15 and 3-19 through 3-22. (IN) 

Reference plane central angle for a slotted-cone 
configuration. See Figure 3-21. (DEGREES) 

Distance from core interface to inert sliver for 
a reference plane. Input not required unless an 
inert sliver is to be simulated. See Figures 
3-19 through 3-21. (JN) 
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TABLE 3-XVIII 

VARIABLE 

ANN 

AN2 


CM 

DE 

DELH 


DT 

ERBAR 


EREXP 


HSUBMG 

NCM 

NCONT 


VARIABLE DEFINITIONS FOR NOZZLE PARAMETERS 

DEFINITION 

Number of nozzles. Assumed to be 1.0 if not 
input. 

Nozzle half-angle. Required for a conical 
nozzle (when NC0NT=O) and for a contoured or 
bell nozzle (when NC0NT=1). See Figure 3-23. 
(DEGREES) 

Nozzle efficiency. Not needed if the TBLCM(I) 
array is input. 

Nozzle exit diameter. (IN) 

Length of the cylindrical portion of the aft 
dome. Assumed to be 0.0 if not input. In 
simulating a submerged nozzle, the nozzle 
entrance must remain within the aft-dome 
and DELH enables the aft-dome region to be 
expanded if needed. When DELH is input, HC0 
should be reduced by the same amount. See 
Figure 3-24. (IN) 

Nozzle throat diameter. (IN) 

Average erosion rate of the nozzle throat 
radius. Used only wh£n the empirical erosion 
rate equation (£=E(P/p)N) is simulated. Also, 
see EREXP and PPNBAR. (IN/SEC) 

Nozzle erosion rate exponent. Used only 
when the empirical erosion rate equation 
(£=E’(P/F)^) is simulated. Generally input 
as 1.0 when used and when more definitive 
data is not available. Also, see ERBAR and 
P0N8AR. 

Distance from the aft tangent plane to the 
nozzle inlet. Used only when a submerged 
nozzle is simulated. See Figure 3-24. (IN) 

Numbar of points in the nozzle efficiency versus 
time arrays. Maximum allowable=50. Also, see 
TBLCM(I) and TIMECM(I). 

Nozzle contour indicator. Input as zero for 
a conical nozzle and as one for a contoured 
or bell-shaped nozzle. Also, see AN2, THETEX, 
SIRATC, and Figure 3-23. 
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TABLE 3-XVIII 

VARIABLE 

NOT 

MPA 

NPSUB 

NSUBMG 


PA 

PPNBAR 


SIRATC 


TA1B{I) 


TA2B(I) 


VARIABLE DEFINITIONS FOR NOZZLe PARAMETERS (Continued) 

DEFINITION 


Number of points in the nozzle throat diameter 
variation versus time arrays. Maximum allow- 
able=25. These arrays will not be used if 
the P0NBAR simulation is used. See TDELDT(I) 
and TIMEDT(I) below. 

Number of points in the ambient pressure versus 
time arrays. Maximum alloW',ble=50. Also, see 
TBLPA(I) and TIMEPA(I). 

Number of points in the TRSUB(I) and TXSUB(I) 
nozzle coordinate arrays. Maximum number=20. 

Submerged nozzle indicator. For no submergence 
simulation, a value of zero is assumed. For 
simulating a submerged nozzle, a value of one 
should be input. 

Ambient pressure, Not required if the TBLPA(I) 
array is input. (PSIA) 

Average nozzle stagnation pressure expected 
during a burn. Used only when the empirical 
erosion rate equation (e=t(P/P)N) is simulated. 
If not input, P0NBAR is automatically given a 
0.0 value and this erosion rate simulation 
will not be used. Also, see ERBAR and EREXP, 
(PSIA) 

Ratio of the theoretical specific impulse for 
the actual nozzle contour to the theoretical 
specific impulse for an equivalent conical 
nozzle. See "Nozzle Contour" in Section 3.5.4, 

Table of part area at the forward end of 
region B in a submerged nozzle simulation. 

Input as: a function of TAUNB(I) and calculated 
as area enclosed by propellant grain minus the 
nozzle entrance area. (IN/'*2) 

Table of projected port area in the aft end 
Of region 8 in a submerged nozzle simulation. 
Input as a function of TAIJNB(I) and caTculated 
as the port area at tlie aft-most propellant 
point in region B minus the external area of 
the nozzle at the submergence point in the 
case. (IN**2) 
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TABLE 3-XVIII 

VARIABLE 

TAUMNA 

TAUMNB 

TAUNB(i) 

TBLCM(I) 

TBLPA(I) 

TDELDT(I) 

THETEX 

TIMEGM(I) 

TIMEDT(I) 

TIMEPA(I) 

TRSUB(I) 


VARIABLE DEFINITIONS FOR NOZZLE PARAMETERS {Continued). 

DEFINITION 

Maximum burn distance in region A. Must be 
input if a submerged nozzle is simulated. 

(IN) 

Maximum burn distance in region B. Must be 
input if a submerged nozzle is simulated, 

(IN) 

Independent array of burn distance used with 
the dependent -arrays TAIB(I) and TA2B(I), 

Maximum number of points=50. (IN) 

Array for nozzle efficiency as a function 
of time. Also, see NCM and TIMECM(I). 

Array for ambient pressure as a function of 
time. Also, see NPA and TIMEPA(I). (PSIA) 

Array for change in nozzle throat diameter. 

Also, see NDT and TIMEDT(I). (IN) 

Nozzle expansion skirt half-angle measured 
at the nozzle exit. See Figure 3-23 and 
"Nozzle Contour" -in. Section 3.5.4. 

(degrees) 

Time array for nozzle efficiency. Also, see 
NCM and TBLGM(I). (SEC) 

Time array for nozzle throat diameter change. 
Also, see NDT and TDELDT(l). (SEC) 

Time array for ambient pressure. Also, see 
NPA and TBLPA(I). (SEC) 

Table of nozzle radii input as a function 
of TXSUB(i) and used to calculate submerged 
nozzle volume and nozzle contraction ratio. 

First point in table must correspond to TXSUB(l) 
0.0 and must be the nozzle entrance radius. 

Last point in table must correspond to TXSUB 
(NPSUB) and must be the nozzle radius at the 
point of submergence. Radii are measured from 
nozzle centerline to nozzle external surface. 

See Figure 3-24. (IN) 


TABLE 3-XVIII 

VARIABLE 

TXSUB(I) 
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VARIABLE DEFINITIONS FOR NOZZLE PARAMETERS (Continued) 

DEFINITION 

Table of X coorcji nates for the submerged nozzle 
radii table, TRSUB{ I). TXSUB( I) is measured 
from the nozzle entrance plane. TXSUB(l) 
must be equal fo 0 . 0 ; TXSUB(NPSUB) must be the 
X coordinate of the submergence point of the 
nozzle in the case referenced to the nozzle 
entrance plane. See Figure 3-24. (IN) 
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TABLE 3-XIX VARIABLE DEFINITIONS FOR GEOMETRY TABLES PARAMETERS 

(If all geometry tables are to be calculated by the internal ballistics 
module, then the terms marked * are not required) 


VARIABLE 


DEFINITION 


ABHD(I)* 

ABN(I)* 


ABNA(I)* 


ABNB(I)* 

AKGYP(I,N)* 

ALPPL(I,N)* 

AP0RT(I}* 

GE0C0N(I,N)* 


KPLANE 

NGE0(I) 

NGE0HD* 


Dependent array of head-end burn area as a function 
of burn distance. Also, see TAUHD(I) and NGE0HD. 
(IN**2) 

Dependent array of aft-end burn area as a function 
of burn distance. Not needed if a submerged nozzle 
is simulated. Also, see TAUN(I) and NGE0MN. 

(IN**2) 

Dependent array for burn area of region A of the 
aft-end as a. function of burn distance. Used 
only when nozzle submergence is simulated. Also, 
see TAUN(I) and NGE0MN. (IN**2) 

Same as ABNA(I) except for region B. (IN**2) 

Dependent array of radius of gyration for reference 
plane N as a function of burn distance. Also, see 
the independent array TAUPL(I,N) and NGE0{I). (IN) 

Dependent array of port perimeter for reference 
plane N as a function of burn distance. Also, see 
the independent array ’TAIJPL(I ,N) and NGE0U). (IN) 

Initial port area for a reference plane. (IN**2) 

Reference plane geometry constants. There are 45 
constants per reference plane, see Figures 3-16 and 
3-17. These inputs are seldom, if ever, hand 
calculated for a run. They can be obtained from 
the punch card output of a previous internal 
ballistics run. 

Number of reference planes. Must always 
be input. 

Number of points in the reference plane geometry 
arrays. Maximum allowable =50. See TAUPL(I,N), 
ALPPL(I,N), and AKGYP(I,N). 

Number of points in the head-end geometry arrays. 
Maximum allowable = 50. See ABHD(I), TAUHD(I), 
PM0IHD(I), RM0IDH(I), and XCGHD(I). 


TABLE 3-XIX 

variable 

NGE0MN* 


NPUTAB 


NWRTAB 

PM0IHD(I)* 

PM0If'l(I)* 

PM0INA(I)* 

PM0INB(I)* 

RM0IHD(I)* 

RMOINA(I)* 

RM0INB(I)* 

RM0IN(I)* 


0256-10020-2 

Sln^df tables parameters 

DEFINITION 

Number of points in the aft-end geometry arrav^ 
Maximum allowable = 50. See ABNIIl 

Geometry array punch flag. 

- l’ Pnn^ff constants output 

: J’ and plane constants^ 

" ’ constants on disk 

- Ti With NAMELIST data for multiple case use 

' ’ wriS “"=tants and 

Geometry array print flag. 

“ ?’ Do^no^ print gpometry arrays. 

= la Prin| geometry arrays in table form. 

Dependent array of head-end roll M.0,1. about 
(SLUG-ir* 2 ) TAUHD(I) and NGE0HD. 

Dependent array of aft-end roll M.O.I. about 
longitudinal axis. Not needed if a submerged nozzle 
(SLUG-IN**2)‘ see.TAUN(I) and NGE0MN. 

Dependent array of aft-end region A roll M.O I 
about longitudinal axis. Used only when nozzle 

submergence is simulated. Also, see TAUNfl) anH 
NGE0MN. (SLUG-IN**2) lAUNUl and 

Same as PM0INA( I) except for region B. (SLUG-IN**2) 

Dependent ahray of head-end pitch M.O.I. about aft 

(SLUG-IN**2f ‘ see TAUHD(I) and NGE0HD. 

Dependent ai'ray of aft-^end region A pitch M.O.I. 
about aft tangent plane. Used only when nozzi®' 
submergence is simulated. Also, see TAUNfll and 
NGE0MN. (SLUG-IN**2) and 

Same as RM0INA{l) except for regioji B. (SLUG-IN**2) 

Dependent array of aft-end pitch M.O.I. about aft 
tangent plane. Not needed if a submerged nozzle is 
simulated. Also, see TAUN{I) and NGE0MN. (SLUG-IN**2) 
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TABLE 3-XIX 

VARIABLE 

TAUHD(I)* 

TAUN(I)* 

TAUPL(I,N)* 

VCHINP* 

VCMINP* 

VCN'IMA* 

VCNINB* 

VFHp* 

V'FN(3* 

VFN^A^' 

VFN0B* 

XCGHD(I)* 

XCGN(I)* 


XCGNA(I)* 


XCGNB(r)* 


VARIABLE DEFINITIONS FOR GEOMETRY TABLES PARAMETERS 
(Continued) 


DEFINITION 

Distance burned; the independent array for 
the head-end dependent arrays ABHD(I), TAUHD(I), 
PH0IHD(I), RM0IHD(I), and XCGHD(I). (IN) 


Distance burned; the independent array for the 
aft-end dependent arrays ABN(I), ABNA(I), 
ABNB(I), PM0IN(I), PM0INA(I), PM0INB(I), 
RM0IN(I), RM0INA(I), RM0INB(I), XCGN(I), 
XCGNA(I), and XCGNB(I). (IN) 


Distance burned; the independent array for the 
reference plane dependent arrays ALPPL(I,N) 
and AKGYP(I,N). (IN) 

Volume of fore-head case. (IN**3) 

Volume of aft-head case. (IN**3) 


Volume of aft-head case. Region A. (IN**3) 


Volume of aft-head case. Region B. (IN**3) 

Initial propellant volume in fore-head. (TN**3) 

Initial propellant volume in aft-head. (IN**3) 

initial propellant volume in region A of the 
aft-end. Used only when nozzle submergence is 
simulated. (IN**3) 


Same as VFN0A except for region B. (IN**3) 

Dependent array of head-end propellant C.G. 
location measured from the forward tangent plane. 
Also, see TAL'HD(I) and NGE0HD, (IN) 

Dependent array of aft-end propellant C.G. 
location measured from the aft tangent plane 
(always positive). Not needed if a submerged 
nozzle is simulated. Also, see TAUN(I) and 
NGE0MN. (IN) 

Dependent array aft-end region A oropelTant C.G. 
location measured from the aft tangent plane (always 
positive). Used only when nozzle submergence is 
simulated. Also, see TAUN(I) and NGE0MN. (IN) 

Same as XCGNA(I) except for region B. (IN) 
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TABLE 3-XX VARIABLE DEFIMITIONS FOR PROPELLANT PROPERTIES AMD 
BURNING RATE MODEL PARAMETERS 


(** - See Section 3.5.3 for discussion of burnina rate inputs) 


VARIABLE 

AINCPL 

AKG(I)** 

AKR(I)** 

AKRH** 


AKRN** 


AKRTAU(I)** 

AKSL0T(I)** 

AKU(I)** 

AMWG(I) 

CSTAR 

CSTR(I) 


DELF 

GAMA 


DEFINITION 

Inprement dividing plane location in cylindrical 
section where anisotropic burning is to be evaluated. 
(IN) 

Critical mass velocity ^GCR) per unit area equation 
constants. 

Burning rate equation constants. 

Head-end burp rate coefficient fpr apisptroDic 
propellant burping simulation. An input vaipe 
for AKRH will override any table value for 
anisotropic burn rate coefficient. 

Aft-end burn rate coefficient for anisotropic 
propellant burning simulation. An input value for 
AKRN will override any table value for anisotropic 
burn rate coefficient. 

Array of anisotropic burn rate coefficients. This 
is the dependent array; see also the independent 
array of distance burned, TAIIAKR(I), and NAKRST. 

Slot interface burning rate equation constants. 

Critical gas velocity (UCR) equation constants. 

Array of propellant gas molecular weight. This is 
a dependent array; see also the independent array 
of total pressure, PRES3(p, and NpSTH. (LPM/MOLE) 

Propellant gas characteristic exhaust velocity. 

Not required if the CSTR(I) arrav is input. 

(FT/SEC) V 

Array of propellant gas characteristic exhaust 
velocity. This is p dependent array; see also the 
independent array of total pressure, PRESS(I), and 
NCSTR. (FT/SEC) 

Solid propellant density. (LBM/IN3) 

Propellant gas specific heat ratio. Not required 
if the GAMAG{ I) array is input. 
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TABLE 3-XX 
VARIABLE 

gamag(i) 

HHR 


NAKENO 


NAKR 

NAKRST 

NCSC0E 


NCSTR • 


NPH 


definitions FOR propellant PROPERTIES and 
BURNING rate MODEL PARAMETERS (CONTINUED) 


DEFINITION 

Array of propellant gas specific heat ratios. This 
? dependent array; see also the independent arrav 
of total pressure, PRESS(I), and NCSTR. 


Reference length of fore-head section. Distance 
along motor centerline from forward tangent 
plane to the propellant gas stagnation point. 
Required only when burn area tables are input 
for the head-end. (IN) 


Indicator which specifies whether or not the burn 
rate adjustment factor will be used for the end 
sections. 

= 0, neither end 
= 1 , head-end only 
= 2, aft-end only 
= 3, both ends 

Also, see TBLAKR(I) and TIMAKR(I). 


Number of points in the burn rate adjustment 
time arrays. Maximum allowable = 25. Also 
TBLAKR{I) and TIMAKR(I)> 


versus 

see 


Number of points in the arrays of anisotropic burn 
CnJ'n coefficient, AKRTAU(I), versus burn distance, 
TAUAKRd). Maximum number of points = 30 . 


Indicator to suppress the c'' 
fit calculation. If NCSC0E 


second order curve 
is greater than zero. 


the second order c* calculation' (CSTAR = CSC0EF(1) + 
CSC0EF(2)*P0N + CSC0EF(3)*P0N**2) is bypassed. 

ThuSj "th6 1 ntsrnsl bsllistics moduls will US6 its 
own input value for c* instead of that determined 
by the reconstruction module or one of the specific 
impulse modules. 


Number of points in tHe propellant gas property 

AMWG{I),'gAMA(I), and 

TC0MB(I). If NCSTR is input, then all five arrays 
must be input and each must contain the same no.* 
of points. Maximum no. of points = 20. 

Number of points in the arrays of fore-head pressure, 
PHST(I)., versus time TIMEPH(I). Maximum no, of points=70, 
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TABLE 3-XX 

VARIABLE 

NTAUT0 

PHST(I) 

PRESS(I) 

R 

RBFLAG 

SLTBpN(I) 

TAUAKR(I)** 

TBLAKR(I)** 

TCj3MB(I) 


VARIABLE DEFINITIONS FDR PROPELLANT PROPERTIES ANP 
BURNING RATE MODEL PARAMETERS (Continued) 

DEFINITION 

Indicator used to suppress anisotropic tail off 
calcula1;ions. If not input, a value of 1 is assumed 
and any anisotropic tailoff simulation which might 
otherwise be simulated will be suppressed (i.e., 
an isotropic tailoff will be simulated). If input 
as 0, an anisotropic tailoff will be unaffected. 

This enables the userto simulate a start transient 
by specifying the buildup pressure trace (and thus 
specifying an ani^otrgpic buildup) without having 
a corresponding anisotropic tailoff. 

Array of head-end total pressure during the start 
transien);. This i^ the dependent array; see also the 

independent array, TIMEPH (I), and NPH. (PSIA) 

Independent array of total pressure for the dependent 
propellant gas property arrays, ,QSTR(I), AHWR(I), 
TCPPCI), and RAMAG(i). (psia) • 

Propellant gas constant. (FT/°R) 

Program control flag to set the burn rate equation 
coeffieien'^5 AKR(^) and ANR(,3^) edVai to tM aniso- 
tropic burn rate coefficient AKRST at the termination 
of the start transient Interval, A nonr-zero value 
will revise the burn rate coefficients. 

Indicator to control total inhibiting of burning 
on the forward and aft interface of a slot, I 
indicates the slot number; 

= 0^0, neither interface inhibited 
= 1 .p, forward interface inhibited only 
= 2.0, aft interface inhibited only 
- 3,0, bpth interfa'ces inhibited 

Independent array of burp distance for anisotropic 
b^jjning simulation. See else A.KRTAPU) and NAKRST. 

Dependent array for radial I'ate variation as 
a function of time, Also, see NAKR, NAKEND, TIMAKR(I), 
and Section 3.5.3. 

Array of head-end total gas temperature. This is 
a dependent array; see also the independent array of 
total pressure, PRESS(I), and NpSTR. (“R) 
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TABLE 3-XX 

VARIABLE 

TIMAKR(I)** 

TIMEPHd) 

Tfi 


VARIABLE DEFINITIONS FOR PROPELLANT PROPERTIES AND • 
BURNING RATE MODEL PARAMETERS (Continued) 

DEFINITION 

. Independent array of time for radial burn rate 
variationi Also, eee NAKR, TIMAKR(I), and Section 
■ 3.5.3. (SEC) 


Array of time durincj the start transient. This is 
the independent array; see also the dependent array 
for buildup pressure, PHST(I), and NPH. (SEC) 


Head-end combustion gas total temperature. Not 
required if the TC0HB(l) array is input. (°R) 

I 
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TABL^ S^XXI VARIABLE pEPINITIQNS FOB NQN-STEADY FLOW, STEApY-STATE, 
AND PROGRAM TIME CONTROL PARAMETERS 

VARIABLE DEFINITION 


AITST 

ANITW 

DELTSS 

DELTST 

DELTT0 

belts? 

PST 

STDYST 

STFLAG 

TIMAX 


Desired raaxlmum valpe of total impulse during the 
start transient. If AITST is input, the start 
transient interval will terminate v/hen the total 
delivered impulse, AIT, reaches or exceeds AITVAC. 
Also see PST and TST and pote that only one of these 
is required to terminate a start transient. (LBF- 
SEC) 

Number of inerpment dividing planes at which web 
burn through must opcur’ before tail off begins. 
Assumed to be 1,0 if not input. 


Computation time increment used during steady-state. 

I udUJi 


Qo(nputat1.on time increment used during the start 
transient. Not required if STFLAG is mt input. 
(SEe) 


Computation time increment used during tail off. 

(SEC) ' ' ‘ 

Special CDrnpptatiprv time increment. Can be 
applied at thq discretion of the program user 
between the input time points TIHpTl and TIMPT2. 

(SEC) ^ 

Desired maximum- value of total pressure during the 
start transient. If PST is input, the start transient 
interval will terminate v/hen the total head-end 
pressure, PH, reaches or exceeds PST. Also-, see AITST 
and TST and note that only one of these is required 
to terminate e ptsrt transient, (PSIA) 

program control flag to indicate steady^state. 
Required as input only when a start transient is 
not simulated (i,e., v/hen STFLAG js not input). 

Should bq input as 1.0 when used, 

Start transient calculation cpntro] flag. If a 
start transient is desired, spi; STFLAG t 1.0. 

Maj<imuro program time Tpr internal ballistiqs 
simulatipn,. (SEE) 
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TABLE 3-XXI 

VARIABLE 

TIMPTl 

TIMPT2 

TST 


VARIABLE DEFINITIONS FOR NON^STEADY FLOW, STEADY-STATE, 
AND PROGRAM TIME CONTROL PARAMETERS (Continued) 

DEFINITION 

If the specie! calculation time increment, DELTSP, 
is used, this is the time at which it begins to 
be used. Also, see TIMPT2 below. (SEC) 

.The time at which DELTSP will stop being used. 

If not input, a value of 0.0 is assumed and 
DELPSP will not be used. If input greater than 
0,0, then DELTSP will be used when TIMPTl < 

TIME ^TIMPTE. (SEC) “ • 

‘Desired maximum value of time during the start 
transient, If TST is input, the start transient 
interval will terminate when the program time, 

TIME-,- reaches or exceeds TST. Also, see AITST 
and PST and note .that only one of these is required 
'to terminate a start transient. (SEC) 



TABLE 3-J{XII 

VARIABLE 

ACCELT(I) 

CKDUMP(I) 

CRP 

CRT 

CRW 

EPCA'I) 

EPCN(I) 

GE0RUN 

KDUMP(I) 


D25G-rl 0020^2 


VARIABLE DEFINITIONS FOR MISCELt,AN^QyE INPUTS 

DEFINITION 

I . - - I — r— r 

Dependent array of vehicle lonpitudinai acceleration 
used in sitnulatino acceleration dependent burning 
rates. Also, see the independent array of time, 
TIMEAC(I), and MACCEL. (G's) 

Times at which the internal ballistics diagnostic 

data dumps are started and terminated. Also, see 
KDUMP(I}' and the descr'|ptiqri of "diagnostic data dumps 
in the internal ballistics part of PROGRAM OUTPUT. 


Convergence value for npn-steady flow discharge 
pressnre at exit of e.agli piass addition region^ 

Assumed to he O.OQj if not input. This is usually 

quite sufficient,^ . 

Qonvergenqe value for ponrsteady f,lov/ discharge 
gas temperature at exit of each mass addition’ region. 
Assunied to bq 0,001 if not input, This is usually 
quite sufficient, ' ' ' 

CpnMergence yalue used tq cQ[npare flov/ rates at 
exit of grain to that which can be discharged tlinouqh 
the no^xle at the same total pressure. ' Assumed to be 
0,001 if not input, This is generally sufficient 
but a smaller yalqe j'dTl give a closer qonvergenpe 
if the need for apquracy is warranted, 

dependent array fop qase strain at the forward tangent 
plane, T^ie case strain is usqd to determine a port 
area change v^hich results from internal prfi^sure. 

Also see EPCN(I), TiMEPS(I), apd MEPS. 

Dependent apray fpr case strain at the aft tangent 
Plane, Also,- see EPPAjl), T|M,EP 5 (I), and NEP 5 . 

■ Indicator used- tp. itermi nat'o an Internal Bal 1 i sti cs 
Module; simUlatibq after referep^^.e plap.e and end 
section gepmetry data h,ave l?,een calculated 'apd 
output, ’If input as r.jQ, only the geomqtfy 
simulation will be performeil. ' jf not input, a 
value of Q.p is assumed and a full Internal 
Ballistics raodule simulation will he performed, 

Indicators useh tp qxepcis.e Internal ballistics 
diagnostic data dumps,’ ^1 So/ see CKDUMPll) and the 
de;§criRt1csn of "diagnostic data dumps" in the 
internal ballistles mpduTe paft pf PROGRAM OUTPUT. 
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TABLE 3-.XXII 
VARIABLE • 
Kf1(aiCG 

NACCEL 

NEPS 

PCTAB 

PHI 

PRTFLG 

RH0T0L 
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VARIABLE DEFINITIONS FDR MISCELLANEOUS INPUTS (Continued) 

' ■ DEFINITION 

Program indicator used to suppress C.G, and M.O.I. 
calculations. If KM0IQG - 1, calculations are 
suppressed, ' -V 

• Number of points in the acceleration versus time 
arrays, A'CGELT(I) and TIMEAC(I), Maximum no. of 
points =' 50. ' / ' 

Number of points In the case strain versus time . 
arrays, EPCA(I), EPCN(I), and TIMEPS(I), ’ Maximum 
no., of points f 50. , 

Percent. of total . burn' apea,’ relative to 1.0, 
required to force ballistic solution for the aniso- 
tropic burn r,a.te^coefficient,AKRST, and head-end 
'pre'ssure, PH, at‘ each time .interval . Required 
only when it is desired to input both head-end 
"'pressure, PHST(T) , and anisotropic burn rate coefficient, 
AKRTAU(I). Seldom used option. 

Ini tiaT estimate for head-end pressure. If STDYST 
is input as 1.0, then PHI is the estimate of head- 
end pressure at TIME =0.0, If STFLAG is input as 
l .O, then oneThalf the value of PHTwill be used 
for the head-end pressure at TIME - 0,0. 

Indicator which controls output of the Internal 
Ballistics Modqle data. 

PRTFLG =0, Motor ballistic data output only. 

= 1, Motor ballistic data and increment 
dividing plane data will be output, 

= 2, Motor ballistic data and grain C.G. and 
M.O.I. data will be output. 

= 3, Motor ballistic data, increment dividing 
plane dataj and grain C.G. and M.O.I. 
data will be output, 

Maximum amount allowed for a density adjustment in 
a propellant mass correlation, The absolu'te value 
of [(adjusted density-input density)/input density] 
is not allowed to exceed RH0T0L. This prevents an 
unreasonable adjustment to density. 



1 

/ 
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?AB^ 3^HX|I 
VARIABLE 

■ ' r ' * 

TIMEAe(iI) 

TIMEPS(I) 

WFACT 

WFT0L 




VARIA84E pEflHlTW FOR MI-^PE^UN^QUS INPUTS (Continued) 

DEFINITION 

Trrrr— 1TT7— 

Independent e)^ray of time .used In s1mulat1nc| 
acceleration dependent burping rates, Also, see 
AqpELT(I) and NApCE4\ ■ {§EC) ' 

Independent array of time used in slaiulating the 
case effects on port area, Also, see NEPS, 

EPCA(I) ■ and EPCN(I); 

Known or insasured prQRplIapt ita^s toward 
a propellanij ma§s pqrrelation Is tU made. If 
not Input, h Value of p,p f? a5fuma?l and no 
Gorrelapon W.1T1 be perfqrjrpd, (tBM') 

The mlnimutn difference between WFACT and tb'P calT 
eplated, prqpei 1 ant mass whleU is tpn?i dered signifi- 
cant enpugh t° Rqrfprm a p,ropeRart a?asf CQr,relat1on. 
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FIGURE 3-7 PROPEUANT MASS CORRELATION LOGIC FLOWCHART 
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AFt 

Tangent 

Plane 


Note: H Is replaced by H for Forehead Section Inputs 

FIGURE 3-11 STRAIGHT THROUGH GRAIN CnNFIGURATION INPUTS 
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star 





Slottcd-Cone 






circular 


Tabular Input } 

Port perimeter versus 


initial port area 


FIGURE 3-13 GRAIN CONFIGIIRATinn OPTIONS 
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FIGURE 3-20 

STAR GRAIN CONFIGURATION INPUTS 


- 3-90 . : V 




FIGURE 
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3-22 GIRCULAR PORT COMFIGURATIflN IfIPIITS 


3-101 


ihioat and Exit Diameters 



Conical Nozzle Half-Angle 


THETEX 



Approximate Bell Nozzle Half-Angl 


figure 3-23 nozzle CONFIGURATION INPUTS 
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■r 


AFT TANGENT PLANE — k— NOZZLE ENTRANCE PLANE 



FIGURE 3-24 SUBMERGED NOZZLE INPUT DATA 
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3.6 SAMPLE INPUT DATA FOR IBM 370 VERSION OF SRB-II 



CONTHOL UATA CARDS 
UNPun 

TITLE-*SAMP' , *LE S» ♦ »HM P» » 'ERFO» t *WHAN» » »CE P» * •W£DI > » *CTIO* • 'N 
NMl, 

NS1«4» , 

NLEXIS*0* 

HlTOT»9770,Of 

AENU 

lUMOZL DATA CARDS 
AONtONZ 
IUNPHT«0t 
AtND 

LEWIS Data 
REACTANTS 


C U 
H 1. 
0 1 . 
N 1. 
AL 1. 
CL It 


lit 15a S 298,15 F 
3.708 S 298.15 F 
39,078 S 298.15 F 
8.783 S 298.15 F 
15,999 S 298,15 F 
21.27A -72000, S 29a, IS F 


INSERT AL203IL1 

END i 

IBM data CARDS 
ilBUATA 
hL0«U58,9f 
DELZ»30. f 

AjNClNtUMO.0«llS.8ilia.8,Al4,6,Al7.6.698.1»701.1*996.9»999,9,1258,9* 
SA(ni«115,8»9l9.6,698.l iV96,9, 

SBa)«lia,a,4l7.6»70i.l,999.9, 

h(= m»lO*71.b, 

ANOtl)«2*U,,b*A,, 

TAUwa)«Z*8,70,36.A.33.S,36,A.,33,5,36,A,33,5,36.A,32,5, 

ALSl (1}>2*39,8, 

ADI (1)«2*17.0, 

AO2U)«2*17.0f 

AD3{U-2*17.0i 

R2{U»2*0.S, 

H5(n»2«2,0f 
DTAOd ) »2*l.,8*2,, 

DTAUWd ) «2*,b.B*2. » 

AN2*17,0, 

Dt»lA0,22, 

DT»b3.0, 

HUNBAk»578,7i 

EKBAR«0.0UAb, 

EREAP»1,0. 

SlhLAO«l., 


TbT=l., 

Dt LIST- 0.25, 
UtLTSS»2,0 , 


DELTTO-l.O, 
ANI rw«b. . 


T I MAX* 135,, 

PA»1A,7» 

PHI *29, A, 

ANR{2)»0,06b, 0.233, 
aKH i 36 J»0. 066 *0,233. 
AKSLOT m ■0,066,0.233 



ORIGINAL PAGE IS 
01 POOR QTJALIW 
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3.6 


(Continued) 


AKU<I)»5000.. 

SLTBRN (I ) «2 . 0 f 0 . 0 f 2 . 0 » 0 . 0 t 
NAK£N0«3* 

NAKH=9, • , 

nMAKR(l>»0.0»7.1»54, .71. 6189 , ilO&.JfHa.aflHA, 6*130, » 

■fBLAKRil)»2».U0105, ,00 132*. 0015.. 00 132*-. 000 ia«-,00l9S*-. 0048 *-,0048* 

AlNCPLa627,6* 

NTAUT0=1* 

NPH=S» 

TIHE.PH(1)»0, 0.0, 25*0. 5. 0,75*1,0* - 

PHST (1) *14, 7. 490, *725, *800,. 815.. 

U£LF=0.064* 

R=55,2* . 

NCSCOE»l* 

NCSTfla6* 

PRESSn)»0,*400,.600, *80 0. . I000,.la00,* 

CSTH U ) *5088,8 *5 125. 0*51 38,1 .5147,9*5154,4.5146,8. 

TCOMB M >*5950 , *6125, .6197. *6248, *6282, *6282. * 
AHwe(l)=28,,2B.27.2a,39*2a,47*28*54*28.54* ' 

GAHAG(l) *1,1612. 1.1615 *1.1616* 1.1618*1. 1619* 1.1619* 

' ULRF=0,01 . . 

0RVRF*0,01* 

AK=O,S0» 

AKK*0,S* 

DH1«36.* 

UNl = 80 , .,. : 

BH=1.0* 

BN=r.O* , ^ ^ ’ 

AOHM=75,42* 

A0NM=S5,62. 

UTAOE*1,0* 

R1G=1B. * 

HHR=35,75* 

TlHtCMa)50,0»31. .32,*3B,0*46,.b4,»58,*76,*110,*116,*120,*124. *140,* _ 

TULGH ( 1 ) *2**965* ,964 * ,962 *,960 * ,959* ,960 * ,961 * ,960 * ,958* ,954* ,949* ,947* 
PRTFLG=1.* 

nwrtab=i* 

tSMOICGsl * . 

VCHINP=764380.* 

VFHU=594380.* 

KPLAN£=10* 

NGEt)HD=12. , W VT n 

TAUMU(1)=0. *0.5*2. 0*4. 0*6, 0>8,7»9,0*10,0*12,0*14,0»16, 253*17.0* 

AbHL) (1 )=55000. *55300, *55400, *54700 , *52700 ,*48200, *44000, *34000 , *21000 , 1 

i0500.,2*0.0* : 

iENU 

CASE END / : : 

all UATA cards for this case have BEEN TRANSFEHED JO piSR 


3-1 Gi5 
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3.7 UNIVAC 1108 CONTROL CARDS 

CARD! @RUN{jSRB , 1 HNTSVXXXXXX . COLSONB IN226 ,2,200 

a b c 

CARD 2 (3ASG,T^SRB01 ,T,29312 

CARD 3 PFREE TPF$. 

CARD 4 0ASGJFPF$,F/1/P0S/3O 

CARD 5 eCfiPIN^SRBOl . JPF$. . 

CARD 6 PFREEJjSRBOl 

CARD 7 PASGJHJ, 23810 

CARDS @ASG,T}d8,F/100/TRK/150 

CARD 9 PASG,T^g,F/T00/TRK/150 

CARD 10 @ASG,Tf«10,F/inn/TRK/T00 

CARD 11 @ASG, 1^520, F/lOO/TRK/1 00 

CARD 12 PASG,T^25,F/100/TRK/100 ^ 

CARD 13 PXQTJiABS 

DATA CARDS 

CARD 14 @PriD,E 

CARD 15 @FIN 

The following is a brief description of each control card. Note that 
each card starts in column 1. 

CARD 1 

This is the job control card which contains the account number a * the engi neer 
or programmer’s name submitting the run b, and the BIN number c. The explanati 
of the other items on this card may be found in the MSFC/Computation Labbra- 
tory Programmer Procedures Manual. 

CARD 2 

This card defines the fiTename (SRBOT) and tape number (29312) which 
contain the SRB-II program. 

CARD 3 

This card sets free the Temporary Program File { TP P) storage area . 
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3.7 (Continued) 

CARD 4 

This card sets up space for the TPF storage area 
CARD 5 


TPF the SRB-II program, 
locatable, and absolute modules are qopfed. 


The symbolic, re- 


CARD 6 


This card sets free the tape unit from which the SRB- II program was read 
CARD 7 

Thermo^datf^nHf f ^ and tape number (23810) for the Lewis 
is used? when the NLEWIS > o option 


Module 


CARD 8 ^ ; ■ . ■ , - 

This card sets up temporary space on disk to be used by the BATES 
CARD 9 ■■■■.". ■ ■ 

■ Same as Card 8. 

CARD 10 

^’o’- *•’" (P^>- 

CARD 11 

IBM module 

CARD 12 

0" diste into which the dtspereton data 


CARD 13 


program*"'^ Ida'iTng and execution of the SRB-II 


3-107 
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3.7 (Continued) 

CARD 14 

This card will cause a postmortem dump to be Generated aftpr thP ny-nnv.=n, 
no^ostortem'X'is"'JakeJ!^ occurred normalTy 

CARD 15 


This card is used to signal that the end-of-run has been reached 
must be physically the last card in the total run deck. 


It 






4.0 


PROGRAM OUTPUT 
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SECTION 4 


A description and identification of SRB-II output is presented in this 
section. Some of the output is standard for a given type of simulation, 
but much of the output capability is optional and at the discretion of 
the program user. A detailed description of the separate outputs is 
presented in the paragraphs that follov/. 

SRB-II can utilize three methods of output: 

(1) Computer printout 

(2) Magnetic tapes 

(3) Punched cards 

These types of output are dispussed in Sections 4.1 , 4.2, and 4.3, 
respectively. Program error messages are dispussed in Section 4.4. 

4.1 COMPUTER PRINTOUT 

4.1.1 Input Data and Initial Values Printout 

The input data and ipitial values ppintop^ are standard forms of output 
which are exercised for every SRB-II run. The input data printout dupli- 
cates on paper the exact card inpup v;hich makps up the data deck. This 

printout is done on a case-by-case basis; i.e. , phe data packages for a 

given case are printed before the other output for that case. "The input 
data printout for a given case is Tdentified with the following heading: 

"CARD IMAGE OF INPUT DATA FOR CASE X" 

The initial values printout lists the values for constants^ indicators, 
dispersion limits and inert mass schedule for the case to be executed. 
Initial values printout format for each case is identified by the 
following heading: 

"THESE ARE THE INITIAL VALUES FOR CASE X" 

4.1.2 Module Data Printout 

Much of the printout for SRB-II is grouped by program modules. Separate 
printouts are given for: (1 ) Bates Specific Impulse Module, (2) Specific 

Impulse Scaling Mqdule, (3) Contractor Data Specific Impulse Module, 

(4) Dispersion Module, (5) Thrust Scaling Module, and (6) Internal 
Ballistics Module. Data from the Inert; Mass Mpdul^ are presented in 
the Thrust Scaling Module prini;op1; or the Internal Ballistics Module 
printout. The Reconstruction l^lodulp h^s no separate printout. Data from 
the Dispersion Module are refl^Ct^d in the Thrust pealing Module printout 
or the Internal Ballistics Module printout as well as' the initial print 
generated for dispersions. Data from the Reconstruction Module and the 
Nozzle Submergence and Contour Effects Module are reflected in the Internal 
Ballistics Module printout. 
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4.1.2 (Continued) 

cases data from a module are transferred to the Output Module 
actual execution of the output procedure. The basic printout 

®a^h parameter, the appropriate 
computed value for both 
International System (SI) and English units. 

4. 1.2.1 BATES Specific Impulse Module Printout 

The BATES _Isp Module calculates the “end item" vacuum Iso and 

velocitv for a given solid propellant at a 

Data Foi^ce devised BATES method. 

Data from this module are printed through the Output Module with the 

format shown in Fijure 4-1. tMs is a standard printout and Ifobtainad 
each time a simulation is made through this module. 

4. 1.2. 2 Specific Impulse Scaling Module Printout 

The Isp Scaling Module calculates the "end item" vacuum Iso and characteristir 

Ste^minini°th^V^ propellant by employing methods for ' 

determining theoretical Isp, divergence losses, and motor inefficiencies 

Data from this module are printed through the Output Module with the 

IZTkJlT ^ standard p?intoS? and Is oSSained 

each time a simulation is made through this module. 

4. 1.2. 3 Contractor Data Specific Impulse Module Printout 

IIIh Module calculates the "end item" vacuum Isp 

exhaust velocity for a given solid propellant using 
test data. Data from this module are printed 
through the Output Module with the format shown in Figure 4-3. This 

through“th1s''mSdiile?''‘,^'’'' “ simuiatlon'is made 

4. 1.2. 4 Thrust Scalino Module Printout 

The Thrust Scaling Module provides motor performance versus time prediction 
by scaling nominal performance data curves to ad.iust for variations 
in burn time and propellant mass for a particular solid rocket motor. 

Data from this module are printed through the Output Module with the 
format shwn in Figure 4-4. Data are printed for multiple time steps, 
but only the first three are shown in the figure. This is a standard 

printout and is obtained each time a simulation is made through this 
module. 
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X ■) 


I 

CJ 


THESE ARE the VALUES FROM THE BATES SPECIFIC IMPl’iSE MODULE 
PARA HFTF R INIEBHAnUNAL. 


bates test ANALYTICAL SPECIFIC IMPULSE? SIBATlH/StSt 


BS71.A392 


_E>d?U.5K_ 

262.213B7 


lAtis'jESf'sPEClFlC IMPULSE? SIBT(H/SiS) . 236A.775I 2A1.1A000 

.BATES- SPZCIELlJHBULSE_aiVEBQEIli:t_LOSS»_SILDyaiH/5i£l_._BajL7'>I33 8. 

bates specific impulse heat loss? SlLOb(M/S?S) 57.355026 B.BAHSBTO 

53 ^ 73 BA 9 “ 


bates specific impulse 2 PHASE FLOWLOSS? SIL2PB (K/S?S) 

bates .NOXOS5.Ji£ECinilJHRliLSEj_5jeNLIH/rStSI_ .... 

STANDARD-ANALYTICAL SPECIFIC IMPULSE? SIAS(M/S?S) 


5.A32A226 

,2S2D?S9BA-.. 257rO?95B„ 

2573.0195 262.37500 


END iTEtT ANALYTICAL SPECIFIC IMPULSE? SIAE(H/S?SI 26BS.76A9 273.B71B3 

_ FNH ITFh -specific IMPULSE divergence LOSS.SlL-OVEiM/S.S) _ 62,1S263A._ .. . . .6.337Bp5T_ 

end item .-specific impulse heat loss? SIL-0EIH/S?S1 IX.729193 I.I96DAA9 


end item SPEC. IMP. 2 PHASE FLOhLOSS? SI.L2PE IH/S iSI 

PELIVEREO EN D ITEM S PECIFIC IHPU LSE.t_SlflE 

divergence CORRECTION FDR NOZZLE? CLAHE 


38,766586 

..2980.6196. 

0,0 


3.9530926 
252.95289 


CHARACTERISTIC EXAUST VELOCITY? CSTAR IH/SiFT/S) 
rSTAH CU RVE FIT Cn EFFICTEHT JlfltUl. ' 


1979.8397 

9795.6172 


9855. lOlA 


CSTAR CURVE FIT COEFFICIENT NO. 2 


0.9S969539E-01 


CSTAR CURVE FIT COEFFICIENT NO, 3 


-0,35990190t-09 


FIGURE 4-1 BATES SPECIFIC IMPULSE MODULE SAMPLE PRINTOUT 
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FIGURE 4-2 SPECIFIC IMPULSE SCALING MODULE SAMPLE PRINTOUT 
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THESE ARE THE values FflOH THE CONTRACTOR data SPECIFIC IMPULSE MODULE 

' ^PTSThYTEIT ~~ ^ ~ I NTEHNAT loN^L " ^ ENCLTSH""* “ 

contractor analytical specific impuls e ^ SIACQNlH/StS l 236B.13B 7 2a 1 .48207 

■ END- ITEM anal ytic al 5PECXFrC~tHpDLSET 'SIAE tH75t5r " '~ ‘ 2605,7640 ’ ' 273.07183" ~~ 

DEnVEfiEd"'END^TWSPEClf rC~fH"fn4rs£iTsiDTt«/‘l3,‘'sr''”‘'2S69.20l'7 ■■ ~ 261". 99407 “ — 

EN D ITE M FLOWRATE. R O E tKG7S»L8H/S) SfA U a477_ , _12B79 . 07fl_ 

TEST MOTOR FLOWRATE. W02 IKG/S.lBh/S) ~ ' ; 4.5350230 lO.OQOflOO"" " V 

— END'TTEM-"NOZZLrEFFICIENCfT'ETANZ““^ — “ — 0,95066221 ' ' “ ’ : 

T EST M OTOR NU ZZLE efficiency. ETANZl 0.903B5562 " . 

divergence CORRECTION FOR' NOZZLE. -CL AHE " 0,97605051) _ 

"~DTVEHfiERCE'CDRRrcf rOB"p7IR'Til02ZCETT:nHl — 9029629 1 " 

E ND i te m CSTAR. CSTAR IM/S.FT/S) 143 6.4460 0712.7 500 . - ■: 

: ; CSTAR CURVE Fit C OEFTICIENT~7lo ; ~l~~ ~'' " " 4674.3672 " ' ^ ' 


CSTAR~caHVE">'rr"rOEFIT:CIEm~NlT7"2 * OLSEAGBSSAE-Ur"^ 

CSTAR CURVE FI T coef fi cient NO. 3 i -0.3S990190E-04 



FIGURE 4-3 CONTRACTOR DATA SPECIFIC IMPULSE MODULE SAMPLE PRINTOUT 
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MOOULR 





SB“ ir.SITin.N TtvB. TP^tSl C.lOCCOOOOe 

01 





V‘B(,(C1= pc^ivii^TPBS 


TNTFRNATIONAL 

ENGLISH 


SPV total Th'.(>;t, rS=w|-BUNiUPF) 
scv !^.^CT rn-llST. Pltf'.LRP) 


0.1B332a7op 

0.25*t'T3S2OE 

09 

OS 

0.2974Q1S0R 

3.57311719= 

07 

04 

SP't total Vict; F(rv;pAT = , tps R V+P IFtKP/S t LPH/SI 
gCM jAiCOT MAIS FLOucArr , mi s { kG /R , LB"/ 5 1 


0.50E2P056E 
0.13F 3213EE 

OA 

0? 

0.112C56SC= 

0.30l440£,aR 

05 

0? 

5f,>- psrPDLL.-AT SPECI=IC I«0JLRF, SHP («/S,S) 
St" RPCflFir I'OiiLS= np IKF6TS. SII(M/S,S1 


Q.2 405A-qAAF 
0.1 »n;2632E 

OA 

04 

0.2656E552E 

1.1°000030R 

02 

0*^ 

SPM c,F=F|cT(- IvpllsF, SR«RI ("/R , El 
S = « PFOPri 1 AVT MASS RPVAININn, Mpp {kg, LBMI 


0.25'3’-954F 

0.75765213F 

04 

04 

0 .36548267 = 
0.15TC7 7R JF 

0" 

07 

SSM IA1--I.T MARR RFMA.rT.r'r., MJU (kg ,LR(a} 

S=A‘ rHAMHik OFESStlRF, PC IN/-ia»2, LPP/IN«-«'l 


0.6R2FPA14F 

0.35170334C 

0 4 
09 

Q.lSij 

3-510in471P 

05 
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0.512352R ir 

05 

02 
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4. 1.2. 5 Dispersron Module Printout 

The Dispersion Module is used to perturb the Internal Ballistics 
Module input data about the nominal values. These perturbed data 
are used to generate a dispersion prediction from the IBM. Output 
for a dispersion prediction consists of the value of the indicator 
NDIp, a dispersion options catalog which lists all available dispersion 
options, the final value of the dispersed variable, and normal IBM 
prints. (Paragraph 4.1 .2.6) . The Dispersion Module print is generated 
in the Output Module with the format shown in Figure 4-5. Formats 
shown in Figure 4-5 are standard for dispersion predictions and appear 
in the output of each dispersion case immediately follovnng the initial 
values print. 

* 

4. 1.2. 6 Internal Ballistics Module Printout 

The InternaT Ballistics Module performs an in-depth analytical simulation 
of solid rocket motor performance as a function of time. This module 
has one standard printout and severaroptfonaT printouts. Each printout 
capability is discussed in the following paragraphs. 

Motor Performance Data Printout ; For predictions in which an internal 
Dal I i sties simulation is performed, a printout of motor performance 
data will be made following each computation time step. A sample 
printout for a large solid rocket motor performance prediction for a 
motor with no nozzle submergence is presented in Figure 4-6. The 
standard printout for a motor with nozzle submergence is identical to the 
one shown in Figure 4-6 with the exception of substitution of data for 
regions A and B for aft-end parameters. A complete printout of motor 
performance data for an internal ballistics simulation, incl udes this data 
for every program computation time step (see Volume III of this document). 

The computation times are controlled by program inputs. This printout 
is executed by the Output Module. 

Increment Pi vi ding PI ane Data Printout : An optional expanded printout 

controlled by PRTFLG (see Table 3-XXII for definition) is available to 
print additional data for the end sections and mass addition regions 
between each increment dividing plane (I.<D.P.). A mass addition region 
is defined as the region between two increment dividing planes (e.g., 
mass addition region 1 would lie betvyeen increment dividing planes 1 
and 2). The values printed are the parameters at the dov/nstream plane 
of the mass addition region except for the generated weight flow, DWD0T, 
which is within the region. 

T.D.P. data are printed with the format ^hown in Figure 4-7 for motors 
without nozzle submergence. When the I,D.P. printput option is exercised, 
the data are written below the motor performance data shown in Figure 4-6 
and corresponds to the same program computation time as the preceding 
motor performance data. This printing is controTTed by subroutine I8M0UT 
of the Internal Ballistics Module. The following list defines the parameters 
given for each T.D.P. : - y ; 
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IfcNITIDN TIHE« TIHE IS) I.OOOODOO 

parameter description tUNITSl 

PPOGRAH 

HOHENCLATURE 

INTERNATIONAL 

ENGLISH 

motor PARAMETERS) 

TOTAL DELIVERED THRUST IN,LflFl 

ftoel 

11329A99, 

2566973.0 

TOTAL VACUUM Thrust iHtLBFt 

FTVAC 

1H3392A9. 

2773976.0 

THRUST CONTRItjuTION OF INERTS IN.LBFI 

FI 

161003.31 

36195.000 

DELIVERED total IMPULSE lM»S,LaF‘Sl 

SRMOTI 

6156I3B.0 

1833573,0 

VACUUM TOTAL IMPULSE (N*S.LaF«S) 

SRHVTI 

9039666.0 

2032198.0 

Thrust coefficient 

CF 

1.6622987 

1.6622987 

DRAIN DISCHARGE HASS FLOMHATE IKD/S.LBH/SI 

WOOT 

M791.57B1 

10563.629 

FLORKATE integral (KG*LUH) 

SWOOTK 

JS23.6AS1 

7767.8672 

Inert mass Flowrate iitG/SiLBH/Si 

mif 

B6.A09332 

190^50000 

INERT HASS remaining IKGiLSH) 

HIR 

M37I.0313 

9636,6766 

TOTAL BURN AREA !M»<atlN»»21 

ABTOT 

365.62231 

535715.96 

TOTAL PROPELLANT VOLUME (H**3.1n»«3( 

VF 

fiS4. 92639 

1S5S666B. 

PROPELLANT HASS REMAINING IKG.LBH) 

wr 

‘*51601.96 

995611,88 

TOTAL GAS HA!iS tf^GiLUK) 

WGTOT 

533,19507 

1175.6966 

RATIO OF specific HEATS 

GAMA 

1.161B061 

1.1618061 

molecular veight of gas (KQ/M0LE»LBH/H0LE> 

AHW 

12.916168 

28.675235 

characteristic exhaust velocity 1K/SiFT/S) 

CSTAR 

1568.5613 

5166,1992 

HaAIMUH CHAHUSS HACH 

AHPK 

0.29865697 

0.29865697 

head end parameters: 
total pressure (NFH«»2,LBF/IN*»2) 

PH 

0.56192256E 09 

SIS. 00000 

PRESSURE INTEGRAL. IN*S/H»»2 ,LBF •S/IN»*2) 

. SPHDT 

0.61B83036E 09 

607.66216 

Burn AREA (M»»a,iN»*2i 

AHH 

35.570129 

55133.816 

BURN rate (H/S, in/si 

RBZIll 

0.T9S65D83E-02 

0,31662958 

DISTANCE BURNED |M,|N| 

TAUzm 

D.56669856E-02 

0.22303106 

PROPELLANT volume iH*»3tIN«*3) 

VFh 

9.5388773 

582098.31 

GAS VOLUME 

VPH 

2.9870596 

182281,69 

GAS STATIC TEMPERATURE lOEO KiOEG R) 

PRNT|lt3) 

3672.3760 

6250.2773 

CYLINDRICAL SECTION PARAMETERS) 
RADIAL burn AHEA IH—2,1N*»21 

A0CYL 

255.87112 

396601.06 

SEOHENT FACE BURN AREA IM**2iIN«»2J 

ABSLOT 

66,632126 

69179,938 

propellant VOLUME 

vfctl 

238.05310 

16526898. 

GAS VOLUME tH»«3,IN»»3) 

VP 

. 89.761520 

5677587.0, 

aft end PARAHETEfiSJ 

TOTAL PRESSUHEINOZ ENT) t«/H*»2,LBF/'IN»«2) 

PON 

0.S2772S57E 09 

765.60137 

PRESSURE integral (N*S/M*»2iLBF«S/1N»«2) 

SPONDT 

0.3B901350E 09 

566.21655 

BUHN AREA . 

AAN 

9.5691009 

16801.137 

BURN RATE (H/S, IN/SEC) 

RBZINI*!) 

0.779S82S6E-02 

0.30692238 

uistance burned (H»1N) 

TAUZINI*!) 

0.SS107139E-Q2 

0.21695727 

PROPELLANT volume tH»*3.TM*»3) 

VFN 

7,3326032 

667650,96 

OA5 VOLUME 

VPH 

6,.b93H73S 

286601,25 

GAS STATIC TEMPERATURE tOEG K.OEG R) 

PRNTlNIt3) 

3651.3865 

6212,6961 

PORT AREA IH»*2aN»»21 

AP 

3.1171980 

6831.6680 

NOZZLE PARAMETERS) 

Throat area iM»«2,iN»»2) 

AT 

1.6237328 

2206.7910 

EXPANSION ratio 

EPR- 

6.9975882 

6*9975882 

PRESSURE RATIO 

repo 

0.2216306SE>01 

0.22163065E-01 

miscellaneous PARAMETERS) 

ANISOTROPIC BURN RATE COEFFICIENT 

AKHST 

0.659 5160 7E-01 

0.65951607E-01 

number of PRESSURE ITERATIONS 

ANLOPS 

3.0000000 

3.0000000 


FIGURE 4-6 INTERNAL BALLISTICS MODULE SAMPLE PRINTOUT OF MOTOR 
PERFORMANCE DATA ■ 
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INCREMENT 

dividing 

PLANE DATA! 












HASS 

PO 

P . 

T 

u 

H 

UP 

AP 

WOOT 

DWDOT 

DW/DT 


TAU 

RBTO 

7 AUTO 

addition 

REGIONS 

P5IA 

PSlA 

DEG. R 

FT/SEC 


INCHES 

so. IN. 

lb/sec 

LB/SEC 

LB/SEC 

IN/SEC 

IN. 

IN/SEC 

IN. 

fOHE 

815.00 

814.75 

6250.28 

82.34 

0.023 

1165.69 

5578.44 

1103.23 

1109,46 

6,25 

0.3144 

0,223 

0.3132 

0.223 

DfliDDO 

814.78 

814,10 

6249.82 

134.65 

0.038 

1165,69 

5578,35 

1802,23 

703.71 

4.70 

0.3144 

0.223 

0.3144 

0.223 

6{]>000 

814.35 

813,04 

6249.15 

187.09 

0.053 

1165.69 

55/8,26 

2501.13 

703.62 

4,71 

0.3143 

0.223 

0.3143 

0.221 

90i000 

813.71 

811.57 

6248.25 

239.76 

0.067 

1165.69 

5578,17 

3199,86 

703.45 

4.72 

0.3143 

0.223 

0.3143 

0.223 

115.576 

812.98 

809,96 

6247.31 

284.93 

0.080 

1165.69 

5578.11 

3795.64 

599.49 

3.70 

0.314] 

0.223 

0.3141 

0.223 

Ua.BDO 

812.61 

805.77 

6243.1B 

429.74 

0.121 

221,94 

3919.72 

4004.65 

472,96 

1.78 

0.314D 

0.223 

0.3140 

0.223 

l^B.BOO 

811.92 


6242,93 

A37*09 

0.123 

223.79 

3985.3] 

4136*69 

134.27 

2.23 

0.3136 

0.222 

0.3136 

0.22? 

178»BOO 

811.21 

803.90 

6242.68 


0.125 

225.63 

4051.42 

4269.70 

135.2B 

2.26 

0.3)35 

0.222 

0.3135 

0,22? 

2QB.BOO 

810.49 

802.96 

6242.1.3 

451,33 

0.127 

227.48 

4118.06 

4403,75 

136.36 

2.30 

0.3134 

0.222 

0.3134 

0.22? 

23a<600 

809.76 

802.00 

6242,17 

458.23 

0.129 

229.33 

4185.23 

4538,84 

137.43 

2.34 

0.3133 

0,222 

0.3133 

0.22? 

368. BOO 

B09tD2 

801. 

6241.93 

^6S« QD 

0.131 

231.17 

42=2,93 

4674.95 

138.50 

2.38 

0.3133 

0.222 

0.3133 

0.22? 

29B.800 

808.28 

800.08 

624], 6B 

471.64 

0.133 

233.02 

4321.17 

4B12.I1 

139,58 

■ 2.42 

0.3132 

0.222 

0.3132 

0.22? 

328,800 

807,52 

799,10 

6241.43 

478.14 

0.134 

234,87 

4389,93 

4950.29 

140.65 

2.46 

0.3)31 

0.222 

0.3131 

0.22? 

35B.BD0 

806.76 

' 798.12 

6241.19 

484.53 

0.136 

236.72 

4459.21 

5089. 5D 

141.72 

2.50 

0.3130 

0.222 

0,3130 

0.22? 

3BB.800 

RD5.99 

797,14 

6240,94 

490.79 

0.138 

238,56 

4529.04 

5229.75 

142,79 

2.54 

0.3)29 

0.222 

0,3129 

n,22? 

Bl^i.376 

805.33 

796.30 

6240.73 

^96.{)6 

0.139 

240.14 

4588.98 

5350.41 

122.58 

1.91 

0.3)28 

0.222 

0.3128 

0,22? 

417,823 

803.40 

780,58 

6234,35 

637,28 

0.179 

221.95 

3919.99 

5820,53 

228.40 

2.82 

0.3127 

0.221 

0,3127 

0.221 

447.600 

R02.IB 

787.19 

6234.14 

641.48 

O.IBO 

223. BP 

3986.64 

5951.21 

132.70 

2.02 

0.3)20 

0.221 

0.3120 

0.22) 

477.600 

ROD. 96 

785.79 

6233.93 

fiO 

0.182 

225,83 

4058.39 

6083.56 

134.68 

2.32 

0,3)19 

0.221 

D.3119 

0.221 

507.600 

799.73 

784.40 

6233.72 

649.65 

0.183 

227.77 

4128.73 

6216.99 

135,79 

2.36 

0.3)18 

0.221 

D.311S 

0.22) 

537.600 

79B.50 

7R3.01 

6233.51 

653.63 

0.184 

229,72 

4199,66 

6351,49 

136.90 

2.40 

0.3)16 

0.221 

0.3116 

0.221 

567,600 

797.27 

781.62 

6?33«31 

657.54 

0.185 

231.67 

4271. IB 

6487.05 

138.01 

2.44 

D.3I15 

0*^B[) 

0.3115 

0.220 

597,600 

796*05 

7fl0.24 

6233.11 

661.39 

0.186 

233.62 

4343.29 

6623.68 

139.12 

2.48 

0.3)14 

0.220 

0.3114 

0.220 

627.600 

794.83 

778.85 

6232.91 


o.iav 

235,57 

4415.99 

6761.37 

140.22 

2.53 

0.3113 

0*220 

0.3113 

0*220 

657,600 

79J*bl 

777.49 

6232,71 

668*91 

o.iaa 

237,51 

4489.27 

6900.13 

141,33 

2.57 

0.3111 

0*220 

0.3111 

0.220 

687.600 

792.39 

776.12 

6232.51 

672.58 

0.1S9 

239, A6 

4583.16 

7039.94 

142,43 

2.61 

0.3110 

0.220 

0.3110 

0.220 

697.878 

791*97 

775.65 

6232.44 

673.85 

0.189 

243.13 

4586.60 

7088.35 

49.05 

0.63 

0.3109 

0.220 

0.3109 

0.220 

701,100 

790.76 

766.57 

6223.56 

822.64 

0.232 

221.92 

3919,15 

7314.94 

467.36 

1.76 

0.3108 

0.220 

D,3lOB 

0.220 

731.100 

789.82 

765.54 

6223.43 

824.62 

0,232 

223.77 

3984,72 

7445.43 

132.81 

2.31 

0.3)00 

0,219 

0.3100 

0.219 

761.100 

788.88 

764.53 

6223.31 

826,56 

0.233 

225.62 

4050.82 

7576.79 

133.70 

2.35 

0.3099 

0,219 

0.3099 

0.219 

791.100 

787.95 

763.52 

6223. IB 

828,47 

0.233 

227,46 

4117.45 

7709.16 

134,76 

2.38 

0.3098 

0.219 

0.3098 

0.219 

821.100 

787.02 

762.51 

6223.06 

830.36 

0.234 

229.31 

4184.61 

7842.55 

135.82 

2.42 

0.3o9T 

0.219 

D.309T 

0.219 

851.100 

786.11 

761. 5i: 

6222.93 

832.21 

0.234 

231.16 

4252.30 

7975.97 

136. BB 

2.45 

0.3096 

0.219 

0.3096 

0.219 

881.100 

785.20 

760.52 

6222.81 

834.05 

0.235 

233.00 

4320.52 

8112.41 

137,93 

2.49 

0.3Q95 

0.219 

0.3095 

0.219 

911.100 

78<4.29 

759.54 

6222.69 

835.86 

0.235 

234,85 

4389.27 

6248.87 

13B.99 

2.53 

0.3094 

0.219 

0,3094 

0.219 

941.100 

783.39 

758.56 

6222.57 

837.65 

0.236 

236.70 


8386.34 

140.04 

2.56 

0.3o93 

0.219 

0.3093 

0,219 

971.100 

782.49 

757.59 

6222.45 

839.42 

0.236 

238.55 

4528,36 

8524.84 

141,10 

2.60 

0.3092 

0.219 

0.3092 

0.219 

996.679 

781.09 

756,10 

6222.30 

841,66 

0.237 

240,12 

4588.30 

8643«99 

121.14 

1.98 

0.3091 

0,219 

0.3091 

0.219 

1000.120 

777.86 

738.71 

6205.75 

1059,83 

0.299 

221.93 

3919.57 

9108,56 

0.0 

2.72 

0.3090 

0.219 

0.309C 

0.219 

1029.900 

776.50 

738.21 

6206.68 

1048.92 

0.296 

224,75 

4019,80 

9237.64 

131,17 

2.08 

0,3073 

0.217 

0.3073 

0.217 

1059,900 

775.18 

737.72 

6207.57 

1038.26 

0.293 

227.59 

4122.02 

9358.69 

133.44 

2,38 

0,3073 

0,217 

0,3073 

0,217 

1089.900 

773.91 

737.22 

6208.41 

1028.01 

0.290 

230,^2 

62if5«4B 

9501.35 

135,09 

2.43 

0.3072 

0.217 

0,3072 

0.217 

1119.900 

772.66 

736,72 

6209.22 

1018,15 

0.287 

233,26 

4330.20 

9635.60 

136,74 

2.49 

0,3072 

0.217 

0,3072 

0.217 

) 149,900 

771,45 

736.22 

6209.98 

1008.65 

0.284 

236.10 

4436,16 

9771.44 

138.39 

2,55 

0.3071 

0.217 

0.3071 

0.217 

1179.900 

770.27 

735, 72 

6210.72 

999.50 

0.282 

236.94 

4543.38 

9908.88 

140.04 

2.60 

0.3071 

0,217 

0.3071 

0,217 

I209»900 

769.12 

735.21 

6211.42 

990.68 

0.279 

241,77 

4651.84 

10047,91 

141.70 

2.66 

0.3q7o 

0,217 

0.3070 

0.217 

1239,900 

768,00 

734.71 

6212,09 

982.18 

0,277 

244.61 

*54 

1C1BB.54 

143.34 

2.72 

0.3o7o 

0,217 

0.3070 

0.217 

1258,900 

757,31 

734.39 

6212.50 

976.95 

0.275 

246.41 

4831,67 

10278,41 

91.64 

1.75 

0.3069 

0.217 

0*3069 

0.217 

AFT 

765.40 

730,42 

6209.95 

1009.1] 

0.284 


4831,67 

10563,63 

290.74 

5.52 

0.3069 

0.217 

• • « * 
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4. 1.2. 6 (Continued) 


P0 

Total pressure 

p 

Static pressure 

T 

Static temperature 

U 

Gas velocity 

M 

Gas mach number 

LP 

Port perimeter 

AP 

Port area 

WD0T 

Discharge mass flow rate 

DWD0T 

Generated mass flow rate 

DW/DT 

Mass buildup in region 

RB 

Isotropic burn rate 

TAU 

Isotropic burn distance 

RBT0 

Anisotropic burn rate 

TAUT0 

Anisotropic burn distance 


If the motor being simulated contains slots, data for the variables 
listed above will also be included for each slot in the I.D.P. printout. 
This data are printed for the slot forward and aft interface locations. 

Center of Gravity and Moment of Inertia Data Printout : When the option 

to calculate propellant center of gravity (C.G.) and moment of inertia 
(M.0.1.) data is exercised, a printout of the data is given below the 
motor performance data for each time step. The C.G. and M.O.I. cal- 
culations are controlled by KM0ICG and the printout is controlled by PRTFLG 
(See Table 3-XXII for definitions). This output is generated by subroutine 
IBM0UT of the Internal Ballistics Module. The following variables are 
printed: 


HEADING 

(FORTRAN NAME) DESCRIPTION 


GRAIN CG Motor center of gravity location from aft tangent 

(XBARIH) plane, (M, IN). 

GRAIN ROLL MO I Motor roll moment of inertia about Tongitudinal 

(AJPP) axis, (KG-M**2,SLUG-IN*’^2)>, 


4-11 


0256-10020-2 


4.1 .2,6 (Continued) 


HEADING 
(FORTRAN NAMEl 

GRAIN XY HOI 
(AJBB) 

head END CG 
(XBH) 

head END ROLL HOI 
(AJPHED) 

HEAD END XY MOI 
(AJBHED) 

CYLINDER CG. 
(CGCYL) 

CYLINDER ROLL HOI 
(AIPCYL) 

CYLINDER XY MOI 
(AIBCYL) 


DESCRIPTION 


Motor pitch monient of inertia about grain center 
of gravity, (KG-M**2, SLUG- IN**2). 


Head end section center of gravity from forv/ard 
tangent plane, (M, IN). 

Head end section roll moment of inertia about 
longitudinal axis, (KG-M**2, SLUG-IN*‘^2) . 

Head end section pitch moment of inertia about 
forward tangent plane, (KG-H**2, SLUG-IN**2). 

Cylindrical section center of gravity location 
from aft tangent plane, (M, IN). 

Cylindrical section roll moment of inertia about 
longitudinal axis, (KG-M**2, SLUG-IN**2). 

Cylindrical section pitch moment of inertia about 

aft tangent plane, (KG-M**2, SLUG- IN**2). 


AFT END CG 
(XBN) 


Aft end section center of gravity from aft 
tangent plane, (M, IN). 


AFT END ROLL MOI 
(AJPNOZ) 


Aft end section roll moment of inertia about 
longitudinal axis, (KG-M**2, SLUG-IN**2). 


AFT END XY MOI 
(AJBNOZ) 


Aft end section pitch moment of inertia about 
• aft tangent plane. (KG-M**2, SLUG-IN**2). 


^Qni|try Tables Printout: An optional printout is available to write 
Hterna] Bal listics Module geometry tables data. Following the 

geometry tables for the reference plane perimeter and 
radius of gyration: versus burn distance, and the end sections burn area 

moment of inertia versus burn distance are completed 
each tab! e^is printed if the input NWRTAB is set to 1 . This prinSut 
IS obtained prior to the printing of the motor performance tiSe-step^data. 

Di agnostic Data Dumps : Several intermediate, diagnostic print statements 
have been permanently built into the Internal Ballistics Module to^a^^^ 

analysis. These user-controlled dumps are 
oJtnu?^ inputs CKDUMP( I) and KDUMP{I). The KDUMP 

L calculation interval unless they are con- 

trolled by the CKDUMP inputs. A description of these inputs follows : 
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4. 1.2. 6 (Continued) 



INPUT 

DEFINITION 

SUBROUTINE 


CKDUHP(T) 

Time at which KDUMP(71) will start 
printing KDUMP 1, 2, and 3 (SEC). 


CKDUMP(2) 

Time at which KDUMP (71) will stop 
printing KDUMP 1, 2, and 3 (SEC). 


CKDUMP(3) 

Time at which KDUMP (5) will start 
printing. (SEC) 


CKDUMP{4) 

Time at which KDUMP(5) will stop 
printing. (SEC) 

— 

KDUMP(I) 

Writes KDP, III, lid, TAUZ(Xn), 
RBZ(IIO), P, SUMDV, VF, ASE, AP, 
WD0T, DWD0T. 

ASESUB, ASTSUB, & 
HNDSUB 

KDUMP(2) 

Writes PIN, PD(III) , PDPR(III-1 ), 
PDPR(III), V(III), VPR(III), APHI, 
AP, U. UTMP, WD0TI, ‘DWD0T, WD0TD, 
R, T, DELT, III, PDY, TEMP,ACCEL. 

AIBST 


Writes IFLAG, PD{III), PDY, WD0TD, 
PDX2, PDY2, El, E2. 

AIBST 


Wr i tes ITER , TXT , TYl , TX2 , TY2 , 
TEl , TE2. 

AIBST 


Writes P, APHI, AP, DELTA, U, R, 
T, WDf3TI, WDPT, DWD0T,' HI. 

AIBSUB 


Writes SLTFLG, IIS, ISl , IS2, AINCX, 
AINCY, AINCW, SCUR(IIS,1), SCUR 
( IIS , 2 ) , TSL|3TF( I IS ) , TSL0TA( IIS ) , 
III, ZCALC(III). 

SEGSUB , 


Writes PH0LD, WD0T, WDB, DELLP, 
PMIN, WD, CL^PS, P0N, DIS, DISB, 
DEED, PMAX, PH, ANLPPS. 

SEVPH 

KDUMP(3) 

Writes IC0UNT, PSA(IIS), PDY, 
WD0TD, PDX2, PDY2, El, E2. 

SL0T 


Writes ITER, TS, T, TX2, TY2, 

TEl, TE2. 
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4.1. 2. 6 (Continued) 
INPUT 


KDUMP(3) 


KDUMP(4) 


KDUMP(5) 


KDUMP(6) 


KDUMP(7) 




definition 


apf(iis), uf 

P0F(IIS), PSF(IIS), TSF 
PD(ISl), WSL0TD(IIS) 
APA(IIS), UA(IIS), P0A(IIS), 
PSAdlS), TSA(IIS), PD(IS2), 
JWfLTFdlS), RBSLTFCIIS), ABSLTF 

ABSLTAdlS), DWDTS(IIS). 


SUBROUTINE 

SL0T 


Writes AFP, ASLVR, BVIM, BVl , 

BV2M, BV2, BV2P, ETA22, ETA2, 
ANGLE, ETAll, ETAl , BVX, BVXX, 

THV, THSLVV, RSLVRV, TAU, RF, III. 


LPT0 


(AL{I),I=1,4), AL(7), AL(8), 
AL7, AL8, OWD0T, DW7, DW8, RB7, RE8, 
AKRT0, TAUT0W. AKRT0X, TAUT0X, 
AKRT0M, AKRT0^/ n, NN, M, III. 


SEGSUB 


SfT NTME, SLTFLG, flMCKl! TWK2. 
T0FLAG (This dump will also be 
written whenever an error situation 
occurs such that ICHN is set to 5}' 


Wn'tes AKRST, DIFAKR, WD0T, DIS, WD, 
DEED. 


Writes PCTAB, DIFAB, WD£)T, DIS, WD, 
DEED* 


C0NV 

C0NV 


Writes AP, ALP, AFF, ALPX, ALPY, APX, SEGSUB 
APy, AINCW,AINCX,AINCY, TEMP, III. 


This flag overrides KDUMP(l), 
KDUMP(2), and KDUMP(3) settings and 
prints these dumps according to the 
simulated burning times specified 
in CKDUMP(l) and CKDUMP(2) 


MNCHN4 


turnout Data Printout: When web burnout occurs in anv increment riivi Hi nn 

WEB BURNOUT INCREMENT LOCATION 
XXXX.XXX 
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4. 1.2. 6 (Continued) 

In segmented motors, each time a slot burns laterally through an 
increment dividing plane, the statement below is automatically printed to 
identify the location of the dividing plane: 

LATERAL BURNOUT INCREMENT LOCATION ZCALC(I) = XXXX.XXX 

4,1.3 Miscellaneous Optional SRB-II Printout 

LEWIS Subroutine Printout ; The LEWIS Subroutine is a modification of 
the CEC71 program developed by NASA Lewis (Reference 3). LEWIS is used 
in all three Isp modules and the Reconstruction Module to calculate 
theoretical rocket performance parameters such as Isp and characteristic 
velocity. Printed output from LEWIS is optional in SRB-II and is obtained 
by setting the indicator IRKTOl = 1 (see Paragraph 3.4.7). A sample 
printout for the LEWIS Subroutine is shown in Figure 4-8. 

IDNgZL Subroutine Printout : The IDN0ZL Subroutine is a modified version 
ot the Une-DImensional Two-Phase Flow Loss Computer Program developed 
by Aerospace Corp. (Reference 4). IDN0ZL is used in the BATES Module 
to calculate the two-phase flow loss for the "epf^ item" rocket motor 
and/or the BATES test motor. Output from IDN0ZL is optional in SRB-II 
and is obtained by setting the indicator IDNPRT = 1 (see Paragraph 3.4.7). 
Extended print for the IDN0ZL subroutine is obtained by setting the 
indicator IDBUG =1. 

The first values that appear in the IDN0ZL printout are the NAMELIST 
input data values, propellant constants, tolerance limits, the results 
of all initial calculations, and the following calculated nozzle 
parameters are listed: 

XT X coordinate at throat (FT) 

X| X coordinate of transition point between linear segment 

and circular segment of the nozzle upstream of the throat 
(FT) 

R-] nozzle radius at X = X-) (FT) 

X 2 X coordinate of transition point betv/een linear segment and 

circular segment of the nozzle downstream of the throat (FT) 

R 2 nozzle radius at X = X 2 (FT) 

XF X coordinate at exit (FT) 

RF nozzle radius at exit (FT) 
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FIGURE 4-8 LEWIS SUBROUTINE SAMPLE PRINTOUT 
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4.1 .3 (Continued) 

The following three parts of tabular data are also printed: 
PART 1 OF TABULATED IDNOZL PRINTOUT: ’ 


The follov/ing quantities associated with the "LAR" case are listed as 
functions of X: 

S Distance dov/n the nozzle from the end of the chamber (FT) 

M Mach number of the gas 

UG Gas velocity (FT/SEC) 

UP Particle velocity (FT/SEC) 

TG Gas temperature (°R) 

TP Particle temperature (°R) 

KD Drag correction factor 

KH Heat transfer correction factor 

.REP Reynolds number of the particles 

ETAG DUg/dx (FT/SEC/FT) 

H Enthalpy of the particles (FT**2/SEC**2) 

FX Surface tension (DYNES/ CM) 

P Pressure at X (PSIA) 

A Area at X (FT**2) 

R Radius at X (FT) 

PART 2 of TABULATED IDN0ZL PRINTOUT: 

The following quantities associated with the "LAG" case are listed as 
functions of X: 

X Distance dovm the nozzle from the end of the chamber (FT) 

M Mach number of the gas 

A/ATH area ratio 

UP/UG Ratio of velocities 
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4.1.3 (Continued) 

TG/TP Ratio of temperatures 

I(VAC) Vacuum specific impulse, no lag (SEC) 

I(0PT) Optimum specific impulse, no lag (SEC) 

WR0STC Ratio of gas density to stagnation gas, no lag 
P/PSTG Ratio of pressure to stagnation pressure, no lag 

I/IREF(V) Ratio of "LAG" to "NO LAG" vaouun specific Impulses 
I/IREFW) Ratio of "LAG" to "NO LAG" optimum specific Impulses 

part 1 Of the tabulated printout win alL be"wHtten™*’°" 


4.2 


MAGNETIC TAPE OUTPUT 


execution of the output is controlled by the Output Module. 

Internal Balli stics Data Tape fI/0 llrn't ip^. Thnc -t--, •t-y t_ 

following variables will be written after each computation time step: 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


( 8 ) 

(9) 


Time from ianition (SEC) 

Vacuum thrust (LBF) 

Vacuum total impulse (LBF-5EC) 

Vacuum specific impulse (SEC) 
Characteristic exhaust velocity (SEC) 
Thrust coefficient 
Head-end total pressure (PSlA) 

Nozzle entrance total pressure (PSIA) 
Head-end static temperature (°R) 
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4.2 (Continued) 

(10) Nozzle throat area (IN**2) 

(Tl) Maxi mum combustion chamber mach number 

(12) Grain discharge mass flov;rate (LBM/SEC) 

(13) Propellant mass consumed (IBM) 

(14) Propellant mass remaining (LBH) • 

(15) Inert mass remaining (LBM) 

(16) Nozzle exit area (IN**2) 

(17) Nozzle expansion ratio 

(18) Propellant grain center of gravity referenced to the aft 
tangent plane (IN) 

(19) . Pitch/yaw moment of inertia (SLUG-IN**2) 

(20) Roll moment of inertia (SLUG-IN**2) 

Internal Ballistics Plot Tape (I/O Unit 13); This tape vnlT he written 
when the indicator NPL0T is set to 1 in the Control Data Package and an 
Internal Ballistics Module simulation is executed. The generated tape 
represents an intermediate step in obtaining hard copy plots. The tape 
can be reformatted for use by existing Benson and Lerner or MSFC SC 4020 
Plot Programs. 

The variables written on the plot tape are the same as those listed above 
for the performance data tape and are written in the same order. 

4.3 PUNCHED CARD OUTPUT 

SRB-II has two optional punched card output capabilities. These are the 
booster/mass property data cards and the Internal Ballistics Module 
geometry table cards. 

Booster Mass Property Cards : Booster mass property cards are punched 

when the indicator NCARD is set to 1 in the Control Data Package 
and an InternaT Ballistics Module simulation is being executed. Cards 
are punched for each computation time step and contain data for (1 ) time 
from ignition (SEC) , (2) mass of propen ant remaining (LBM) , (3) mass 
of inerts remaining (LBM), (4) propellant center of gravity referenced 
to the aft tangent plane ( IN) , (5) propel 1 ant moment of inertia 1n pitch 
or yaw (SLUG-IN**2), and (6) propellant moment of inertia in roll (oLUG- 
IN**2). The execution of this output is controlled by the Output Module. 
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4.3 (Continued) 

^ternal B allistics Module Geometry 4rrav Cards : After the Internal 
Bain sties Module has calculated the end section and reference plane 
geometry arrays and plane constants, these data may be punched on cards 
for use as input in subsequent program runs. , If these data are available 
as mput^in the Internal Ballistics Data Packaqe, then the proqram 
calculations for the data are bypassed and orogram execution time is 
reduced. NPllTAB is the program indicator which controls output for these 
Table 3-XIX for definition). When this ootion is exercised, 
the foil owing data can be punched in NAMELIST format. Definitions for 
these terms can be found in Table 3-XIX. 

REFERENCE PLANE DATA 

NGE0(I) 

GE0C0N(I) 

AP0RT(I) 

TAUPL(I) 

ALPPL(I) 

AKRYP(I)* 

FORE-HEAD AND AFT-HEAD DATA 

NGE0HD NGE0MM XCGNA(I)** 

TAUHD(I) TAUN(I) ‘ XCGNB(I)** 

ABHD(I) ABN(I) PM0rNA(I)** 

PM0IHD(I)* PM0IN(I)* PM0INB(J)** 

RM0IHD(D* RM0IN(D* RM0INA{I)** 

XCGHD(I)* XCGN(D* RM0INB(I)** 

* - These arrays will not be punched if propellant center of 
gravity and moment of inertia calculations are not being 
made (KM0ICG = 1). 

** - These arrays will be punched only when a submerged nozzle 
is simulated (NSUBMG=1). 

4.4 ERROR MESSAGES ^ ^ ^ ^ ^ ^ . 

The majority of the error messages in SR8-II are self-explanatory. The 
following list of error messages require clarification. Messages are 
referenced to the subroutine which generated that message. 
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4.4 (Continued) 


SUBROUTINE EOLBRM ERROR MESSAGES 

1. "DERIVATIVE MATPaK SINGULAR" - This situation occurs very rarely 
inasmuch as, if sinqularities in the matrix sol^Jtions occur, they 
generally occur first in the iteration matrices and the program does 
not get as far as the derivative matrices. However, it is possible 
for the iteration matrix to just barely avoid being singular and for 
the derivative matrix to be singular. When this occurs the above 
message is printed, DLVPT is det equal to -1 , DLVTP is set equal to 
1, and the program continues. 

2. "LOW TEMPERATURE IMPLIES SPECIES SHOULD HAVE BEEN INCLUDED ON AN 
INSERT CARD, RESTART" - this message can occur only for an HP or 
UV problem. It occurs only when the omission of an important 
condensed product of reaction causes the program to seek a combustion 
temperature that is unrealistically Ipv/ (T < lOOK). When this 
occurs the program skips to the next case, 

3. "PHASES OF A CONDENSED SPECIES ARE OUT OF ORDER" - This statement is 
written when the thern:iodynamic 4ata for three or more condensed 
phases of a species are not in a permitted order as discussed in 
Appendix A. After the statement is written, a table of output 

for all completed points is written and program goes to next case. 

4. "SINGULAR MATRIX" - Singularities rarely occur in the. course of 
obtaining a matrix solution in GAUSS. When they do occur the program 
restarts with new estimates. However, if the singularity occurs 
again even with new estimates, the above error message is printed 
and control is returned to the main program Which starts the next 
case. One possible v/ay to get through this difficulty is to assign 

a si ightly modified equivalence ratio or o/f. 

5. "THE TEMPERATURE ~ (Degrees K) IS QUT OF RANGE FOR POINT (Number)" - 
This message is printed whenever the converged temperature for the 
indicated point is outside the temperature range read in on the card 
following the THERMO control card. This temperature range,. which at 
present is 300 to 5000 K, is the one over which the gas phase thermo- 
dynamic data have been fitted. Generally the thermodynamic data can 
be extrapolated a short distance \'/ithout much loss in accuracy. 
However, to prevent large errors due to extrapolation if TLOW/l,5 > 

T > 1.25 THIGH, then after the message _ has been printed, the program 
writes the output for all completed points and then skips to the 
next case i ■ 
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4.4 (Continued) 

6. "35 ITERATIONS DID NOT SATISFY CONVERGENCE REQUIREMENTS FOR THE POINT 

{Number}" - Compositions are typically obtained in 3 to 12 iterations. 
The number 35 was arbitrarily selected to indicate that if convergence 
has not been reached by that number the problem probably will not* 
converge at all. This situation occurs rarely. All completed po.ints 
up to this point are printed and program goes to the next case. If 
the cause of nonconvergence is not obvious from the output, it is 
recommended that the problem be rerun with intermediate output, 

SUBROUTINE ROCKET ERROR MESSAGES 

1. "DID NOT CONVERGE FOR AREA RATIO - (Value of area ratio)" - The 

program permits a maximum of 1.0 iterations to converge to the pressure 
ratio corresponding to the assigned area ratio. The usual number of 
iterations required is 1 to 5., Present experience indicates that the 
only time the number of iterations will exceed is for an assigned 
area ratio very close to 1 such as 1. <Ag/A|^ <1.0001. This is due 
to the fact that the converged throat conditions do not correspond 
exactly to an area ratio of 1. 

SUBROUTINE SEARCH ERROR MESSAGES 

1. "INSUFFICIENT STORAGE FOR FOLLOWING 4Number of) SPECIES" - This 
statement shovys that for the chemical system under consideration, 
the program found more possible species in-THERMO data than can be 
accommodated by storages reserved for the thermodynamic data in 
labeled common/SPECES/, This excess number of species is given in 
the error message. When this situation occurs, the names of the 
possible species are printed and control is returned to the main 
program where the next case is read in. 

This situation can be resolved in two ways. First, the program can 
be recompiled with dimensions increased to accommodate the excess 
species. Secondly, flMIT cards can be used to eliminate the required 
number of excess species. 

The program is currently dimensioned for 115 species which has been 
found adequate for all problems calculated to date. 

SUBROUTINE HCALC ERROR MESSAGES 

1 . "REACTANT (Number) IS NOT IN THERMO DATA" - Enthalpies can be calculated 
by the program only for those reactants that are also included as 
reaction species in the. thermodynamic data. The error message is 
printed when this condition is not met. The program skips to the 
next case. 
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4.4 (Continued) 

2- "REACTANT TEMPERATURE OUT OF RANGE OF THERMO DATA" - Subroutine 
HCALC permits calculation of thermodynamic data for reactants for 
temperatures that extend up to 20 percent beyond the temperature 
range over which the data have been fitted. 

SUBROUTINE LEWIS ERROR MESSAGES 

-1. "ERROR IN LEWIS INPUT CARD. CONTENTS BELOW" - This message indicates 
a keypunching error in a control card, a missing control card, or 
an extraneous card. The entire contents of this card are ignored 
by the program. 

2. "ERROR IN REACTANT CARDS" - This is due to an error in a chemical 
symbol s ch as the symbol not being left-adjusted or not included 
in BL0CK DATA. Program skips to next case, 

SUBROUTINE ENDCSB ERROR MESSAGES 

1 . "ALPHA E MAX TOO SMALL" 

A0EM is input negative or zero. 

2. "REl TOO LARGE" 

Nozzle opening diameter is greater than case diameter. 

3. "INVALID BETA E" 

Ellipse ratio must be 1.0 + .01. . 

SUBROUTINE SCI ERROR MESSAGES 

1. "TDMAX GREATER THAN TAUM" Maximum burn distance in fore-head 
is greater than that for forward tangent plant. 

SUBROUTINE GAMA2S ERROR MESSAGES 

1. "NEGATIVE HEAD END WEB FRACTION Increase size of BTA0E. 

SUBROUTINE AIBST ERROR MESSAGES 

1. "ICOUNT - X^ ITERATION FOR INCREMENT DISCHARGE PRESSURE. 

This comment usually indicates a data error; however, it 
is possible for convergence failure in particular cases 
where the port- to-throat ratio os less than 1,10. 

2. "ITERATION FOR INCREMENT X X DIGCHARGE TEMPERATURE..,." 

See (1) abovei Also, frequently indicates a c* input value 
which is not sufficiently consistent with burn rate equation 
constants, or an incorrect input value f or VFH0 or VCHINP. 
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4.4 (Continued) 

SUBROUTINE SETPH ERROR MESSAGES 

1. "AT = xxxx.xxx, THROAT AREA GREATER....." The throat area 

is too large to maintain sonic flow; decrease throat diameter. 

This comment is fallov/ed by the relevant nozzle and discharge 
parameters. 

2. "ITERATION FOR NOZZLE PRESSURE RATIO....." Failure at 
this point generally indicates an input error. 

SUBROUTINE TISUB ERROR MESSAGES 

1. "SLOT NUMBER HAS BURNED INTO ...." With anisotropic 
burning it is possible for the slot face to burn into 
the fore-head section; relocate slot or modify burn rate. 

SUBROUTINE LIN ERROR MESSAGES 

1. "THE INDEPENDENT ENTRY IS ABOVE OR BELOW THE TABLE" "TABLE NO. = N" - 
The dependent variable is out of the range of XTABLE, The subroutine 
LIN is called from a subroutine with the table number N in the 
call statement. The program will stop when this error occurs. 

SUBROUTINE XLIN ERROR MESSAGES 

1 .. "X IS OUT OF RANGE IN XLIN" - This error statement may be generated 
only from calls in the IBM module. The program will continue after 
this message. The independent variable that is returned from XLIN 
after this message is issued will be unpredictable. • 

SUBROUTINE IDNjaZL ERROR MESSAGES 

1. "THE VALUE COMPUTED FOR GAMMA IS NOT WITHIN THE LIMITS (1.0, 1,7). 
CASE REJECTED" - Check- input data for ID NOZZLE control section. 
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APPENDIX A 


THERMO DATA {FORMAT AND LISTING) 
MSFC TAPE 2381 0 


ORDER 

CONTENTS 

FORMAT 

CARD COLUMNS 

1 ■ 

THERMO* 

3A4 

1 to 

6 

2 

Temperature ranges for 2 sets of 
coefficients: lowest T, common T, 
and highest- T* 

3F10.3 

1 to 

30 

3 

Species name 

3A4 

1 to 

12 


Date 

2A3 

19 to 

24 


Atomic symbols and formula 

4(A2,F3.0) 

25 to 

44 


Phase of species (S, L, or G for 
solid, liquid, or gas, respectively) 

A1 

45 



Temperature range 

2F10.3 

46 to 

65 


I nteger 1 

115 

80 


4 

Coefficients a^-(i = 1 to 5 ) in equa- 
tions (90) to (92) (for upper temp- 
erature interval) 

5(E15.8) 

1 to 

. 

75 


Integer 2 

15 

80 


5 

Coefficients in equations (90) to 
(92) (ag, 37 for upper temperature 
intervaT, and a-i , an, and an for 
lower) 

5(E15.8) 

T to 

60 


Integer 3 

15 

80 


6 

Coefficients in equations (90) to (92) 
(sij ag, ag, ay for lov/er temperature 
intervaT) 

4(ET5.8) 




Integer 4 ^ 

120 

80 


(a) 

Repeat cards numbered 1 to 4 in cc 80 
for each species 




(Final 

card) 

END (Indicates end of thermodynamic 
data) 

, 3A4 

1 to 

3 


3 

Gaseous species and condensed species with only one condensed phase 
can be in any order. However, the sets for two or more condensed 
phases of the same specTes must be adjacent. If there are more than 
two condensed phases of a species, their sets must be either in in- 
creasing or decreasing order according to their temperature intervals, 

*M0T INCLUDED ON THERMO TAPE. 
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300.000 100 

ALtS) 

0 . 

0 . 

-O.AA.923'J93E-Oa 

AL(L) 

0.381R5052F 0! 
-0.961l6HA4f 02 

U. 

AL 

0.2545n6SOE 01 
0.3Hi*9P9S7F 05 
-0. l54a77f.9F-08 
AL* 

0.?813R51ftF 01 
0.10BSq616F 06 
-0.97‘}lR340E-ll 
ALBn2 

D.7172P995F 01 
-0.676B36B2K 05 
0,1025132<fF-07 
ALCL 

0.A3963254E 01 
-0.7S3^264flE 04 
0.42909I63E-08' 
ALCL* 

O.A3784Sfl5E*01 

D.101H1202E+06 

0.4S906177F-0H' 
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E2C3 J12/65 

0.S5015760Ct0l 
-C. 10323555E+06 
0.1l‘l27358e-07 


P 2. 0 3. 

0.573920036-02 
-o.isessicoMoa 
-0.2561991SF-1 1 


no n. 00 0 . g 
-0. 195238151-05 
G. 367BC531HC.I 
-1. 102030075*06 


3q 0.000 

0.36115921C-0O 
0.201 38 3861: -0 1 
0.O0209 6068*01 


5000.000 

-n. 258552 llE-13 

-0, 21637125E-05 


B3C3CL3 J 3/65 

0. 19?82565E*C2 
-0.2C320GB3E *06 
0*222231286-07 


K 3. II 3. 
0.631725816-02 
■0, 678983326*02 
-0.151812956-11 


r.L 3. 00 0. G 

-0,272529266-05 
0.505559836*01 
-0, 199516326*06 


3C0.000 

0.520579106-09 

0.552605976-01 

0.905355276*01 


5000.000 

-O. 366777906-13 

-0.557507616-05 


63C363 J 3/65 

0.168586166*02 
--0.290931056*06 
0. 255380396-08 


8 3. n 3. 

0. £86857556-02 
-0.598719586*02 
0.553587616-11 


6 3.’ 00 0. G 

-0.373810586-05 
0.307988616+01 
-0.207122136+06 


300.000 

0.71 £705016-09 
0.556365926-01 
0.U5621966+0H 


5000,000 

-0.503769176-13 

-0.330980266-05 


BEISI J 9/61 

0.163862576*01 
-8,555308760,03 
0. 139335576-07 


KE 1 . n 0 . 
0. 139810706-02 
■0.871561656+01 
•0.351155536-11 


00. 0 0. S 

-n.S191b095£-08 
-0.53E86L976+00 
-0. 278182370+03 


300.000 

-O. 503575516-09 

0.132509586-01 

0.103511006+01 


1556.000 
Q. 155765526-12 
-0.190602 676-05 


E61U J 9/61 

0. 315599976+01 
C.30556750E+03 

C. 


R6 1. 0 0. 

-0.565733566-05 
-0.165565356+02 

0 . 


no. 0 0. t. 
0.277335866-06 
0 . 

0. 


1556.000 

-0.895316776-10 

0 . 

0 . 


5000.000 
0. 01533L236-L5 

0 . 


BE J 9/61 

0,250603506+01 
0.386661556+05 
0.355653296-10 


HE 1. 00 0. 

0,138727516-03 
Q.265555C26+01 
-0.135381076-13 


00 0. 00 0. C 

-0, 119132066-06 
0.259825606+01 
0.386331506+0‘J 


300.000 

0.259337276-10 

0.129100086-05 

0.215011156+01 


5000.800 

-0,565236506-15 

-0.331323836-07 


36+ J 6/65 

8.251013626+01 
0.157555766+06 
0,110802316-09 


06 1. 6 -1. 
■0.216B37506-05 
0.277075906+01 
■0.517393986-13 


00 0 . 00 0 . c 

0, 157566856-07 
0.2595273HE+01 
0, 157559856*06 


300.000 

-0.577082326-11 

0.513921016-05 

0 . 285986916*01 


5000.000 
0.525655966-15 
-0. 105 T2B F76-06 


BE0C2 J 6/66 

8.691033766*01 
-0.605057156+05 
0.60B65373E-08 


36 1. n 1. 
0.326636356-02 
-□.917585816+01 
“0.172762856-12 


n 2. on 0. G 
-0,136781206-05 
0.200691206+01 
-0.592351976+05 


300.000 

0,25 5762116.-09 
O'. 180553256-01 
0. 157923656 + 0? 


5800.000 

+0* 177Z7751E+13 

-0.1691 75316-85 


6ECL J 9/66 

0.510S2B7R6+01 
0,597530603+05 
0.158525076-08 


HE 1, r,L 1. 
0.5T5617C16-03 
0.255136556+01 
0.552068956-15 


no 0. 80 0, C 

-0.1799652UE-IJ6 
0.283219876101 
0.7i29062586»P5 


300.000 

0.32563903F+10 
0.555667656-02 
0.118 711513751*01 


5000,000 
-0.206528501-15 
-0,55502 1616-05 


BECl* J 6/68 

0.538275006*01 
Q, I 15997176106 
0.356326506-08 


86 1 . CL 1 . 
-0.105711906-02 
-0.507537596101 
-0.659250806-12 


6 - I . 0 p . C 

8.111236036-05 
0.209659851 *01 
0.116715661*06 


300.000 

-0, 1695299&3-09 
0, ‘.12675921-02 
n. 711251 3591 nil 


5000.000 

n.6l007C’!lF-15 

-0.655274111-05 


A-7 
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EFCLF j 6/65 

I'F 1 . Ct. 

1. 

F t, 00 ij. 

h 

K'O.OCO 

0.6A402791E+01 

0,114636938-02 


-0.4H545 I6OF-06 


1). <1287;i65t -If) 

-0.71O59771EtO5 

-0, I741h5S6E»C-1 


0 , 4102a 18 ir+nl 


n? 

0.AOO76232E-OU 

-0.5If27939E-l2 


-0. /i)46i»71t>F + 07 


n.'^iui(600A'»i ►ni 

EECt.2(5) j 6/65 

R8 1 . CU 

2. 

00 0. 00 11. 

S 

iOO^OGO 

Q • 

0. 


0. 


0. 

0* 

0« 


0.300657458+01 


0. 1‘#*J . 

-0.29604l5flE-07 

0.23‘:34861F-IC 


-O.60 7221C0F+05 


-0. 12579 772f +02 

2ECL2(U J 6/55 

lie 1. CL 

2. 

00 0. 00 0. 

L 


O.U603719F.+02 

0. 


0. 


0. 

-Q.54698417Et05 

-0.76448784EF02 


0.146037198+02 


0. 

0 ■ 

0. 


-0.64498417F+05 


-0. 76445 7H4F+02 

BeCU J 6/65 

E8 1. CL 

2. 

00 0. 00 0. 

G 

3(*o*noo 

0.6 7043191E+01 

0.H7166468E-03 


-0. 372S5053E-06 


0. 705671)061-10 

-0.5549455BEt-05 

-0.H43516418VO1 


0,44927125E+ai 


0.<J05 155451-02 

0.5i;a9 ;o49E-o8 

-0.534980928-12 


-0.449528H1F+C5 



REF J J/63 

HS 1 , F 

1. 

CO 0* 00 0* 

G 

jon^ooo 

0.3/lt.5220F*-0l 

0.EB446143F-03 


-0.3571 36l2F-0ti 


0,t*700/4Ff>tJ -IC 

-0.2/)229OS2E+O5 

0.3123976IJ8F01 


0.12 725419F.+01 


n.28il4-l431F-03 

-0.522892548-08 

0.196293208-11 


-0.260126856+05 


0.586737168+01 

FEF2tS) J 6/70 

RE 1. F 

2. 

n 0 . 0 <v. 

S 

3*:o.noc 

n* 

0. 


0. 



c « . 

0. 


0.20593770F+02 


-n.Ati 

a.89a00555n-06 

-0.96669264E-09 


-0.1269370BF+C6 


-0.91 785 1 1 3F+02 

EEF2IS1 J 6/70 

RE 1. F 

*2. 

00. 0 C. 

S 

•^no^onQ 

0* 

0. 


0. 


0. 

Q« 

0, 


0.569655768 +01 


O.4O25+358F-02 

0« 

0. 


-0.1252HUH4F+06 


-n./» ?r*n TH«?Hfn2 

BEF2(U J 6/70 

B8 1. F 

2. 

on. n 0 . 

L 

H/*>.noo 

0.609B9639E+01 

0.43320498E-C2 


0.1875440JE-06 


~C. 1601 7482 L-O') 

-0.l25U36lEtO6 

-0.29C262488>n2 


0.774233618+01 


-0.696800651-0 1 

0.31262592E-08 

-0.254562798-11 


-0.125465306+06 


-0. 17440 3046+02 

EEF2 J 6/70 

PE 1. F 

2. 

• 0 0. 0 0. 

G 

IQOpOOO 

0.6045763ie<-01 

0. I5629374E-02 


-0.6610B197F-06 


0. 1244/551E-09 

-C.97779127E+05 

-O.793104556tOl 


0.352342748+01 


0.9 1i1902]14F-ii7 

• a.429H09a9E-08 

-0.57 t51ll3E-l2 


-0.971304618+05 


n.Ad7nHi25F+oi 

PEN J 3/63 

8E 1 . H 

1. 

00 0. 00 0. 

G 

300.000 

0.30570218£*01 

0.14977223E-02 


■ -0.568 729638-06 


0. 10260517F-09 

C.37639513F+05 

0.33‘!71052Et01 


0.371123058+01 


-0.191415481-02 

-0.329258835-08 

0.66638562E-12 


0.375655608+05 


0.17541895F+00 

FEH+ J 9/66 

HE 1. H 

1. 

8-1. 00 0. 

G 

300.000 

Q.2'J0l599 2Ei-01 

0.16751761E-02 


-0.66805503E-06 


0. 1251095 lE-09 

0.1 3816812E+06 

0.3542S41SetOl 


0.370957128*01 


-0.15H52031F-02 

-0.189332218-08 

0.I7173264E-12 


0.133028668+06 


-0.2y*>9«407e>no 

HEN J 6/63 

HE 1- H 

1. 

00 0. 00 0. 

c 

IGO.OOO 

0.37BSS93VEt01 

0.823865758-03 


-0.3"U6n2E-0& 


0.61 'JSIOHHF-IO 

0.50066180Ft05 

0.309241948+01 


0,316842868+01 


Q, l02H2<.a3r-02 

-0.434810996-08 

0.17534453E-11 


Q.50310451F+0S 


0,6649359lf-f0l 

EECtS) J 9/63 

88 1 . n 

1. 

00 0. 00 0. 

S 

300.000 

Q.3 31aO4SOE+01 

0. 559172128-02 


-0.463012718-05 


0,191QH779C“0H 

-0.7345 15536F05 

-0. 19 304223E+02 


-0.197233488+01 



0. 1954003 lE-07 

-0.492163198-11 


-0.722700698+05 


0.6711631196+Ot 

EECIU J 9/63 

EE 1. 0 

1- 

00 D. 00 0. 

L 

2820.000 

0.a05l2e67Ef01 

0. 


0. 


0. 

-O.71290168E+O5 

-0.46972j466+02 


Q+6 


0. 

0. 

0. 


0. - . 


Q* 

8E0 J 9/63 

OE 1. 0 

1. 

no 0 . 00 0 . 

G 

300.000 

0,3504a549EF0L 

Q.t0967BB4B-02 


-0.439908778-06 


Opa20H8<)67F-lO 

0, 144350 lOE+OS 

0.331673478+01 


0,352506208+01 


-0.130442818-02 

-0.63674976E-08 

0.20911849E-11 


0.145621928+05 


0.38 1214568+01 

EECH J 9/63 

88 1. 0 

1. 

HI. 00 0. 

G 

300.000 

0.36351960E+01 

0.28532436E-02 


-0.102303268-05 


0.1706'19)4F-09 

. -0.13B56067F*05 

0. 512962948+01 


0. 37;»703'j66+0L 


-n. 34554324F-01 

-0,922690468-08 

0.325183028-11 


-0.13732272F+05 


0.53H99404Ff0l 


‘ lOT.ocn 

J'.'.'lt -l'( 
■i!. -ri VI yif It -O') 


000 

n. 

-n.f,ftOI W.C'iF-O') 


•)OOfi.''CO 

0 . 

0 . 


•inoii.'JOu 
-r».<V'M 1'>»*0F-U 
-n.art II I? I'lF-C) 


5000. not) 

-P.t.6t.?ld30E-lti 

n.‘.oit)t3':tF-05 


500.000 

0. 

-0 . 12002 lone-oi 


025,000 

n. 

0 . 


2000.000 
0.nU3iJH22E-l3 
0.?57‘) JACAE-05 


5 COO. 000 

-n.fl57ttiOA3F-15 

-0.n5t)3f)2CBE-05 


5000.000 
-0.59lft6'3 70S-l4 
0.4691O325F-O5 


5000.000 
-o.m 7t.lt,66E-l4 
n.36?2i)76Oe-0S 


SOOOiOOO 
-i).52i<OO051E-15 
G.27 176C17E-05 


2020.000 
-<'.27HA30 57E-12 
-0. J1A9HJ 3OF -04 


5000.000 

0 . 

0. 


5000.000 

-0.5677HHC0E-14 

n.62H3S237E-05 


5000.000 

-0.l0iK0'7’JAg-l3 

0.fi21fl8520e-05 
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PECt*+ 

J 6/6U 

l!6 1 . n 

1. 

Hi. f -1. 

0 

300.0 GO 

9COO.OCC 

1 


0.375<iqn5f + 01 

0.339395H3E-02 


-0.125062638-06 


0.2159!.*l27r-O9 

-0.1511* 7 108-13 

f. 


C.q6l<966C<)l *196 

0. I905 70995 *01 


0.2961'?2?Cf *01 


0, 1697?'Jil9[-02 

0.99557i661 + -05 

n 


-0.l03R2')32E-07 

0.933877338-11 


0.966937535*05 


0, <183*1 5035*01 

4 

EEC2H2 j 3/6? 

PE 1. I) 

2. 

H 2, 00 0. 

0 

300.000 

5<*00.000 

L 


0.6<.H63231F*OL 

0.583883391-02 


-0.218392038-05 


0,3766'j29nr-<?9 

-0.255609 15fc-13 

7 


-0.ail2'.?3RE*fj6 

-0,936801 19F *01 


0.?96Mi0?2* *01 


O.UM25l*8 'OI 

-0.165555 17,8-05 

3 


O.62732150E-OH 

-0.306975378-12 


-o.iini37n«r*o6 


0. 101155*171 *C? 


H 

EE2C 

J 9/63 

P6 ?. 0 

1. 

CO 0. 00 c. 

fi 

30n.0<30 

5001 . 000 

1 


0.56<i6<)r36C*lJl 

0,219703858-02 


-0.929195788-,?'. 


0. 1 /5965l0f -09 

-0. 12189') 52/ -13 

7 


-C.9A95R9fj'iPtO<t 

-0.56R358£OE*0l 


n.2752 7B975*< , 


0.ii9/i*li679f -f./ 

-n.')5iirri2 575-09 

3 


-0i367691 8fiE-09 

0.11C15972E-11 


-0.871797095*09 


0.653875565*01 



l'E2f.F2 ,1 6/66 

lie 2. 0 

1. 

1 2. 00 0. 

G 

360.0C0 

5009.000 

1 


0. 10311 3<i3E*02 

0.29 2581 618-02 


-0-129819875-05 


0.2365?l69r-0<) 

-0, 1655<)lf08-13 

V 


-0.l<>H9662 3f*06 

-0.295008965*02 


0.98600026.5 *01 


n, 199399B2*-01 

-0. lH81Fj7fcOf -05 

S 


0.7l009‘i03f-OH 

-0.372252588-12 


-0. 19 7039598*06 


0. 12285 70*5*01 


<r 

EE2r2 

J 9/63 

K8 2 . n 

2. 

00 0. 00 0. 

G 

300. OCD 

5000.000 

1 


0.7!7B3652E*0l 

0.3a7<36926f-C2 


-0.131622735-05 


0.259735155-09 

-0. I 7597,1 390-13 

A 


-0.5 19B4B76E*05 

-0.125387588*02 


0. 171027395*131 


0-182559395-01 

-0.153772538-05 

3 


0.21.260in6E-08 

0.196919938-11 


-0.505123665*05 


0.15201 3571*02 


A 

lit JCl 

J 9/63 

PE 3. 0 

3. 

00 0. 00 C. 

C 

300, OOG 

■.nor . 000 

I 


0.91907322E*01 

0.736237C1E-C2 


-0.312927295-05 


0. 591625891 -09 

-0.513601958-13 

? 


-n.l306lB99E*06 

-0.233300938*02 


0.200026925*01 


0. 2000‘j 1 725,01 

0.575178575-06 

3 


-0. 17092B05E-07 

0.898627858-11 


-0. 128268678*06 


0.156077905*0? 


A 

RE6CA 

J 9/63 

88 9. n 

9, 

CO 0. 00 C. 

C 

300.000 

5000.000 

1 


C,195A7030E*02 

0.819037305-02 


-0.351627898-05 


0- 66923*5 7! -C9 

-0.570059538-13 

2 


-0. l970<iS95P*06 

-0.515099298*02 


-0.138185385. *01 


0.523n*H2RF-<31 

-.9.508930 185-05 

i 


D.<i7379707E-0a 

0.999591698-11 


-0. 192783568*06 


0.303999038*02 


A 

pH 

J 9/61 

HR 1. 0 

0. 

0 0. 0 0. 

G 

300.000 

5000.000 

1 


0.210781095*01 

0.673570075-03 


-0.295059318-96 


0.35C9*37*F-10 

-0. 175591965-15 

?. 


0.128528195*06 

0.895518095*01 


0.296595658*01 


0.2901160*5-03 

-0.91H59575E-06 

i 


0, in07991£-OB 

-0,906696075-12 


0,127139598*05 


0.692937876*01 


A 

BR2ILJ J <9/61 

KR 2. 0 

0. 

0 0. 0 0. 

L 

265.900 

1900,000 

1 


0. 

0. 


0. 


0. 

0. 

2 


0. 

0. 


0.155019198*02 


-0. *58068268-01 

0. 119769865-03 

3 


-0.12886307P-06 

0.999909108-10 


-0.336675975*09 


-0.603336068*02 


A 

HR2 

J12/61 

GR 2. 0 

0. 

o 

o 

o 

o 

G 

300,000 

57)00.000 

1 


0.<i9979995E + Ql 

0.10C512C88-03 


-0.163938168-07 


0.226856218-1 1 

-0. 102367 758 - 15 

? 


0.236599915*09 

0.90e8S93ie*<31 


0. 389695805*01 


0.261U89IC-02 

-0.509351578-05 

3 


0.281206895-08 

-0.732562C26-12 


0.298969895+09 


0.6969698 55+01 


A 

CIS) 

J 3/61 

C 1. 0 

Q. 

0 0. 0 0, 

S 

. 300.000 

5000.000 

1 


0. 13609992E *01 

0. 191822378-02 


-0.890903898-06 


0.165*87076-09 

-0. 117,726 798-13 

2 


-0,6‘>7!387HE*03 

-0.800702078*01 


-0.997780535*00 


0.531.910028-02 

-0. 397755716-06 

3 


-0.909592985-08 

0.21 1399 395-11 


-0.992806R8E+02 


0.168*07918*01 


A 

C 

.1 3/61 

Cl. 00 

0. 

00 0. 00 0. 

C 

300.000 

5000.000 

1 


0.25B10663F*01 

-0.196962025-03 


0.793880695-07 


“C. 795810798-11 

0.589009778-16 

2 


0.852162995*06 

0.931288795*01 


0.253207056*01 


-0. 1911876518-0 3 

0,306820626-06 

3 


-C.26T70Q69E-09 

0.879888278-13 


0.85290*225*05 


0.560623755*01 


A 

C + 

L12/66 

C 1. 5 - 

•1, 

00 0. OO 0^ 

G 

300.000 

5000. COO 

1 


0.251182798*01 

-0.172597898-09 


0.950926765-08 


-3.221885188-11 

0. 186218928-15 

2 


0.216677218*06 

0.928612988+01 


0.259530505*01 


-0.506866555-03 

0.689236698-06 

3 


-G.52669B7BE-09 

0.15C83377E-1Z 


0.216662B1P+06 


0.389572988501 


A 

c- 

J 9/65 

Cl. 6 

1.' 

00 0, OO 0. 

G 

300.000 

5000.000 

1 


0.299705915*01 

0,112869288-03 


-0,78591562E-37 


0.197786155-10 

-0. 111055558-15 

2 


0.699729695*05 

0.923569925*01 


0.259256506*01 


0.931530688-n* 

-0.133079951 -06 

3 


0. 13951 3798-09 

-0.521509928-13 


0.699557578*05 


0. 398.116578 + 01 


A 

CCL 

J12/69 

C 1 . CL 

1. 

no. no. 

G 

300.000 

5000. COO 

L 


0.939B9 ?278*0i 

0. 500778958-03 


-0.200128336-06 


0,386809928-10 

-0,255511 131- -15 

2 


0,590765995*06 

0.333702015*01 


0.319525575*01 


0,2ii077.31 Bl-02 
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CURVE FITTING EQUATIONS FOR DERIVING SRB-II INPUT DATA 


A substantial savings in computer time may be realized by using the 
program options to input curve fitted data for output from the analytical 
simulation models rather than executing the programs themselves. Ex- 
amples of these input data are CSC0EF(1), (2), and (3) for the 
theoretical c* simulation derived from the LEWIS Subroutine; and AK(1), 
(2), and (3) for theoretical vacuum Isp versus expansion ratio at 1000 
psia derived from the LEWIS Subroutine. The c* curve fit coefficients 
are derived automatically in the Isp modules but the vacuum Isp curve 
fits are not. Thus, it is necessary for the user to exercise a curve 
fitting procedure in order to obtain the necessary coefficients for the 
input data. 

Least Squares Second Order Fit 

If we choose to represent a set of data (x^, y^}, where i = l,2,3,,..,n, 
by the same relationship y = f(x), containing r unknov^n parameters a-i, 
a 2 ,a 3 , . . . ,a^ , and form the deviations as 

Vi = f(xi) - yi 


Then the sum of the squares of the deviations 
n . . n „ 

S = s Vi = s Ci^(xi) - yi]’ 
i=l i=l 


is clearly a function of a-j,a 2 ,a 3 ,...,a^. We can then determine the a's 
to that S is a minimum. 

Now if S(a],a 2 ,a 3 ,...,aj.) is a minimum, then at a point in question 


8S 

aa-j 


= 0 , 


— = 0 = 0 - g 

aa2 ’ 333 aay. 


0 ) 


The set of r equations, (1), called normal equations, serve to determine 
the r unknown a's in y = f(x). This particular "best-fit" criterion is 
known as the principle of least squares. 

Now for a second degree polynomial fit. 

y = a-| + 32X + a 3 X^ 


and 

v^. = a-] + a 2 X^- + 33X^.2 - y. 


R-1 



D256-1 0020-2 
APPENDIX B (CONTINUED) 


then 


3v. 3v. av. 

1= 1 __L= X. ^ = X 2 
I ’sag 333 ^i 


3a 


Since the normal equations for a polynomial take the form 
2 V . ^^i 

i=l SaT " k = 1 ,2,3 


The normal equations are 


n 


£ (a, -f a^x. + 


i=l 

n 


1 “2^i 


a^x 


- y,.) -1 = 0 


.2^ (a^ + agX. + a3x/ - y^) x. = 0 

“ •i'i) =‘0 

Collecting the coefficients of a. 

J 

na, = ^E^y. 

= ^,.y, 

Solving (2), (3) and (4) for a., 

a = Uxzy-nzxy)(£xzx^-nzx^)-(Ex^i:v-ni:x^)r(zx)^-nsx^l 
(sxLx^-nEx^)2-C(Ex^)2-nzx'^][{Ex)^-nEx2] 
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32 = 


Uxzy - nzxy) - (zxzx^ - nzx^) 


C(zx)^ - nrx^] 


( 6 ) 


"1 = 


ly - zx ag - EX 83 


(7) 


NOTE: A11_e in equations (5), ( 6 ), and (7)- are from i=l to n and 
variables x and y are subscripted^.. 

Equati ons (5), ( 6 ) and (7) may be used to derive the reouired 
curve fit coefficients. 7 ’ 


require a semi-log to the base e curve fit (AK(10, (11), 
value of log (y.) should be substituted in equation 
(5), ( 6 ) and (7) for y.. ^ 
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